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FORBVORD 


This  report  v;as  prepared  by  the  Design  Criteria  Branch,  Materials  Laboratrry, 
Directorate  of  Laboratories,  Wright  Air  Development  Division.  Preparation 
of  the  report  was  initiated  under  Project  No.73Sl»  “Platerials  Application," 

Task  No.  73SI2,  "Data  Collection  and  Correlation,"  with  Sidney  Allinikov  act¬ 
ing  as  project  engineer. 

The  compilation  of  papers  presem;ed  herein  represent  a  significant  effort 
in  the  development  of  design  data  applicable  to  cylindrical  sandwich  constructions. 
The  purpose  of  this  report  is  to  provide  wide  distribution  of  these  particular 
papers,  which  has  not  been  done  previously. 
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AEbTRACT 


This  report  is  a  compilation  of  papers  which  present  a  comprehensive 
treatment  of  the  theories  and  parameters  associated  with  the  design  of 
cylindrical  sandwich  constructions.  ^lany  of  the  formulas  developed  are 
applicable  to  a  wide  variety  of  core  and  facing  combinations.  Experimental 
data  on  flat  and  curved  sandwich  sections  are  furnished  to  support  the 
theoretical  solutions  related  to  the  design  of  these  structures. 
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ELASTIC  STABILITY  OF  CYLINDRICAL  SANDWICH 
SHELLS  UNDER  AXIAL  AND  LATERAL  LOAdI 


By 

EVERETT  EUGENE  HAFT,  Engineer 

Forest  Products  Laboratory,  —  Forest  Service 
U.  S.  Department  of  Agriculture 


Summary 

A  linear  solution  for  the  determination  of  the  loads  under  which 
a  cylindrical  sandwich  shell  will  buckle  is  presented.  The  facings  of 
the  sandwich  cylinder  are  treated  as  cylindrical  shells  and  the  core  as 
an  orthotropic  elastic  body.  The  method  of  solution  is  of  interest  in  that 
it  is  of  sufficient  generality  to  be  applied  to  many  problems  in  sandwich 
analysis.  The  characteristic  determinant  that  represents  the  solution  to 
the  problem  is  solved  numerically.  Curves  that  show  how  the  buckling 
load  changes  as  the  parameters  of  the  problem  change  are  given. 

—This  report  is  one  of  a  series  prepared  by  the  Forest  Products  Laboratory 
under  U.  S.  Navy,  Bureau  of  Aeronautics  Order  No.  01593. 

—Maintained  at  Madison,  Wis.  ,  in  cooperation  with  the  University  of 
3Vis  consin. 

Manuscript  released  25  February  X9oC  for  publication  as  a  VkDu  Technical 
Hepcrt. 


Introduction 


Sandwich  construction  is  a  resxilt  of  the  search  for  a  strong,  stiff, 
and  yet  light  weight  material.  It  is  usually  made  by  gluing  relatively 
thin  sheets  of  a  strong  material  to  the  faces  of  relatively  thick  but 
light  weight,  and  often  weak,  material.  The  outer  sheets  are  called 
"facings"  and  the  inner  layer  is  called  the  "core." 

Such  a  layered  system  presents  difficult  design  problems.  What 
is  offered  here  is  a  straightforward  method  for  dealing  with  some  of 
these  problems. 

The  problem  to  which  the  method  ’  ipplied  is  that  of  the  elastic 
stability  of  a  sandwich  cylinder  under  uniform  external  lateral  load 
and  uniform  axial  load. 
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Notation 

- 

CD 

radial,  tangential,  and  longitudinal  coordinates. 

respectively 

*  a 

radius  to  middle  surface  of  outer  facing 

b 

radius  to  middle  surface  of  inner  facing 

t 

thickness  of  each  facing 

i 

length  of  cylinder 

£ 

modulus  of  elasticity  of  facings 

Poisson's  ratio  of  facings 

G 

modxilus  of  rigidity  of  facings 

£ 

modulus  of  elasticity  of  core  in  direction  normal  to 

facings 

^re 

modulus  of  rigidity  of  core  in  r9  plane 

Grz 

modulus  of  rigidity  of  core  in  rz  plane 

q 

intensity  of  uniform  external  lateral  loading 

k 

1 

£t  log  b 
b  = 

1  +  -  ^ 

^  E  a 

c 

<r 

r 

normal  stress  in  core  in  radial  direction 

r9  rz 

transverse  shear  stresses  in  core 

U,  V,  w 

radial,  tangential,  and  longitudinal  displacements. 

y 

respectively 

n 

number  of  waves  in  circumference  of  buckled 
cylinder 

.  H- 

H 

* 

m 

number  of  half  waves  in  length  of  buckled  cylinder 

s/aLL  TR  bO-133 

j 

X 

mira 

i 

fine 

c  -  n 

2G  -  2 

r0 

fiz 

/ 

^  i. 

N„,  N  ,  N 

0  z  6z 

normal  forces  and  shear  force  per  unit  length  of 
facing  .  '  ■  .  ' 

D 

CD 

D 

N 

transverse  shear  forces  per  unit  length  of  facing 

Mg. 

bending  moments  per  unit  length  of  facing 

^20’  ^0Z 

twisting  moments  per  unit  length  of  facing 

R,  0,  Z 

surface  forces  per  unit  area  of  facing 

P 

E^a  (1  - 

Et 

qa  (1  - 
Et 

‘*’2 

(1  - 

Et 

a 

t2 

12a^ 

a* 

t^ 

12b^ 

log 

natural  logarithm 

A.  B.  C,  D,  K,  U  A' 

arbitrary  constants 

note  --  any  of  the  above  terms  that  appear  with  a  prime  (as  N .' ) 
refer  to  the  inner  facing. 
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Mathematical  Analysis 


As  previously  stated,  the  core  is  relatively  weak.  Because  of 
the  high  strength  of  the  facings  the  core  need  carry  little  tension  or 
compression  except  in  a  direction  perpendicular  to  the  facings.  The 
facings  are  able  to  resist  shearing  deformation  in  their  plane  and  it 
is  necessary  only  that  the  core  be  able  to  resist  shear  in  the  radial 
direction  in  planes  perpendicular  to  the  facings.--  In  this  analysis  the 
core  is  considered  to  be  an  orthotropic  elastic  body.  It  is  unable  to 
resist  deformations  other  than  those  just  mentioned.  This  assumption 
makes  it  possible  to  determine  explicitly  how  the  stresses  vary  through¬ 
out  the  thickness  of  the  core. 

The  facings  arc  treated  as  shells. 

^  '  Intcr.depcndn  nee  of  the  c  ore  and  the  facings  is  gained  by  equating 

the  ir,disipla'c  erne  Ills  .at  the  interfaces.  To  simplify  the  analysis  the 

'  ‘’Of  •'  - 

core  is  a.ssunictl  to  extend  to  the  middle  surface  of  each  facing. 

Figure  1  ‘show-s  the  cylinder  and  the  coordinates  that  are  used. 

Prebuckling  Stresses 

Before  buckling  occurs  the  cylinder  is  in  a  state  of  uniform 
compression.  The  axial  load  is  carried  by  tlie  facings  since  the  core 
material  is  assumed  to  be  incapable  of  carrying  load  in  this  direction. 
Vfith  facings  of  like  material  the  stress  i.'s  the  same  in  both  facings. 
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If,  in  addition,  the  facings  have  the  same  thickness,  then  the  loading 
per  unit  length  of  facing,  Njj  or  N^' ,  will  be  the  same.  This  means 
that  for  a  total  load  P, 

2ira  N,  +  2irb  N,'  =  P. 

z  z 

The  calculation  of  stresses  due  to  the  lateral  pressure  is  a 
problem  in  rotational  symmetry.  Differential  elements  of  the  core  and 
of  the  facings  are  shown  in  figure  2. 

Summing  forces  in  the  radial  direction  gives  for  the  core 

^  +  =  0. 

Sr  *■ 


N„  -  Et  (  +  ^  )  .and 
y  r  =  a 

N  '  =  Et  (  +  ^)  ,  these  equations  can  be  solved  for  a  ,  N 

0  b  r  s=  b 

and  Ng' .  The  results  are* 

cr  =  q  ^  k 
r  r 

N.  =  qa  (1  -  k),  and 

y 

*For  a  more  detailed  derivation  of  these  terms  see  Reference  1. 
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N  '  =qak  where 

e 


b  £t  log  ^ 

1  + - ± 

a  £  a 
c 

As  P  and  q  increase,  N  ,  N  N„,  N„' ,  and  a  also  increase. 

z  z  8  0  r 

£ventually  a  condition  may  be  reached  where  a  slight  increase  in  load 
causes  the  cylinder  to  lose  its  state  of  uniform  compression  and  buckle 
as  a  result  of  elastic  instability.  This  buckling  is  assumed  to  cause 
only  a  small  change  in  the  stress  distribution.  These  small  changes 
will  now  be  considered. 


Buckling  Stresses 
The  Core 

A  free  body  diagram  of  an  element  of  the  core  is  shown  in 
figure  3. 

Neglecting  terms  which  are  products  of  more  than  three  differ¬ 
entials,  a  summation  of  forces  in  the  r,  0,  and  z  direction  gives 


cr  +  r 


^'^r  ,  9'*’  r0 


30 


+  r 


rz 


3^ 


=  0 


ra 


3r 


+  2 1 


re 


=  0 


(1) 

(2) 


T  +  r  ^’’’rz  =  0 

rz  - 

3r 


(3) 
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and 


(20) 


u  =  (Aa  +  Br  +  C  +  Da  log  ^  )  cob  Hq  cos  ^  z 

r  a'^ 

V  =  [  -  Ana  +  Bnr  log  a  +  C  ^7  6^^  lian  (log  I.  +  l)  + 

Fr  ]  sin  n0  cos  ^  z  (21) 

2  g 

w  =  [  AX.r  +  BaX  (^^—5-  -  _E  log  —  )  +  CXa  log  —  +  DXr 
Za'^  ^  ^ 

(log  —  -  1 )  +  La  ]  cos  n6  sin  —  z  (22) 

cl  & 

The  Facings 

Free  body  diagrams  of  a  facing  element  showing  the  forces  and 
moments  are  shown  in  figures  4  and  5.  It  is  necessary,  in  this  type 
of  problem,  to  include  components  of  forces  which  result  from  elastic 
deformation  of  the  element.  The  geometry  of  the  situation  is  such  that 
•»  it  is  difficult  to  write  equations  of  equilibrium.  It  is  safest  to  use 

results  obtained  from  a  mathematical  theory  of  thin  shells.  Such 
theory,  as  developed  by  Osgood  and  Joseph  (ref.  2),  when  applied  to 
cylindrical  shells  yields,  for  the  outer  facing  at  r  =  a,  the  following 
equations 


..aZjij  TR  oC-1_>3 


10 


2; 

N 

^  1 

X 

X 

X 

X 

X  1 

US 

^  /o 


0 

X 

a 

0 

X 

Ori 

N 

fsl 

u 

2 

X 

X 

X 

(M  N 

r< 

/t:  1  /r: 

a. 

X 

2  ^ 

/C 


o 

> 

rsj  N 

/X)  X) 


3  M 
(NJ  ® 

»X  /X 


Sw 
X 
X 


VaJX'  TR  oG-1_ 


ze  '  "*0  ^z^e  ^  30 


00 

fsj 


^1 


CD 


I  nj 


3  <M 
^1  fO 
»-H  I  rt 


I 


3 

rsj 

/C 


rsl 

/t) 


2 


>  N 
/C  /T) 


^  O 
/TD  ^ 


<NJ 


CD 

/O 


>  N 

/ol  /t) 


>1 

1 

N 

1  /D 

1 

rt 

CD 

1 

SJ 

(<q 

/Z) 

fXi 

<D 

N 

% 

1 

® 

II 

N 

< 

=> 

M 

ct 

s 

w 


..iiDI;  TR  6C-lj)3 


12 


As  is  customary  in  such  problems  the  stretching  of  the  mitlfl/e.  surface 

r  y  •*  -  , 

is  taken  into  account  by  substituting  in  equations  (23)  to  (28)  , 

(1  +  eg)  for  N^. 

Ne  (1  +  e  z)  ^0’ 

and  multiplying  the  surface  forces  by 

(1  +  *0  )  (1  +  € 

In  these  expressions 

1  u 

‘0  =1^+1 

*  z  • 

p 

N  and  No  of  equations  (23)  to  (28)  are  replaced  by  ( - + 

z  0  ^  ^  '  2Tr  (a  +  b) 

AN^)  and  [qa  (1  -  k)  +  ANg],  and  in  the  corresponding  equations  for  the 

,  P 

inner  facing  N^'  and  Ng'  are  replaced  by  ’  2tt  (a'  +  "b)  ^  ^z'  ' 

(qak  +  ANg').  This  is  necessary  because  the  forces  in  the  buckled  shell 
are  the  prebuckling  forces  plus  the  forces  due  to  buckling.  The  AN^, 
ANj.',  AN0,  and  ANg'  are  the  forces  due  to  buckling  which  are  later  to 
be  eupressed  in  terms  of  displacements. 

All  forces,  moments,  and  twists  other  than  the  prebuckling  forces 
are  considered  to  be  small  quantities  resulting  from  the  buckling.  The 
displacements  u,  v  and  w,  and  their  derivatives,  are  also  small  quantities 
resulting  from  the  buckling.  In  equations  (24)  to  (28)  products  of  any  two 

■v.ilJD  PH  bC-lj)3 


such  small  quantities  are  neglected.  Equation  (26)  is  solved  for  Qq 
and  equation  (27)  for  Q^.  The  results  are  substituted  into  equations 
(23),  (24)  and  (25).  This  gives: 


Ih 
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(These  equations  are  for  the  outer  facing.  A  aimilar  set  is  obtained  for 
the  inner  facing.)  Into  equations  (29),  (30),  and  (31)  expressions  for 
the  forces,  moments  and  twists  in  terms  of  the. displacements  (ref.  3) 
are  substituted.  The  surface  forces 

R  =  qa  -  (q^-h  +  Ao-j.)  r  =  a  outer  facing)  where  q  —  k  is  the 

prebuckling  stress  and  A  tr^  the  stress  due  to  buckling, 

0  =  -  (t  q)  (for  outer  facing),  and 

■u  r  =  a 

Z  =  -  (t^^)^  _  ^  (for  outer  facing), 

are  also  expressed  in  terms  of  u,  v  and  w  and  substituted  into  the  three 
equations. 

This  leads  to  three  equations  in  terms  of  u,  v,  and  w  for  the  outer  facing 
and  three  similar  equations  for  the  inner  facing.  The  equations  for  the 
outer  facing  are: 
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To  achieve  proper  interaction,  between  the  core  and  the  facings,  the 
displacements  of  the  middle  surfaces  of  the  facings  are  set  equal  to  the 
displacements  of  the  core  at  r  =  a  and  r  =  b. 

Thus  displacements  u,  v  and  w  in  equations  (32),  (33),  and  (34) 
are  replaced  by  equations  (20),  (21),  and  (22)  with  r  made  equal  to  a. 

In  this  manner  three  equations  in  six  arbitrary  constants  (A,  B,  C,  D, 
L,  and  F)  are  written  for  the  outer  facing.  In  a  like  fashion  three 
equations  are  written  for  the  inner  facing.  The  coefficients  of  the  six 
arbitrary  constants  are  shown  in  the  form  of  a  determinant  on  the 
following  page. 
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It  is  possible  to  find  simultaneous  values  of  <()  ^  and  <j)  ^  ior  which 

these  six  equations  will  be  satisfied  for  any  values  of  the  arbitrary 
constants.  Mathematically  this  means  that  for  such  a  combination  of 
loads  the  deflections  are  indeterminate.  The  shell  becomes  elastically 
unstable  and  the  loads  that  bring  about  this  condition  are  called 
critical  loads. 

Numerical  Computations 

A  literal  solution  of  the  sixth  order  determinant  for  the  eigenvalues 
is  not  feasible.  A  numerical  solution,  from  which  curves  may  be  drawn, 
is  possible  if  a  digital  computer  is  used.  A  CPC  Model  2  was  available 
to  make  computations.  Even  with  the  CPC  the  task  seemed  over¬ 
whelming.  If,  however,  E  is  made  infinite,  some  of  the  terms  of 

c 

determinant  vanish.  The  assumption  that  E^  is  infinite  is  common  in 

work  with  sandwich  construction  and  has  been  found  to  give  satisfactory 

results  in  most  cases.  The  sixth  order  determinant  with  E  made 

c 

infinite  is  represented  below: 
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^1 

®1 

^1 

0 

^1  ^1 

®2 

^2 

0 

^2  ^2 

0 

0 

^3  ^3 

^4 

0 

*•4  ^4 

A 

■5 

B 

5 

C 
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This  determinant 

is  then  reduced  to  a  fourth  order  determinant  shown 

below: 
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^  a 
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(F5k+F6)L3  - 

(C3I. 

c^)  A3 

^6 

(L5I+  L^)  C3 

(Ls^+L^)  F3 

The  determinant  is  then  programmed  for  the  CPC.  A  trial  and  error 
solution  is  made  by  substituting  values  of  4>  ^  or  ^  until  a  value  is 
found  that  will  make  the  determinant  zero.  This  was  done  by  finding 
values  on  each  side  of  zero  and  interpolating  to  find  the  eigenvalue. 

TH  bC— 


I 
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Discussion  of  Results 


Since  the  problem  is  solved  by  numerical  methods  the  results 
are  presented  by  the  curves  shown  in  figures  6,  7,  8  and  9.  The 

b  a 

values  of  _  =  0.  97  and  _  =  1,  000  were  used  for  all  of  the  curves, 
a  t 

Figure  6  is  a  family  of  curves  in  which  -(j>2  is  plotted  against 

i  E 

-  for  different  values  of  n.  In  these  curves  the  values  -  =  10,  000 

ma  Gj.0 

and -  =  1,  000  were  used.  Such  a  set  of  curves  is  used  in  the 

Grz 

following  manner. 

Knowing  ^  of  the  cylinder  one  picks  a  value  for  m  and  n.  A 

value  of  -4>,  is  determined  by  reading  above  ^ -  on  the  corresponding 

4  ma 

n  curve.  This  procedure  is  repeated  until  the  lowest  possible  value  of 
-<t)2  is  found.  The  axial  load  under  which  the  cylinder  will  buckle  can 
then  be  determined. 

The  curves  of  figure  7  differ  from  those  of  figure  6  as  a  result  of 
making  a  and  a'  zero.  This  is  equivalent  to  neglecting  the  bending 
stiffnesses  of  the  facings.  A  comparison  of  the  curves  of  figure  7  with 

i 

those  of  figure  6  shows  that  for  values  of  -  greater  than  0.  15  there  is 

ma 

little  difference.  It  can  be  concluded  that  only  for  very  short  cylinders 

need  the  bending  stiffnesses  of  the  facings  be  considered.  For  - - 

ma 

less  than  0.  15  the  curves  of  figure  7  approach 

Tn  oO— 2  < 


-‘t’2  = 


a  (1  -  tx^)  (1  -  1)^ 
Et  log  (1  +  |.) 


This  value  is  obtained  by  making  n  and  i  zero  and  expanding  the 
determinant.  Solving  for  N^,  and  replacing  log  ^  by  the  first  term 
of  its  series  expansion  shows  that 


N  =  -  G  . 

z  u  rz 

1  +  - 
a 


The  curves  of  figure  8  are  the  result  of  increasing  G^^  and  C^q 


tenfold.  The  value  of  -4>2  corresponding  to 


N 


z 


(a  -  b) 
1  +  ^ 


G 


rz 


appears  as  a  flattening  of  the  curve  in  the  region  of  i _ =  0.  01.  For 

ma 

f 

smaller  values  of  ~~  the  curve  rises  due  to  the  stiffness  of  the 
ma 

facings.  For  values  of  - -  greater  than  0.  1  the  curves  show  a 

ma 

considerably  lower  buckling  load. 

From  a  comparison  of  figures  6  and  8  it  appears  that  as  G^.^  is 
decreased  the  buckling  load  for  all  cylinders,  except  those  long  enough 
to  fail  as  an  Euler  column,  will  approach 


Nz  =  - 


=  -  -  b) 

b 


1  + 


rz 
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This  limit  has  been  recognized  (ref.  4)  as  the  critical  load  for  shells 
with  a  low  value  o£  It  should  be  noted  that  this  load  depends  only 

upon  the  thickness  and  the  modulus  of  rigidity  of  the  core. 

Figure  9  shows  curves  of  *4)  plotted  against  i _  for  different 

ma 

values  of  these  curves  was  taken  to  be 

X  -  ““I 

2(14|) 

This  represents  the  case  for  an  end  load  equal  to  qna^.  The  situation 

is  like  that  of  a  cylinder,  with  ends,  under  uniform  pressure.  The  ends 

of  course  stiffen  the  cylinder,  but,  if  the  cylinder  is  not  too  short, 

reasonable  results  can  be  expected.  Since  <j)  ,  decreases  as  i _  is 

^  ma 

increased  it  must  be  concluded  that  the  cylinder  will  buckle  with  m  =  1. 
The  critical  pressure  can  be  determined  by  reading  4>j  from  the  lowest 
n  curve. 

Conclusions 

Although  only  a  few  curves  were  drawn  it  is  apparent  that  this 
analysis  is  helpful  in  understanding  the  effect  produced  by  a  variation 
of  the  parameters  that  enter  the  problem.  Further  study  is  reqvdred 
before  it  will  be  known  whether  the  actual  buckling  load  may  be  predicted. 

It  is  felt  that  the  method  by  which  this  problem  is  solved  can  be 
applied  with  advantage  to  many  problems  of  sandwich  construction. 
Unfortunately  in  most  cases  a  numerical  solution  will  be  required. 
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Summary 


This  report  presents  a  theoretical  analysis  for  the  behavior  of  long,  circu¬ 
lar,  cylindrical  shells  of  sandwich  construction  under  axial  compressive 
loads.  The  analysis  is  designed  to  evaluate  the  effects  of  the  relatively  low 
shearing  moduli  of  sandwich  cores  on  buckling  stresses.  Families  of  curves 
are  presented  for  use  in  designing  shells  of  sandwich  construction  having  iso 
tropic  facings  and  orthotropic  or  isotropic  cores. 

The  results  of  the  theoretical  analysis  were  compared  with  those  obtained 
from  tests  on  a  series  of  curved  panels.  It  was  found  that  the  theory  applied 
reasonably  well  to  curved  plates  of  sizes  sufficient  to  include  at  least  one 
ideal  buckle.  Application  of  the  theory  thus  is  not  limited  to  long,  complete 
cylinders. 


—This  progress  report  is  one  of  a  series  prepared  and  distributed  by  the 
Forest  Products  Laboratory  under  U.  S.  Navy,  Bureau  of  Aeronautics 
Order  No.  NAer  01237  and  01202,  and  U.  S.  Air  Force  No.  USAF  18 
(600)  -70.  Results  here  reported  are  preliminary  and  may  be  revised 
as  additional  data  become  available.  Original  report  published  June  1952. 

2 

—Maintained  at  Madison,  Wis.  ,  in  cooperation  with  the  University  of  Wis¬ 
consin. 
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Introduction 


In  the  design  of  aircraft  and  guided  missiles,  it  was  found  necessary  to  de¬ 
vise  a  method  of  determining  the  stress  at  which  curved  sandwich  panels 
subjected  to  axial  compression  become  elastically  unstable.  It  is  known 
that,  for  thin,  homogeneous  materials,  a  curved  form  greatly  increases  the 
critical  load  as  compared  to  a  flat  sheet  of  the  same  approximate  size.  A 
similar  increase  may  be  expected  for  curved  sandwich  panels.  Although 
this  '.report  applies  primarily  to  sandwich  construction  for  aircraft,  the  re¬ 
sults  are  general  and  apply  to  any  structures  of  the  type  considered. 

This  report  presents  a  theoretical  analysis  of  the  behavior  of  long,  circular, 
cylindrical  shells  of  sandwich  construction  under  axial  compressive  loads 
and  an  experimental  confirmation  of  this  analysis  by  tests  on  curved  panels 
of  sufficient  size  to  include  at  least  one  ideal  buckle.  Thus,  these  panels 
are  assumed  to  simulate  the  action  of  complete  cylinders. 

The  buckling  of  a  homogeneous,  isotropic,  thin-walled  cylinder  was  treated 
by  von  Karman  and  Tsien  (16)—  and  by  Tsien  (13,  14)  in  related  papers. 

These  authors  assumed,  in  addition  to  the  wave  form  of  the  classical  theory, 
inward  buckles  of  diamond  shape  to  represent  the  characteristic  buckles  that 
are  actually  observed.  They  used  an  energy  method  to  determine  the  criti¬ 
cal  compressive  stress.  This  method,  in  which  only  diamond-shaped  buckles 
are  used,  was  applied  by  March  (7)  to  cylinders  made  of  plywood,  an  ortho¬ 
tropic  material.  Particular  attention  was  paid  to  the  effect  of  initial  irregu¬ 
larities  that  contribute  to  the  observed  scatter  of  experimentally  determined 
critical  stresses  of  both  isotropic  and  orthotropic  cylinders. 

In  this  report,  the  effect  of  shear  deformation  in  the  core  of  a  sandwich 
cylinder  is  taken  into  account  by  employing  an  approximate  "tilting"  method. 
This  method  was  used  by  Williams,  Leggett,  and  Hopkins  in  their  analysis 
of  flat  sandwich  panels  (18)  and  by  Leggett  and  Hopkins  in  their  analysis  of 
flat  sandwich  panels  and  cylinders  (4).  It  amounts  essentially  to  assuming 
that  the  transverse  components  of  shear  stress  are  constant  across  the  thick¬ 
ness  of  the  core.  The  form  of  buckles  assumed  by  Leggett  and  Hopkins  (4) 
in  the  cylinder  is  different  from  that  assumed  in  this  report. 

The  core  and  facings  are  taken  to  be  orthotropic,  with  two  of  their  natural 
axes  parallel,  respectively,  to  the  axial  and  circumferential  directions  of 
the  cylinder.  The  facings,  which  may  be  equal  or  unequal  in  thickness,  are 


—Underlined  numbers  in  parentheses  refer  to  Literature  Cited  at  end  of 
report. 


assumed  to  be  thin,  but  their' flexural  rigidities  are  not  neglected,  as  these 
may  be  of  importance  in  certain  cases.  All  stress  components  in  the  core 
are  neglected  except  the  transverse  shear  components.  It  was  pointed  out 
by  Reissner  (12)  that  the  stress  component  in  the  core  normal  to  the  facings 
may  be  of  importance  in  the  analysis  of  sandv/ich  shells.  Preliminary  calcu¬ 
lations  indicate  that  the  effect  of  this  component  is  small  in  the  problem  under 
consideration. 

As  was  done  in  work  described  by  Forest  Products  Laboratory  Report  No. 
1322-A  (7),  initial  irregularities  are  assumed  to  be  present  and  to  grow  un¬ 
der  increasing  compressive  load  until  buckling  occurs.  For  a  discussion  of 
this  important  matter,  the  reader  is  referred  to  that  report  and,  in  particu¬ 
lar,  to  the  observations  of  the  growth  of  artificially  produced  initial  irregu¬ 
larities.  Also  as  described  in  report  No.  1322-A,  a  large  deflection  theory 
is  used  to  take  into  account  the  nonlinear  support  associated  with  the  curva¬ 
ture  of  the  shell,  as  discussed  by  von  Karman,  Dunn,  and  Tsien  (17).  The 
derivations  of  the  differential  equation  for  a  stress  function  and  of  the  ex¬ 
pression  for  the  energy  of  deformation  are  extensions  of  the  analysis  used 
by  von  Karman  and  Tsien  (16)  for  the  homogeneous,  isotropic  cylinder  to  the 
sandwich  cylinder  composed  of  orthotropic  materials.  Suitable  modification 
is  made  for  the  effect  of  shear  deformation  in  the  core  of  the  sandwich. 


Theoretical  Analysis 


Choice  of  Axes  Notation 


The  choice  of  axes  is  shown  in  figure  1,  the  coordinate  y  being  measured 
along  the  circumference.  The  notations  for  stress  and  strain  are  those  of 
Love's  treatise  (5).  The  components  of  displacement  in  the  axial,  circum¬ 
ferential,  and  radial  directions,  respectively,  are  u,  v,  and  w,  the  latter 
being  positive  inward.  Since  initial  irregularities  of  the  cylindrical  surface 
are  assumed,  the  symbol  Wg  is  used  to  denote  the  initial  distance,  measured 
radially,  of  a  point  of  the  middle  surface  from  a  true  cylindrical  surface  of 
radius  r,  and  the  symbol  w  to  denote  the  corresponding  distance  at  any  stage 
of  the  deformation.  The  thickness  of  the  core  is  denoted  by  c  and  that  of  each 
facing  by  f ^  and  f^,  respectively. 


Extensional  Strains  and  Stresses 


Expressions  can  now  be  written  for  the  extensional  strains  uniform  across 
the  thickness  of  the  cylindrical  shell  and  for  the  corresponding  mean  mem¬ 
brane  stresses.  On  these  will  be  superposed  a  system,  of  flexural  strains, 
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and  the  energy  of  deformation  associated  with  each  system  of  strains  will 
be  found. 

The  extensional  strains  are  expressed  by  the  equations: 
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®xx  ■  ax  '''  2  V  ax  ) 
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®yy  “ay  z  Kdy 
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In  each  facing,  the  corresponding  stress  components  are; 
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where  E„  and  E  are  Young’s  moduli,  p  is  the  modulus  of  rigidity  for 

shearing  strains  referred  to  the  x  and  y  directions,  cr  and  cr  are  Poisson's 
=>  _  i  ’  xy  yx 

ratios,  and  X.  =  1  -  cr  <r  .  All  of  these  quantities  are  elastic  properties  of 
xy  yx 

the  facings.  Because  the  stress  components  Yy,  and  Xy  are  neglected 

in  the  core,  the  mean  membrane  stress  components  for  the  cylinder  are: 
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By  introducing  (10)  in  (8)  and  substituting  the  results  in  (11),  the  following 
differential  equation  for  F  is  obtained: 
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where; 
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It  is  readily  established  that  the  following  expression  represents  the  energy 
of  extensional  deformation  of  a  rectangular  portion  of  the  shell  with  edges  of 
length  a  and  b; 


W, 


-1/  J 


_  2  —  2  2(r  _  _  1—2 

BX  +  AY  -  — ^  Y  +  -  X^^ 

X  y  Ea  y  y 


dy  dx 


(14) 


Form  of  Buckles  and  Initial  Irregularities 

The  stress  components  X^^,  Yy,  and  Xy  in  (14)  are  derived  from  a  stress 
function  F,  satisfying  the  differential  equation  (12),  which  involves  deriva¬ 
tives  of  Wq  and  w  representing  the  initial  and  deformed  middle  surface  of  the 
shell.  For  w,  the  inward  radial  deflection,  the  following  form  will  be  chosen: 


—  =  g  +  6  cos^  (0y  -  ax)  cos^  Oy  +  ox) 
r 


where 

IT 


a  -  "  _ 

p  K’  “  r 
b  a 


(15) 


(16) 


The  nodal  lines  of  the  trigonometric  portion  of  equation  (15)  are  shown  in 
figure  2.  The  displacement  w  is  positive  inward.  In  equation  (16),  a  and  b 
represent  the  length  and  width,  respectively,  of  a  diamond.  The  initial  ir¬ 
regularities  will  be  assumed  to  have  the  form  (15).  This  is  done  for  the  pur¬ 
pose  of  simplifying  the  calculations.  Then  Wq  is  chosen  in  the  following  form: 
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(17) 


o  2  2 

~ '  So  ■*■ 

An  initial  flat  spot  on  the  surface  of  the  cylinder  could  be  described  roughly 
by  equation  (17).  An  initial  irregularity  was  introduced  here,  as  it  was  in 
report  No.  1322-A  (7_),  to  obtain  a  qualitative  description  of  its  influence  in 
causing  an  isolated  buckle  to  develop  in  its  vicinity.  If  the  initial  depth  of 
an  irregularity  of  the  form  (17)  is  very  small,  the  dimensions  of  the  area 
that  it  occupies  are  not  very  important.  For  this  reason  in  order'to  simplify 
the  calculations,  the  dimensions  a  and  b  in  equation  (17)  are  taken  to  be  the 
same  as  those  in  equation  (15).  From  the  qualitative  description  that  is  ob¬ 
tained  of  the  development  of  an  isolated  buckle,  conclusions  were  drawn  in 
report  No.  1322-A  that  led  to  the  derivation  of  the  final  formulas  from  the 
analysis  for  the  case  w^  =  0. 

The  details  of  substituting  (15)  and  (17)  in  (12),  of  obtaining  the  stress  func¬ 
tion  F  and  the  stress  components  Y  ,  and  X^.,  of  substituting  these  stress 

components  in  equation  (14),  and  of  related  operations  are  identical  with  the 
corresponding  operations  performed  in  report  No.  1322-A  (7).  Reference 
is  therefore  made  to  equations  (21)  and  (31)  of  that  report.  The  following 
differences  in  notation  should  be  noted: 

Notation  of  Report  No.  1322-A  Notation  of  R-esent  Report 
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X,  Y,  X, 


From  equation  (14),  the  energy  of  extensional  deformation  Wj  is  then  found 
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This  equation  is  equation  (31)  of  report  No.  1322-A  (7)  with  the  proper  values 
inserted  for  the  former  abbreviations  M  and  S.  The  quantities  p  and  cj 

represent  the  mean  compressive  stress  and  the  mean  circumferential  stress, 
respectively. 


If  n  is  the  number  of  buckles  in  a  circumference,  the  width  b  of  an  individual 
buckle  and  n  are  related  by  the  equation: 


b  = 


2Trr 


(19) 


Then: 


a  -  1  ^ 
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(20) 


It  will  be  convenient  to  denote  the  ratio  “  of  the  dimensions  of  a  buckle  bv: 

a  ^ 


a 


(21) 


In  the  expression  obtained  from  (18)  by  using  equations  (20)  and  (21),  let: 
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Equation  (18)  then  becomes: 
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FI exu ~  ^1  Energy  of  the  Shell 

To  determine  flexural  energy  of  the  shell,  the  following  simplified  expres¬ 
sions  for  the  changes  in  curvature  and  unit  twist  are  used: 


8  (w  -  w^) 


9^  (w  -  w^)  8^  (w  -  Wq) 


A  discussion  of  the  approximations  involved  will  be  found  in  a  paper  by 
Donnell  (1).  These  expressions  were  used  by  von  Karman  and  Tsien  (16) 
and  by  March  (7).  The  expressions  (30)  are  exactly  those  used  in  calculat¬ 
ing  the  flexural  energy  of  a  flat  sandwich  plate.  The  approximate  flexural 
energy  of  such  a  plate  was  found  by  March  (9)  and  by  Ericksen  and  March  (2) 
by  using  the  "tilting"  method  of  Williams,  Leggett,  and  Hopkins  (4,  18). 

In  this  method  it  is  assumed  that  any  line  in  the  core  that  is  initially  straight 
and  normal  to  the  vindeformed  plate  will  remain  straight  after  the  deforma¬ 
tion,  but  will  deviate  in  the  x  and  y  directions  from  the  normal  to  the  de¬ 
formed  plate  by  amounts  that  are  expressed  by  the  parameters  k  and  k'. 
These  parameters  are  determined  by  an  energy  method.  These  "tilting" 
factors  k  and  k'  are  introduced  as  well  as  two  quantities  q  and  q'  that  deter¬ 
mine  the  positions  of  the  surfaces  in  which,  respectively,  the  components  u 
and  V  of  the  displacement  in  the  core  vanish.  The  letters  k'  and  q'  replace 
h  and  r,  respectively,  of  report  No.  1583-B,  because  h  and  r  have  already 
been  used  in  the  present  report.  The  following  derivation  of  the  expression 
for  the  flexural  energy  follows  closely  that  used  for  the  flat  sandwich  panel 
in  Forest  Products  Laboratory  Report  No.  1583-B  (2),  to  which  reference 
is  made  for  further  details.  For  the  sake  of  simplicity  in  writing,  the  initial 
irregularity  w^  will  be  for  the  present  taken  equal  to  zero.  It  will  then  be 

introduced  in  the  final  steps  by  replacing  w  by  w  -  w^. 

The  components  of  displacement  in  the  core  (fig.  3)  are  taken  to  be: 

=  -  k  (^  -  q) 

v^  =-  k’  (^  -  q')|^  .  (31) 

W^  =  w  (x,  y) 
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Thus  ^  =  q  denotes  the  surface  in  which  the  components  of  displacement  in 
the  X  direction  vanish  and  k  is  the  parameter  describing  the  inclination  in 
the  X  direction  of  the  respective  plane  sections  to  the  normal  to  the  deformed 
surface.  Similarly  q'  and  k'  are  related  to  the  displacements  in  the  y  direc¬ 
tion.  These  four  quantities  are  to  be  determined  in  such  a  way  that  the  flexur¬ 
al  energy  associated  with  a  prescribed  deflection  w  is  a  minimum. 

To  arrive  at  expressions  for  the  components  of  displacement  in  the  facings, 
it  is  noted  that  the  continuity  of  the  displacement  at  the  facing-to-core  bonds 
requires  that  the  components  (31),  evaluated  at  ^  =  0  and  ^  =  c,  shall  be 
those  at  the  inner  surfaces  of  the  facings  fj  and  f2,  respectively.  Within 
each  facing,  the  components  of  displacemeTit  are^ssumed  to  be  such  that 
a  straight  line  initially  normal  to  the  undeformed  surface  of  the  plate  will 
be  straight  and  normal  to  the  deformed  surface.  Accordingly,  the  compon¬ 
ents  of  displacement  in  the  facings  fj  and  f^,  respectively,  are: 

/  , 

v,  =  (k'q'  -  t)  (32) 


Wj  =  w  (x,  y) 
and 

uz  =  -[k  (c  -  q)  +  C  -  c]|^ 
v^  =  -  [k'  (c  -  q')  +  ^  -  c]  I 


W2  =  w  (x,  y) 

The  components  of  strain  in  the  core  c  and  facings  fj  and  f2  will  be  denoted 
by  the  superscripts  c,  1,  and  2,  respectively.  ~ 

From  (31),  the  transverse  shear  strains  in  the  core  are: 

(3.) 


The  effect  of  the  remaining  strains  in  the  core  is  assumed  to  be  negligible. 


In  finding  the  strain  energy  of  the  facings  in  the  bending  of  the  plate  (or  shell), 
it  is  convenient  to  consider  the  components  of  strain  in  the  facings  to  result 
from  the  superposition  of  two  states  of  strain.  The  first  of  these  consists  of 
the  membrane  strains  in  the  facings  associated  with  flexure,  that  is  the 
strain  in  their  middle  surfaces.  From  (32)  and  (33),  these  strains  are  found 
to  be: 
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(1) 


=  (kq  +  7-) 


3^w 


-XX  -  --  2  'g^2 


e--/  =  (k'q’  +iL)^ 

yy  2%y2 


M) 


=  (kq  +  k  q  +  £^) 

and 

'k(c-q)  +  7 


9^w 

8x9y 


9^w 


9x‘ 


(35) 


yy 

J2) 

-xy 


-  l^k  (c  -  q)  +  k  (c  -  q')  + 


9  w 
9x9y 


(36) 


The  second  state  of  strain  in  the  facings  is  that  associated  with  their  bending 
about  their  own  middle  surfaces.  This  state,  in  either  facing,  has  the  com¬ 
ponents: 


I  „ '  9  w  I 

e  =  -  ?,  - e 


XX 


9x" 


yy 


,<  9^w  '  -,,1  9^w 

^  e  =  -  2^  - 

,2  xy 


(3  7) 


9y‘ 


9x9y 


where  ^  is  measured  from  the  middle  surface  of  the  facing  under  considera¬ 
tion. 

The  strain  energy  in  the  core  or  facings  is  given  by  the  expression  (6,  8). 


U 


MI 


e^xx  -t-  E  e^yy  -I-  2E  cr  e  e 

X  y  X  yx  XX  yy 


+  ^i^xy  e^xy  +  +  ^^^Cx  ®\x  '^y 


(38) 


where  for  the  material  under  consideration  (core  or  facing),  k  =  1  -  cr  cr  ; 

'  xy  yx 

Ejj  and  E^  are  Young's  moduli;  h^y’  ^yC’  *^t.x  moduli  of  rigidity; 

and  cr^y  and  Cyx  are  Poisson's  ratios.  Primed  letters  will  denote  the  elastic 
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constants  of  the  core  material  and  unprimed  letters  will  denote  those  of  the 
facing  material.  The  integration  indicated  in  formula  (38)  is  to  be  carried 
out  over  the  area  OABC  of  figure  2  and  the  thickness  of  the  core  or  facings. 


The  energy  in  the  core  is  obtained  by  substituting  expressions  (34)  into  (38), 
the  remaining  strains  in  the  latter  formula  being  neglected  as  previously 
stated.  After  integrating  with  respect  to  ^  over  the  thickness  of  the  core, 
the  expression  for  the  energy,  denoted  by  U^.,  is 


U 


a  b  _ 

If' 


M-  (1  - 


+  y^  (1  -  k  ) 


/9w\ 

\9y/ 


dy  dx 


(39) 


The  strain  energy  in  the  facings  associated  with  the  membrane  strains  is  the 
sum  of  the  energies  obtained  from  (35)  and  (36).  With  the  substitution  of 
these  expressions  into  (38)  one  obtains,  after  integration  with  respect  to 
the  following  expression  which  is  denoted  by  Uj^. 


The  strain  energy  in  the  facings  associated  with  the  flexural  strain,  is 

obtained  by  substituting  expressions  (37)  into  (38)  and  integrating  over  the 
volume  of  each  facing.  After  integrating  with  respect  to 


Tr(  c 


Now  in  all  of  the  expressions  (39),  (40),  and  (41)  replace  w  by  w  -  w^.  The 

flexural  energy  W2  of  the  region  OABC  of  the  shell  is  the  sum  of  Uj^,  Up, 

and  U  • 
c 

For  equilibrium,  the  "tilting"  factors  k  and  k'  and„_the  ordinates  q  and  q'  of 
the  neutral  surfaces  are  to  be  chosen  so  that  the  total  energy  is  a  minimum. 
But  these  factors  appear  only  in  the  flexural  energy  ^2.’  Hence,  they  must 
be  chosen  to  satisfy  the  conditions 

aw,  aw,  aw,  aw, 

-—±-  =  0,  ^  =  0,  - 4  =  0,  - ^  =  0 

a(kq)  a(k'q')  Bk  Bk' 

By  proceeding  exactly  as  in  report  No.  1583-B  (2),  the  quadratic  form  (A14) 
of  that  report  with  k'  and  q'  replacing  h  and  r,  respectively,  is  obtained  for 
W^.  The  coefficients  B^^  in  equation  (A14)  are  defined  by  equations  (A15)  in 

terms  of  the  quantities  Aj,  which  are  defined  by; 


o  o 
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On  factoring  out  the  comnnon  factor,  it  is  found  that: 

=  r^  (6  -  6^)^  abj^  Bj’  (kq)^  +  2B2'  (kq)  {k'q')  +  B3’  (k’q')^ 

+  2B^'  (kq)  k  +  2B5'  (kq)  k’  +  ZB^'  (k'q')  k  +  2B^’  (k’q')  k' 

+  B^’  k^  +  2Bg'  kk'  +  Bg'  k’^  +  2Bjq'  (kq)  +  ZB^j’  (k’q’) 

+  2BJ2'  k  +  2BJ3'  k'  +  Bj^'  +  Bjg'j  (47) 

t  f 

where  the  quantities  B^  are  defined  in  terms  of  the  quantities  Aj  by  equa¬ 
tions  (A15)  of  report  No.  1583-B,  each  Aj  replacing  the  corresponding  Aj 

in  those  equations.  The  quantity  in  brackets  in  equation  (47)  corresponds  to 
2u'  in  report  No.  1583-B.  (Note  that  equation  (A22)  of  that  report  should 
read  P  =  2U'). 

It  is  easy  to  see  that  the  steps  of  imposing  the  conditions 


aw 


a(kq) 


^  =  0. 


aw. 


a(k'q' 


aw, 

0,  — 2  _ 

ak 


0, 


aw 


ak' 


2  _ 


=  0 


and  of  determining  kq,  k'q',  k,  and  k'  and  substituting  their  values  in  the  ex¬ 
pression  (47)  for  2W2  are  identical  with  those  taken  in  report  No.  1583-B  (^) 

and  that  2W2  is  equal  to  the  right-hand  member  of  equation  (A25)  of  that  re¬ 
port  multiplied  by  r^  (6  -  Sg)^  ab.  It  is  concluded  from  equations  (A26),  (A27), 
and  (A28)  of  report  No.  1583-B  that: 


W^  =  1/2  r^  (6  -  6^)^  ab 


Aj  +  2A2  +  A3'  +  (Aj  A3'  -  (-^  + 


A  4>  A  -4> 

1  +  — —  +  — —  + 


4,^  (Aj’A3’  -  A^’^) 


+  (A^'  +  ZA^'  +  A3’)j 

where  ,  and  b  are  defined  by: 


(48) 
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(j)  = 


^1  ^2 
'1  ^'2 


and  Aj'  by: 


^1  ={  (3  E^a'^  + 


A^  -  ^  (^x°’yx  ■*■  ^Kxy^  ^ 


A3'  =  ^  i'i 


3C.2  k’: 


Note  that 
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^  =T^ 


(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


wh  ere 

K.  =  3E  z"*  +  3E„  +  2(E  cr  +  2\ii  )  z^ 

*+  X  q  '  X  yx  ^xy' 

I  t 

and  introduce  the  following  abbreviations  in  the  expressions  for  Aj  ,  A2  , 


» 

and  A3  ; 


wh  ere: 


.'•'  rpt  C-*  i_5 


(56) 


(57) 


d  =  3E  +  \|JL  z^ 
1  X  ^xy 


*^2  ~  (^x  ^yx  '*' 

d3  =3Ey+  K^^yZ^ 
Also, 


o  2  t  o  q2  I 

3of 


(58) 


8 
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Substituting  these  expressions  for  A  ^  equation -(48)  becomes 


^2  =  ;^  r2  (6  -  ay  ab/3^ 


f  h "  ^  '■‘1-3  -  "  jT^)] 

|l  ^  8/32dj<i,  ^  8/3^d3<},  ^  640'^4>^(djd3  -  d^^) 


9\^)jl'„  p  -z^ 
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+  IfK, 


(59) 


The  following  transformations  are  made  by  using  equations  (20)  and  (22): 
8/3^4)  _  2  n^4)  _  2q4)  _  f)  c  ' 

I  ”  I  2  ~  '  ”  ®x 

3Xp  3\p  r^  3Xp  rh  E 


;x 
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(60) 


and 
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3Xp  ^x 
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where 

2E  4> 

c - ii _ 

y  ~  3Xp'  rh 

y^ 


(61) 
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The  coefficient  of  the  expression  in  brackets  in  equation  (59)  is  also  trans¬ 
formed  by  using  equations  (20)  and  (22).  The  expression  for  becomes: 

(5  -  S  ab  h  *  t  'S  " 

w  = - - -  < - 

32Kr  ^djS^  qd3S  (d^d3  -  d^^)  S^S 


+  (Ij/h  ) 


E  ^  z‘ 
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W,  =  e.  ab  (e  -  iy  ^ 


where 
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[(l/h^)[K4+(did3-d/)^(^+S  )] 
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qdjS^  qd3S  q2(djd3  '  ^2^)  ^  S 
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Virtual  Work  of  the  Compressive  Load 

Exactly  as  in  equation  (35)  of  report  No.  1322A  (7),  the  virtual  work,  V/3, 

of  the  compressive  load,  calculated  for  the  region  OABC,  figure  2,  is - 

found  in  the  notation  of  the  present  report  to  be: 

W3  =  abh  I^Bp^  +  pc^  +  ^  r^  (6^  -  6^^)  pj 

=  abh  ^Bp^  +  pcj  +  eg  q  (|^  -  4^^)  p  “j  (6E 


where 


J.  2 
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It  will  be  convenient  to  consider  the  mean  energy  per  unit  volume  of  the 
cylindrical  shell.  Hence: 


W  =  (Wj  +  -  W3)/abh 


(67) 


In  accordance  with  equations  (29),  (63),  and  (65): 
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9"  (e  -  e/)  -  e^r,  (e"  -  e,")  (e  -  e„)  +  e3  (e  -  e,)' 


+  e 
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(68) 


The  Buckling  Stress 

For  equilibrium,  the  derivatives  of  W  with  respect  to  the  various  parameters 
8.’  ^1’  i.’  H’  vanish.  From  the  condition 


aw 

-  =  0,  it  follows  that  c,  =  0 

9c,  1 


(69) 


Now  Cj  denotes  the  mean  circumferential  stress.  The  parameter  appears 

only  in  the  expression  for  c^  as  given  by  equation  (30)  of  report  No.  1322-A 

(7^).  The  fact  that  c^  vanishes  implies  that  g^,  which  describes  a  uniform 

radial  expansion  of  the  cylinder,  takes  on  such  a  value  that  the  mean  circum¬ 
ferential  stress  vanishes.  Further  consideration  of  the  parameter  g  is  not 
necessary. 

9W 

From  the  condition  — —  =  0,  it  follows  that: 

95 


P  = 


-  io)  h 
e,-  r 


(70) 


where  p,  as  previously  noted,  is  the  mean  compressive  stress. 
Let: 


'1 


li  =!i  =1± 

E^’  ^2  E^’  ^3  E^’  ^4  E^ 


Then  (7  0)  can  be  written: 
p  =  E 


24 


^4  4 


(4  -  4,)  h 


(71) 


(7  2) 


The  mean  compressive  strain  €  is  expressed  by: 


.'•^vTD  TR 


-  € 


As  in  report  No.  1322-A  (7), 

In 

64 


„  ,  3ti  2  /t2  t  2,  h 

'  =  Bp  +  z  (I  -  )  T 


it  is  found  that: 


The  further  equilibrium  conditions 


(73) 
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0  and 


9W 

9z 


=  0 


are  to  be  satisfied.  The  first  of  these  is  associated  with  the  number  of 
buckles  in  a  circumference  or  with  the  width  of  an  individual  buckle,  and  the 
second  with  the  ratio  b/a  of  the  width  of  a  buckle  to  its  length.  It  is  not 
analytically  feasible  to  use  these  conditions  in  connection  with  equation  (68). 


The  following  method  of  arriving  at  the  critical  value  of  p  is  based  upon  an 
extended  discussion  in  report  No.  1322-A  (7).  Briefly,  it  was  considered 
that  an  isolated  initial  irregularity  would  increase  in  size  and  depth  with  in¬ 
creasing  mean  compressive  stress  p.  It  was  therefore  considered  that  the 
load-mean  compressive  strain  curve,  with  p  as  a  function  of  *,  for  a  given 
small  initial  depth  of  irregularity  would  be  the  envelope  of  the  family  of 
curves  for  p  as  a  function  of  t  ,  drawn  for  a  series  of  values  of  ri  by  combin¬ 
ing  equations  (71)  and  (72).  On  taking  into  consideration  the  possibility  of 
jumps  from  one  energy  level  to  another,  it  was  concluded  that  the  critical 
values  of  p  would  scatter  considerably,  as  they  actually  do  in  test,  depend¬ 
ing  upon  the  depth  of  the  initial  irregularity  and  the  characteristics  of  the 
loading  process.  It  was  noted  that  the  value  of  p  at  the  relative  minimum 
point  on  the  envelope  of  the  curves  for  p  as  a  function  of  t,  drawn  for  =  0, 
was  intermediate  among  the  possible  critical  values  of  p.  This  minimum 
was  accordingly  chosen  as  the  "theoretical"  critical  stress,  because  it  could 
be  conveniently  determined  by  finding  a  relative  minimum  of  p  as  a  function 
of  ^  and  r|.  It  is  necessary  to  employ  numerical  methods  to  determine  the 
relative  minimum  value  of  p. 

In  report  No.  1322-A  (7),  the  aspect  ratio  z  of  the  buckles  was  assumed  on 
the  basis  of  experimental  observations  before  the  minimization  of  p  was 
undertaken.  Here,  because  of  the  influence  of  shear  deformation  in  the  core, 
a  suitable  value  to  assign  to  z  can  not  be  estimated. 

Equation  (7  1),  with  =  0,  can  be  written  in  the  form; 

p  =  KE^|  (7  4) 

'Ir?  c.:: 


where 


The  mean  compressive  stress,  pf,  in  the  facings  is  related  to  the  mean  com¬ 
pressive  stress  in  the  shell  by  tKe  equation 

p .  „6) 

On  recalling  the  definition  of  it  is  seen  that  equation  (7  4)  can  be  written 

Pf  =  KE^7  (77) 

For  a  relative  minimum  of  p^,  the  condition 
=  0  must  be  satisfied. 


From  this  condition,  it  follows  that: 


The  substitution  of  this  value  of  ^  in  (75)  yields 

In  equation  (7  9).  K  is  a  function  of  n  and  z,  which  occur  in  the  definitions  of 

■qj  7,1  v_,  y,,  and  y  .  By  using  the  definitions  of  the  quantities  E  ,  E,  , 

—  1  ’  4  3  4  _ ^  P 

and  that  appear  through  the  symbols  A,  B,  and  C  in  the  equations  (26), 
(27),  and  (28),  the  following  expressions  are  obtained  for  y^,  y^,  and  y^ 
(see  equations  (71),  (26),  (27),  and  (28): 
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Ey  ^x^y 

In  obtaining  the  expression  for  from  equations  (6  3)  and  (70),  it  is  con¬ 
venient  to  introduce  the  notation — 

T  =  3z  +3  (Ex‘’‘yx  ^ 

so  that 


K4  =  T 


E  f,  +  f? 
a  1  <2 


32\h  (fj  +  f^) 
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Buckling  Stress  of  Sandwich  Constructions  with  Isotropic 
Facings  and  Orthotropic  or  Isotropic  Core 

For  isotropic  facings,  considerable  simplifications  can  be  made.  In  this 


Ex  =  Ey  =  (E^<r^^  +  2Xp^y)  =  E,  =  (x  =  E/2  (1  +  a) 


O’  =0*  =  O’ 

xy  yx 


yx  ^xy 


Then 
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V 


(87) 


V,  = 


1  +  z 


17  z' 


—  — . -.I.  -j-  '■■■■■  .  X  . . -  ■  '  f  - . . .  ■■ 

1  4096  512  (z^  +  9)^  512  (9z^  +  1)^  2048  (1  +  z^)^ 

4  -  —  - 


1 


y  “ 

2  512  32  (1  +  z^)^ 

~  256  32  (1  +  z^)i. 


4  -a  -3  \  u  2 


32\  (f^  +  f^) 


1  (d,d,  -  T1 

I  +  - i-LZl—  (JL  +  s  J 

21  ?  ?  2  y 

(3^  +  3  +  2z-^)  - 


TidjSx^  -ndjSy  ^  n^SxSy  (djd3  -  d^^) 

1  +  ;>  ?  7 
Ez^  E 


+  h 


(88) 

(89) 

(90) 


where 

T  =  3z'^  +  3  +  2z^ 


2  4)  E 

~  3\  IJ.'  tdT 


2  4i  E 
3X  u'  „  rh 


(91) 

(92) 


After  some  manipulation  involving  substitution  of  expressions  for  Yj,  y^, 

y^,  and  Y4,  formula  (79)  for  K  for  sandwich  construction  with  isotropic 
facings  and  orthotropic  core  can  be  written  as; 


M, 


K  = 


21 


^  3X  h^  (fj  +  f^) 


M2T1  +  M3T1  Sj 


Ll  +  +  Mg-n^Sx^ 


-  ^29 


1 


where 

64  , 

"  3^2  ^^3  ■  32y 
T 

^2  =  ;2 


(dids  -  )  1 


(93) 

(94) 

(95) 

(96) 
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O'. 


rh 

\  \ 


V 


M4  =  (-4+  ^3®)  ^ 

z 


(d^d3  -  d/)  e 


Sy  H 

0  =  -^  or  0  = 


d,  =  3Ez  +  ”  (1  -  <t)  Ez^ 


d,  =  E(rz  +  —  (1  -  O')  Ez‘ 


d3  =  3E  +  -  (1  -  (t)  E2 


If  O'  is  taken  to  be  “  then: 

4 


dj  =  3Ez^  (z^  +  4 


d3  =  3E(~+1) 


If  also  Vj,  Y2'  Y3.  and  are  expressed  in  terms  of  z  and  0  (  Eq.  87,  88,  89, 

and  90)  then  the  following  expressions  can  be  used  in  formula  89: 

r  ,  3  -1  2 


2  3  -^  + - 

M  =_J_+_i5 _ Lfe4z  4(1 

^  12z^  3(1  +  z^)^  1  +  z^  4  z^*  4 


16(z^  +  9)^  16(9z^  +  1)“^  64(1  +  z^) 


2  .  1  \2 


2  3 

M,  =  3z‘^  +  2  + 


M3  =  (9z^  +  70z^  +  9)  (^  +  0) 


I 


(107) 

1 

(108) 

fl 

i 

h 

i 
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(109) 


M4  =  g  [sz^  +  1  +  {z2  +  8)  e] 


=![■ 


9z  +  70z  +9  0 


(110) 


For  constructions  for  which  the  shear  deformation  in  the  core  is  negligible, 
as  it  is  when  |jl'  is  very  large  and  0  is  finite,  S  may  be  taken  equal  to  zero. 

Then  expression  (93)  can  be  minimized  with  respect  to  tj,  resulting  in: 


Ko  =  2 


'2M1M2  (I  +  If) 

3X  (fi  +  f^)  h^ 


(111) 


Thus,  Kg  is  a  function  of  and  M2  and  the  stiffness  of  the  sandwich.  It 
was  found  by  computation  that  a  relative  minimum  of  MjM2  =  0.24  occurs  at 
z  =  0.95.  The  minimum  buckling  stress  is  then  proportional  to: 


*^0  =  5^ 


(112) 


where 


'X  (fi  +  f2) 
I+If 


(113) 


By  letting  N  =  7^,  the  following  expression  can  be  written  from  equation  (89) 


for  constructions  having  any  value  of  S„: 


5M.Q, 

N  =  — : - + 


M2II  + 


4ii  6Q,  (1  +  Q2)  ,  ,  w  c  4.  XX  2c  2  ^  ^2^2’! 

'  1  '  2'  1  +  M4TiS^  +  M5TI 


where 


(114) 


Q2=T 


(115) 


It  was  found  in  Forest  Products  Laboratory  Report  No.  1505  (10)  that  the 
values  of  I  +  l£  and  can  be  expressed  as  follows: 
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(116) 


I  +  I. 


1  r  3  3  12cd^  ' 

=  n  ■  TT^_ 


(h  -  c)^  (h  -  c)d^ 
“  48  4 


(117) 


where 


d  = 


fi  -  fz 


(118) 


After  substituting  these  expressions  in  the  formulas  for  and  and 
simplifying  Qj  and  become: 


1Z\  (1  -  c/h)^ 


(1  -  c^/h^)  (1  -  c/h) 


IZcd 


(119) 


^2  = 


d  ^ 

i  (1  -  c/h)^  +  3  (.1  -  c/hf  ^ 


(1  -  c^/h^)  (1  -  c/h)  -  ^  (1  -  c/h)^  -  3(1  -  c/hf  ^ 

^  h^  h 


(120) 


In  equation  (114)  Mj^,  M3.  M^,  and  depend  upon  z_  and  0^,  and  Qj^ 

and  Q2,  depend  upon  c/h  and  d/h.  Formula  (1 14)  can  then  be  written  with 

appropriate  values  of  0,  c/h,  and  d/h  and  then  a  relative  minimum  value  N 
found  by  choosing  a  series  of  values  of  z  and  r|.  The  facing  stress  at  which 
buckling  will  occur  is  then  given  by: 


Buckling  Stress  of  Sandwich  Constructions 
with  Isotropic  Facings  of  Equal  Thickness 
and  Orthotropic  or  Isotropic  Core 

Factors  in  formula  (114)  can  be  simplified  for  sandwich  constructions  hav¬ 
ing  facings  of  equal  thickness.  Then  d  =  0  and  after  simplification 


c  /h 


3X _ 

■  +  c  /h  +  1 


(122) 
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and  finally  formula  (114)  becomes 


5Ql  J^l  ,  (1  +  c/hf 

^  “  4  I  "H  24\ 


M2il  + 

1  +  M,tiS  +  ^ 

4  X  5  X 


(1  -  c/h)^ 
72X 


(123) 


(124) 


The  buckling  load,  which  is  proportioned  to  N  is  obtained  by  finding  the  lowest 
relative  minimum  of  expression  (124)  with  respect  to  t|  and  z.  Expression 
(124)  can  be  minimized  by  taking  a  derivative  with  respect  to  t|  and  setting 
the  derivative  equal  to  zero.  This  leads  to  a  sixth  power  equation  in  q. 
Minimum  roots  of  q  with  respect  to  z  and  for  various  values  of  i  c/h,  and 
0  were  determined  by  means  of  a  digital  computer.  Minimum  values  of  N 
at  various  values  of  and  for  c/h  equal  to  0.9,  0.8,  0.7,  and  6  equal  to 
0.4,  1.0,  and  2.5  are  given  in  Table  1  and  shown  as  functions  of  Sjj  in  figures 

4,  5,  and  6.  Also  included  in  the  table  and  figures  are  values  of  N  for 
c/h  =  1.  These  values  represent  sandwich  constructions  for  which  the  stiff¬ 
ness  of  the  individual  facings  are  assumed  to  be  zero.  Although  no  actual 
constructions  can  be  made  of  this  type,  the  values  can  be  considered  as 
representing  the  limit  for  constructions  having  extremely  thin  facings.  These 
values  of  N  were  obtained  as  follows.  Substitution  of  c/h  =  1  in  equation  (124) 
for  N  leads  to 


N>-vr 


l2q+  M3q 


1  +  M^qS^  +  Mgq^S^^^ 


(125) 


which  has  one  relative  minimum  value  for  q  =  oq, 
given  by 

N  =  1Z^S^  M5 


This  minimum  value  is 


(126) 


Substituting  in  this  equation  the  values  of  M3  and  Mg  given  by  equations  (108) 
and  (1 1  0)  yields 


5  ^  +  ®) 

12/\  Sj,  0 


(127) 
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which  is  minimum  for  z  =  oa  This  minimum  value  (for  cr  =  1/4;  hence 
K  =  15/16)  is  given  by 


N  = 


0.431 


(128) 


Substitution  of  this  value  of  N  in  equation  (118)  and  using  the  value  of  for 
ff  =  f2  and  c/h  =  1  leads  to  the  following  limiting  expression  for  p^; 


4Eh 

"  10^  r 


5  _  h  • 

12jX‘  Ef  "  2f 


(129) 


For  values  of  ranging  from  0  to  about  0.6  it  was  found  that  equation  (128) 
did  not  give  loWest  minimum  values.  In  this  range  of  the  minimum  values 

were  obtained  from  equation  (124)  by  use  of  a  digital  computer. 

The  value  of  N  given  by  equation  (128)  and  the  value  of  the  stress  given  by 
equation  (129)  are  independent  of  the  radius  of  the  cylinder  and  are  the  usual 
critical  values  associated  with  shear  instability  of  the  core  (15). 

The  value  of  6  of  0.4  and  its  reciprocal  2.5  were  used  in  the  calculations  be¬ 
cause  they  apply  to  honeycomb  cores  oriented  with  the  weak  direction  and 
the  strong  direction  parallel  to  the  length  of  the  cylinder.  It  has  been  noted 
from  figures  4  and  6  that  in  the  range  of  small  values  of  where  N  is  inde¬ 
pendent  of  c/h,  the  orientation  of  the  core  makes  little  difference  in  the 
value  of  N. 


Application  of  Theoretical  Results 


The  compressive  facing  stress  at  which  buckling  of  cylinders  of  orthrotropic 
sandwich  construction  occurs  is  given  by  equation  (77). 


Pf 


=  KE  - 
X  r 


where  Ej^  is  modulus  of  elasticity  of  facings  in  axial  direction,  h  is  sand¬ 
wich  thickness,  r  is  mean  radius  of  curvature,  and  K  is  given  by  formula 
(79)  as  “  — 

>3  9^2^ 


K  = 


64 

3z2 


q  32Y^q 


Y4TI 
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where  values  of  "y^i  S-nd  functions  of  z  according  to  equations  (80), (81), 

and  (82)  and  y^  is  a  function  of  t|  and  z  according  to  equation  (86)  and  K  is 
taken  as  the  least  relative  minimum  with  respect  to  t|  and  z. 


For  sandwich  constructions  having  isotropic  facings  of  unequal  thickness 
and  orthotropic  core  K  is  given  by 


K  = 


4N 

5Qi 


where  N  is  given  by  equation  (114)  as 


5M,Qj 

N  =  + 


4ti  6Qj  (1  +  Q2) 


M2T1  + 


1  + 


+  Q2M2q 


where  Qj^  and  Q2  are  given  by  equations  (119)  and  (120)  and  Mj,  1^2, 

and  Mg  are  functions  of  z  according  to  equations  (106),  (107),  (108), 

(109),  and  (110)  and  N  is  taken  as  the  least  relative  minimum  with  respect 
to  q  and  z. 


For  sandwich  constructions  having  isotropic  facings  (Poisson's  ratio  1/4)  of 
equal  thickness  and  orthotropic  core  such  that  0  =  0.4.^  or  0  =  2.5— or  isotropic 
core  (0  =  1.0)  equation  (1  14)  for  N  has  been  solved  for  c/h  =  1.0,  0.9,  0.8, 
and  0.7.  Values  of  N  for  various  Sj^  values  are  given  in  Table  1  and  in 

graphs  in  figures  4,  5,  and  6.  Then  the  critical  facing  stress  is  given  by 

4N  h 
Pf  “  5Q^  ^  r 


where 


^  ^yc^/h^  +  c/h  +  1 

and  N  is  given  in  terms  of  S  where 

_  =  X 


4 

“These  ratios  for  0  were  chosen  as  representative  of  honeycomb  cores  such 
as  were  evaluated  in  Forest  Products  Laboratory  Report  No.  1849. 
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s 


X 


l6cfE 

45rhu.', 

bx 


and  c  is  core  thickness,  f  is  facing  thickness,  E  is  modulus  of  elasticity  of 


facings, 
and  LL  , 


h  is  sandwich  thickness,  (c  +  2f),  r  is  mean  radius  of  curvature, 
is  modulus  of  rigidity  of  core  associated  with  shear  strains  in  the 


axial-radial  plane. 


The  graphs  can  be  used  with  little  error  for  determining  N  for  constructions 
having  facings  of  unequal  thickness,  provided  is  calculated  using  formula 

(60)  and  Qj  is  calculated  using  equation  (119). 


The  analysis  may  be  extended  to  apply  at  stresses  greater  than  the  propor¬ 
tional  limit  stress  of  the  facings  by  use  of  an  appropriate  tangent  or  reduced 
modulus  of  elasticity  for  the  facings.  This  entails  a  "trial-and-error ''  solu¬ 
tion  involving  use  of  the  tangent  or  reduced  modulus  in  the  quantity  S  and 
elsewhere  until  the  resultant  facing  stress  is  compatible  with  the  stress- 
modulus  curve. 


Results  of  the  theoretical  analysis  fall  approximately  into  three  zones,  de¬ 
pending  upon  whether  there  is  no  shear  deformation  in  the  core  (Sj^  =  0), 
some  shear  deformation  in  the  core  (small  values  of  S^),  or  considerable 
shear  deformation  in  the  core  (large  values  of  S^). 

For  no  shear  deformation,  the  buckling  stress  is  determined  essentially  by 
means  of  the  isotropic  or  orthotropic  theory  (depending  upon  facing  proper¬ 
ties)  with  the  stiffness  determined  by  considering  the  spaced  facings  of  the 
sandwich. 


For  large  shear  deformations,  the  critical  stress  is  associated  with  insta¬ 
bility  of  the  core  in  shear.  This  has  been  observed  for  sandwich  construc¬ 
tions  in  general  (15),  and  it  has  been  found  that  the  mean  critical  stress 
thus  determined  is  the  same,  regardless  of  the  original  assumption  of  the 
buckled  shape.  The  smallest  value  of  at  which  the  critical  stress  is 

determined  by  shear  instability  of  the  cere,  however,  is  greatly  affected  by 
the  assumed  form  of  the  buckled  shape.  The  inclusion  of  the  stiffnesses  of 
the  facings  If  gives  rise  to  the  family  of  curves  for  different  values  of  c/h, 

as  shown  in  figures  4,  5,  and  6,  instead  of  a  single  curve.  K  the  stiffnesses 
of  the  individual  facings  had  been  neglected,  one  curve  only,  that  for  c/h  =  1, 
would  have  resulted.  The  percentage  increase  in  buckling  stress  due  to  the 
stiffnesses  of  the  individual  facings  increases  as  the  shear  deformation  in¬ 
creases.  For  small  shear  deformations,  the  increase  is  negligible. 
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From  the  reasoning  involved  in  the  theory  leading  to  equation  (77),  consider¬ 
able  scatter  in  the  experimental  values  of  the  critical  stress  is  to  be  expected 
because  of  the  effect  of  initial  irregularities.  Similar  scatter  is  exhibited 
by  homogeneous  cylindrical  shells,  for  the  same  reason. 

The  possibility  of  failure  by  wrinkling  of  the  facings  at  a  stress  lower  than 
that  predicted  by  equation  (77)  should  be  considered. 

The  analysis  in  this  report  involves  a  number  of  approximations  and  assump¬ 
tions.  Such  procedures  are  necessary  until  a  more  rigorous  treatment  of  the 
problem  is  developed.  A  completely  rigorous  treatment  of  the  buckling  of  a 
homogeneous  cylindrical  shell  is  still  lacking,  in  spite  of  the  noteworthy 
contributions  of  von  Karman  and  Tsien. 


Tests  of  Curved  Panels 


The  large  size  of  complete  circular,  cylindrical  shells  having  realistic  fac¬ 
ing  and  core  thicknesses  and  curvatures  could  not  be  adapted  to  the  available 
testing  apparatus.  Therefore,  axial  compressive  tests  were  conducted  on 
rectangular  panels  curved  to  various  radii.  The  dimensions  of  these  panels 
were  chosen  so  that  their  widths  and  lengths  were  large  enough  to  include  at 

least  one  buckle  of  a  size  predicted  by  theory  and  a.»  )  as  shown 

in  table  2. 

It  was  then  assumed  that  the  curved  panel  would  behave  approximately  as  a 
complete  cylinder.  The  type  of  edge  support  (described  later)  was  such  as 
to  produce  no  clamping. 


Test  Specimens 


The  test  specimens  were  essentially  of  isotropic  construction  having  facings 
of  clad  24ST  aluminum  alloy  on  cores  of  either  balsa  wood,  oriented  so  that 
the  grain  direction  was  normal  to  the  facings,  or  of  corkboard  of  three  dif¬ 
ferent  densities.  Corkboard  cores  were  chosen,  because  their  low  moduli 
of  rigidity  afforded  means  of  exploring  shells  in  which  sizeable  reductions 
of  buckling  stresses,  caused  by  large  core  shear  deformations,  could  easily 
be  obtained.  These  corkboard  cores  had  shearing  moduli  of  1,500,  950,  and 
320  pounds  per  square  inch,  as  compared  to  15,000  pounds  per  square  inch 
for  the  end-grain,  balsa-wood  core. 
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Dimensions  of  the  specimens  are  given  in  table  2.  The  panel  sizes  ranged 
from  approximately  70  inches  square  to  panels  12  inches  wide  and  30  inches 
long.  Mean  radii  of  curvature  ranged  from  approximately  90  inches  to  10 
inches.  The  sandwich  constructions  had  facings  of  0.012  inch,  0.020  inch, 
or  0.032  inch  thickness  on  cores  of  approximately  1/8  inch,  1/4  inch,  or  1/2 
inch  thickness.  All  constructions  tested  had  facings  of  equal  thickness. 

The  specimens  were  manufactured  by  the  bag-molding  process.  Detailed 
description  of  techniques  and  bonding  adhesives  used  in  this  process  are 
given  in  Forest  Products  Laboratory  Report  No.  1574(3).  The  curvature 
was  attained^at  the  time  of  molding  by  using  a  steel  mold  curved  to  the  de¬ 
sired  radius.  A  strip  of  aluminum  1  inch  wide  and  0.032  inch  thick  was 
bonded  to  the  facings  at  each  end  of  the  specimen.  This  was  done  to  facili¬ 
tate  machining  of  the  specimen  ends  and  also  to  prevent  local  end  failure 
during  the  test.  The  ends  of  the  specimens  were  machined  square  and  true 
in  a  milling  machine. 


Testing 

The  vertical  edges  of  the  specimens  were  held  straight  by  loose-fitting  wood 
guides.  These  guides  were  approximately  2  inches  by  2  inches  in  cross  sec¬ 
tion  and  of  lengths  1/4  inch  shorter  than  the  test  specimen.  They  were 
grooved  in  the  lengthwise  direction  with  grooves  approximately  1/4  inch  deep 
and  wide  enough  to  allow  the  guides  to  be  slipped  onto  the  edges  of  the  test 
specimen.  No  attempt  was  made  to  clamp  the  vertical  edges  by  fitting  the 
guides  tightly. 

The  lower  ends  of  specimens  not  wider  than  30  inches  were  placed  on  a 
heavy  flat  plate,  which  was  supported  by  a  spherical  bearing  placed  on  the 
lower  head  of  a  hydraulic  testing  machine.  The  heads  of  the  testing  machine 
were  then  brought  together  until  the  specimen  just  touched  the  upper  platen 
with  no  load  indicated.  Adjustments  were  made  on  the  spherical  base  until 
no  light  could  be  seen  between  the  ends  of  the  specimen  and  the  loading  heads. 
Screw  jacks  were  then  placed  under  the  lower  loading  plate  to  prevent  tilting 
of  the  plate  while  the  load  was  being  applied  to  the  specimen.  A  single  thick¬ 
ness  of  blotting  paper  was  inserted  at  the  ends  of  the  specimen  to  help  pre¬ 
vent  local  end  failures.  The  load  was  then  applied  slowly  until  failure  oc¬ 
curred. 

Specimens  wider  than  30  inches  were  tested  between  the  heads  of  a  four- 
screw,  mechanically  operated,  testing  machine.  No  spherical  bearing  was 
used.  The  specimens  were  cut  as  true  as  possible.  If  light  could  be  seen 
between  the  ends  of  the  specimen  and  the  heads  of  the  testing  machine,  shims 
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of  paper  or  brass  were  inserted  until  the  gap  was  closed.  These  wide  speci¬ 
mens  were  also  very  long;  therefore,  small  irregularities  in  the  end  bearing 
were  absorbed  early  in  the  test  without  causing  large  variations  from  uni¬ 
formity  in  the  stresses  in  the  facings. 


Results  of  Tests 


The  facing  stresses  at  the  failing  loads  of  the  curved  panels  are  given  in 
table  2.  For  later  compari^son  with  theoretical  values,  the  parameter  N 
was  calculated  for  each  test\  specimen  by  using  the  formula 


N 


5Q 


iPfi 


4Eh 


where 

E  =  10,000,000  pounds  per  square  inch  (modulus  of  elasticity  of  facings) 

The  visible  failures  of  the  specimens  were  of  a  type  caused  by  buckling. 
Large,  thin  specimens  actually  showed  large  buckles,  which  disappeared 
after  release  of  load.  The  appearance  of  these  buckles  always  caused  a 
sudden  drop  in  the  load.  Small,  thick  specimens  showed  buckling,  followed 
immediately  by  a  crimping  appearance  at  the  edges  of  the  buckle.  This 
crimping  was  undoubtedly  due  to  shear  failure  of  the  core  caused  by  high 
stresses  induced  in  the  sandwich  by  the  buckle.  Many  of  the  thick  speci¬ 
mens  exhibited  no  visible  signs  of  buckling  but  showed  similar  crimping. 
The  rapidity  of  failure  occurring  immediately  upon  buckling  undoubtedly 
prevented  visual  observation  of  the  buckle  itself.  Similar  behavior  was  ob¬ 
served  for  cylindrical  shells  of  plywood  (11). 


Comparison  of  Theoretical  and  Experimental  Results 


The  theoretical  and  experimental  values  of  N  are  given  in  table  2.  A  com¬ 
parison  between  them  may  be  obtained  by  referring  to  figure  7  which  shows 
the  experimental  values  plotted  against  the  theoretical  values.  The  scatter 
of  points  about  a  line  representing  equality  between  experimental  and  theoreti¬ 
cal  values  shows  that  theory  and  experiment  agree  within  approximately  ^  30 
percent. 

In  view  of  the  inevitable  scatter  of  experimentally  determined  buckling 
stresses  that  is  associated  with  initial  irregularities  of  shape  and  variations 
of  material  properties,  it  is  concluded  that  the  agreement  between  results  of 


hi 


f.4 

j 

IV 


th  6  0— 133 


tests  and  of  theory  is  satisfactory.  The  scatter  is  not  as  great  for  shells  of 
sandwich  construction  as  that  observed  for  thin,  homogeneous  shells  (fig.  44, 
report  No.  1322-A  (7))  or  for  plywood  shells  (fig.  3,  Forest  Products  Labora¬ 
tory  Report  No.  1322  (18)).  This  reduction  in  scatter  may  be  attributed  to  a 
greater  total  thickness  of  shell.  Thus,  irregularities  that  depart  from  the 
true  cylindrical  surface  of  the  order  of  the  thickness  of  the  shell  are  less 
likely  to  occur  in  sandwich  shells  than  in  thin,  homogeneous  shells. 

Conclusions 

The  buckling  stress  of  long,  thin-walled,  circular  cylinders  of  sandwich  con¬ 
struction  in  axial  compression  can  be  found  with  satisfactory  accuracy  by  the 
formulas  and  curves  of  the  approximate  theoretical  analysis  of  this  report. 

Curved  panels  of  sizes  large  enough  to  include  at  least  one  ideal  buckle 
2irr  2irr 

(b  =»  ^  and  a  y  ^^  )  buckle  at  stresses  approximately  equal  to  those  of  a 

long,  complete  cylinder. 
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Notation 


a 

length  of  buckle. 

A 

f/Eb- 

Aj,  ■^2*  "^3*  "^4*  "^5 

defined  by  equations  (42)  to  (46). 

a'  a'  a'  a'  a’ 

■^1’  ^Z’  3’-  4’  5 

defined  by  equations  (52)  to  (55). 

b 

width  of  buckle. 

B 

1/E,. 

c 

thickness  of  the  core. 

Cj  mean  circumferential  stress. 


I^m 

fl  -f2 


^2-  <^3 
®xx>  ®xy’ 
ep  e^, 

^4 

^5 

^x’ 

Ea 

Eb 

fp  f, 
g 


2 

defined  by  equation  (57). 
components  of  strain. 

defined  by  equations  (26),  (27),  and  (28). 
defined  by  equation  (64). 
defined  by  equation  (66). 

Young's  moduli  of  the  facings 

Ex  (fl  +  fz) 
h 

Ey  (fl  +  fz) 
h 

thicknesses  of  the  facings. 

quantity  proportional  to  mean  radial  expansion. 
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defined  by  equation  (49). 
defined  by  equation  (50). 
parameters  introduced  in' equations  (31). 
see  equations  (74),  (75),  and  (79). 
defined  by  equations  (23),  (24),  and  (25). 
defined  by  equation  (84). 

2  IT  r /b  . 

mean  compressive  stress, 
compressive  stress  in  the  facings, 
introduced  in  equations  (31). 

radius  of  middle  surface  of  the  cylindrical  tjhell. 

defined  by  equations  (92). 

defined  by  equation  (83). 

axial  component  of  displacement. 

circumferential  component  of  displacement. 

strain  energy  of  the  core  in  the  bending  of  the 
sandwich  shell. 

strain  energy  of  the  facings  in  the  bending  of  the 
sandwich  shell. 

radial  component  of  displacement, 
extensional  strain  energy. 


flexural  strain  energy. 


w 

X^,  Xy,  etc. 

z 


Yi-  Y2.  Y3.  Y4 
6 


C 

X 

^^m 

e 

°^xy’  °^yx 
<)> 

0 


virtual  work  of  the  compressive  load. 

(Wi  +  -  W2)/abh. 

components  of  stress, 
b  /a  . 

Tr/a  . 

■tr/b  . 

defined  by  equations  (80),  (81),  (82),  and  (86). 

a  parameter  that  is  proportional  to  depth  of  a 
buckle. 

initial  value  of  6. 
mean  compressive  strain, 
coordinate  shown  in  fig.  3. 
h/r. 

-  ‘^xy 

modulus  of  rigidity  of  the  facings. 

moduli  of  rigidity  of  the  core. 

“xy  (‘l  “  ‘2' 
h 

6  r  /h. , 

6  r/h. 
o 

Poisson's  ratios  of  the  facings, 
defined  by  equation  (51). 
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Table  2. --Experimental  and  theoretical  data  for  curved  panels  of  aandvlch  conetructlon  In  axial  compreagion. 
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Figure  2.  --Nodal  lines  of  assumed  diamond 
buckling  pattern. 
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gure  4.  --N  values  for  buckling  of  cylinders  of  sandwich  constructions  in  axial  compression. 
Isotropic  sandwich  facings  of  equal  thickness,  orthotropic  core. 
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SECTION  III 


FLEXURE  Ai'il;  TORSION  OF  COMPOSITE  CYLINDERS 


By 

W.  S.  ERICKSEN 

Air  Force  Institute  of  Technology 
For  Materials  Laboratory 


Sunmary 


A  solution  to  the  problem  of  determining  the  components  of  stress 
and  displaconent  in  composite  cylinders  supported  as  cantilever 
beams  and  subjected  to  flexttral  and  torsional  loads  is  given  in 
this  rq)ort.  The  type  of  cylinder  considered  is  that  composed 
of  three  circular,  co-oxial,  layers  of  different  materials  that 
are  bonded  at  their  junctures.  Two  main  cases  are  considered; 
one  in  which  the  material  in  the  center  layer  is  cylindrically 
aeolotropic,  the  other  in  which  it  is  isotropic.  The  material 
in  the  inner  and  outer  layers  is  assumed  to  be  isotropic.  By 
tedcing  the  thickness  of  one  of  the  layers  equal  to  zero,  the 
results  are  applicable  to  a  two-layer  cylinder  and,  by  making 
the  inner  radius  zero,  they  are  applicable  to  a  solid  cylinder. 
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LIST  OF  ABBREVIATIONS 


The  following  is  a  partial  list  of  the  symbols  used  in  this  report: 

E^:  Young's  modulus  of  an  isotropic  material. 

E^,  E^,  Eg:  Young's  modiilil  of  a  cylindrically  aeolotropic 
material. 


El  =  +  ^2^2  ^3^3  intermediate  layer  is  isotropic 

=  ^3^3  Intermediate  layer  is  cylindrically 

aeolotropic . 

e^^,  .  ,  ,  ,  e^Q*  components  of  strain. 

2(1  -  -A<-) 


h- 


{cr  -  c5-  )E 
z9  zr  © 


H  2(  +  3  cr,e) 


F  =  (n  +  l)J  ^  +  (n  +  2)  „  i 

n  I  zr  ze  j 

Gjjj  Shear  modulus  of  an  isotropic  material. 


n  =  1,  j 


^r©'  <^zr  :  Shear  moduli!  of  a  cylindrically  aeolotropic 

material . 


GI  =  *^2^2  '^3^3  intermediate  layer  is  isotropic 


Hi  = 


=  Gt I,  +  G„  lo  +  if  intermediate  layer  is  cylindrically 

aeolotropic . 


+  12 


bo 
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List  of  Abbreviations  (Cont'd) 

\i~cr_cr 


rz  zr  ^6r  ^  ^ 9z  ^zr  1 


Hjj  =  (n  +  1)  (n  +  2) 


^  ■  '^ez  ^z9  1 

•-  (n  +  1)| 
1  ' 

r  J 

1  - 

"e  ,  ■  ! 

, 

E 

r 

*^9*  “^ze  .  1 

1 

I 

re 


\  -  Vi) 


1  =  (1  +A)'^  -  1 


n  - 


1,  J 


j  =-(l  +/«.)2  .  1 


k  r 


Oz 


zr 


1 

5 


^z: 


-  (n  +  l)(n  +  2) 


C5~0z 


1  /  ^  zr 

=  g-  Fjj  -  (n  +  1)  -  V-  •  -/V..  •  -g— 


=  -  -  (n  + 1)(„  *  2)  % 

z  r  ^9 

<^ze  ^  (n  +  l)^r9  .  (n  +  l)(n  +  2) 

—  -a  4 


^  *n 

:  length  of  cylinder 


m  :  index  denoting  1,  2  or  3 
M  :  applied  torsional  lacraent 
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M 


List  of  Abbreviations  (Cont'd) 


n  Index  that  ranges  over  the  values  1,  1,  J. 


applied  radial  force. 


r_  outer  radius  of  layer  ». 

B 


u ' ,  u ' ,  Ug  :  coBQionents  of  dlsplaceaent  associated  vlth  derived  stress 
^  con^nents . 

u*,  u*,  u^  :  conponents  of  rigid  body  displacement. 

Uj  =  ^  conponents  of  displacement. 


z,  r,  0  :  cylindrical  coordinates. 


^  J  ii£k£k-2 


•'r© 


% 


^ '  ‘’ez  ‘^ze 


A 


^m-l 


CT'jjj  :  Poisson's  ratio  of  an  Isotropic  materleil. 


^  :  The  value  of  Poisson  s  ratio  when  this  ratio  Is  assvned  the 

same  for  each  layer. 


<^zr'  ^re'  ' 


>  :  Poisson's  ratios  of  a  cyllndrlcally 


9z 


aeolotropic  material. 


:  A  symbol  denoting  sunmatlon  over  the  Indices  n. 
z>  ’  ’  ’  ’  ^0z  ■  Coii^>onents  of  stress. 
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INTRODUCTION 


This  report  deals  wltb  the  problem  of  determining  the  canponents  of 
stress  and  displacement  in  canposite  circular  cylinders  under  flexure  and 
torsion.  The  composite  cylinders  under  consideration  are  ccmposed  of  three 
co-axial  layers  of  different  materials.  These  layers  are  bonded  together  at 
their  junctures  so  tluat  at  the  junctures  the  components  of  stress  and  displace¬ 
ment  are  continuous.  The  materials  of  which  the  inner  and  outer  layers  of 
the  cylinders  are  coi'iposed  are  assumed  to  be  isotropic  while  the  materiel  of 
the  intermediate  layer  is  treated  in  one  case  as  isotropic  and  in  the  second 
as  cylindrically  aeolotropic.i 

The  results  developed  here  are  applicable  to  cylinders  that  are  long  in 
ccmparison  with  their  diameters.  This  restriction  Is  usual  for  solutions  of 
the  present  type.  There  are  no  further  restrictions  on  the  dimensions  of  the 
cylinder  or  on  the  relative  thicknesses  of  the  various  layers.  Thus  the 
results  cover  the  range  from  a  cable-like  cylinder  with  a  solid  core  to  a 
hollow  cylinder  with  walls  that  are  thin.  Also,  by  taking  the  thickness  of 
one  of  the  extreme  layers  as  zero,  the  results  apply  to  a  cylinder  composed 
of  two  layers. 

The  solutions  of  the  differential  equations  of  equilibrium  and  compati¬ 
bility  that  are  used  as  a  basis  for  the  analysis  have  been  previously  derived.^ 
For  the  case  of  three  isotropic  layers  results  could  have  been  obtained  by 
methods  developed  by  Muskhelishvili  and  others..  The  explicit  results  given 
by  Muskhelishvili,  those  for  a  cylinder  of  two  isotropic  layers  with  equal 
Poisson's  ratios,  have  been  compared  with  the  present  results  for  this  case 
as  a  means  of  checking  the  formulas  developed  here. 

By  way  of  arrangement,  this  report  is  divided  into  four  parts,  the  first 
of  which  contains  an  enumeration  of  the  equations  of  equilibrium  and  compati¬ 
bility,  the  boundary  conditions  and  the  conditions  imposed  upon  the  stress 
resultants  and  resultant  moments  over  sections  of  the  cylinder  perpendicular 
to  the  axis.  The  solution  of  the  problem  for  the  case  of  an  isotropic  inter¬ 
mediate  layer  is  given  in  Part  2  and  Peu:t  3  contains  the  solution  for  the 
case  of  a  cylindrically  aeolotropic  intermediate  layer.  In  Part  4  the  results 
obtained  under  the  assumption  of  equal  Poisson's  ratios  for  the  materials 
of  the  three  layers  are  given  for  the  cases  considered  in  Parts  2  and  3» 

These  results  are  considerably  simpler  than  those  for  the  general  cases  and 
may  serve  for  approximate  use. 


i  Bending  and  Torsion  of  Circular  Cylinder  Cantilever  Beams  of  Cylindrically 
.aeolotropic  t^laterial  by  W.  h,  Ericksen,  Journal  of  Applied  Mechanics,  June  195^* 

p 

=  Some  Basic  Problems  of  the  Theory  of  Elasticity  by  N.  J.  Muskhelishvili, 

P.  Noordhoff  Ltd  Groningen  Holland  1953* 
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The  writing  of  this  report  is  motivated  mainly  by  the  applicability 
of  the  results  to  sandwich  cylinders.  In  this  applicatiou  it  is  possible 
that  simplified  formulas  may  be  obtained  if  the  intermediate,  or  core 
layer,  is  thick  and  of  low  density  in  comparison  with  the  outer  layers. 
Simplifications  should  also  result  if  the  thickness  of  each  l^yer  is  small 
ccmpared  with  some  standard  radius,  say  the  radius  of  the  mid-surface  of 
the  cylinder.  These  simplifications  are  left  to  be  made,  where  possible, 
in  Individual  applications. 
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1.  FORMULATION  OF  THE  PROBLEM. 

Let  the  axes  of  reference  for  a  cylindrical  coordinate  system 
(z,  r,  0)  he  fixed  in  the  cylinder  with  the  longitudinal  axis  in  coinci¬ 
dence  with  the  axis  of  syianetry  of  the  cylinder  as  shown  in  Figure  1. 

As  indicated  in  the  figure  the  syntools  r^  and  r^  designate  the  inner  and 
outer  radii  of  the  cylinder,  respectively,  while  r^  and  rg  denote  the 
radii  at  the  junctures  of  the  layers.  The  layers  are  numbered  1,  2,  and 
3,  according  to  the  convention  that  the  one  for  which 


m-1 


£-  r  ^  r 


is  layer  m.  These  layer  numbers  are  used  as  subscripts  with  various 
symbols  to  associate  the  quantity  designated  by  the  symbol  with  a  given 
layer • 


In  the  preceding  notation  the  equations  of  equilibrium  are  given 
as  follows: 

5  ^  rrm  ,  1  P  r0m  P  ^  rzm  ^  rrm  -  00m 

— pT—  +  r  — - r —  =  ° 


1  P  ©0in  _  (P  ^  r0m  ,  6zm  ,  2 

r  Tr  Tl  F  ^  rem  =  ° 


■  ^  zrm  ^  1  V  t  ezm  ^  c-  zzm  ,  1  'V  _ 

pr  "^r  P  Q  ■*'  p  z  r  ^  rzm  ' 


t] 

CO 

[3] 


m  =  1,  2,  3 


WADD  IE  60-133 


92 


The  compatibility  conditions  are 


Let  the  cylinder  be  fixed-  at  the  base  z  =  0,  and  be  subjected  to 
a  force  P  directed  along  the  radius  0  =  0  eind  a  torsional  inoment  M 
at  z  .  Hie  conditions  loosed  upon  the  stress  resultants  and 
couples  are  then 


z 


f ' 

^m-1 


zrm 


cos 


e  -  a 


ezm 


sin  6)r  dr  d©  =  P 


m=l 


J  J'  ("^^j^sln  9  t'^-'e^mCos  e)r  dr  de  =  0 


r  ^  0 

m-1 


^ith  limitations  discussed  in  Section  4. 
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Ihe  boimdary  conditions  to  be  Iniposed  are  that  the  outer  and 

Inner  curved  surfaces  of  the  cylinder  are  free  of  stress  and  that  the 

components  of  stress  and  displacement  are  continuous  at  the  Jxmctures 

of  the  layers.  Thxis  If  the  synfcol  r,  ©)  designates  any  of  the 

stress  components  T. _ ,  %  “id  if  r,  e)  designates 

1*1111  rbn  zra  ^ 

any  one  of  the  components  of  displacement  u^  or  Uq^^  the  conditions 

Imposed  are 

"^1  V  =  0  I 


"^1  (z>  ®)  =  \  (z>  0) 

^2  ^2'  ^  ^3  ^2’ 

(z,  r,,  e)  =  0 


and 


u^  (z,  r^,  e)  =  Ug  (z,  r^,  ©) 
Ug  (z,  Tg,  8)  =  Uj  (z,  T^,  e) 


WADD  TR  60-133 


95 


2.'  RESULTS  FOR  THREE  ISOTROPIC  LAYERS. 


The  stress-strain  relations  for  the  case  of  three  Isotropic  layers 
are  written 


zzm  E^ 
1 


rrm 


'0©m  ~ 


'r0m 


e  = 

ezm  2Gjjj  ^  92m 


e  =  ^1. 

zrm  2G  zrm 
m 


- 

-  1 

■  a —  i- 

zzm 

m  rm 

m  0 

-(T 

+  "2-  ^ 

-  cT-  a 

zzm 

rrm 

m 

-  ^ 

-  ^  1 

. 

m 

zzm 

m 

rrm 

(T) 

m 

1 

‘■r9m 

06m 


M 


m  =  1,  2,  3. 

The  components  of  stress  given  by  the  following  expressions  satisfy 
the  equations  of  equilibriimi  and  j^3  j/  and  by  the  tise  of 

relations  they  also  satisfy  the  compatibility  equations  M  to 

[0  inclusive . 

Tw] 

r^o] 

H 


^  -(8crA+C)  -  z)  cos  0 

^ zzm  mm  m 

=  (2A  r  -  2  (^  -  z)  cos  9  • 

^rrm  m  ^3 


1 


2B 


96m 


(6Ajjjr  +  -  z)  cos  0 

r'^ 


^rem  '  (2A^r  -  -  z)  sin  9  - 
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T.  =j+^+A  jscr-  +  — (1  Ir^ 

znn  j.2  m  |  m  m  J 

C 

m  ,  _  .2 

+  TO-  (3E_  -  2^5"  G  )r  cos  8  •  ■  . 

^  OE  m  mm 

m  -J 

7  .  fij  .  ^_3a 

L  m  2  ®3  E  m'' 

ezm  L  r  I  m 

-3C_  E  2" 

(  -3  -  2  (f  G  )r  sin  e  +  D  G  r  .  . 

Se^  T  m  m''  J  ^ 


[.J 


H 


lliese  particular  components  of  stress  prove  to  be  satisfactory  for 
the  present  loading  conditions.  They  may  be  derived  by  the  methods  of 
Section  3  of  Reference  1.  However,  the  fact  that  they  satisfy  the  con¬ 
ditions  of  equilibrium  and  compatibility  is  readily  verified  by  direct 
substitution. 

When  the  components  of  stress  are  of  the  preceding  forms  the  com¬ 
ponents  of  displacement  are  given  by  the  expressions  that  follow.  In 
these  expressions  each  component  is  separated  into  two  parts,  the  first 
of  which  is  primed  and  designates  the  component  directly  associated 
with  the  given  stress  components,  while  the  second  part  is  starred  and 


represents  a  rigid  body  displaceiaent .  Uaus 


where 


u  =  u  +  u 

y.TH  y.TTt  7.Tn 


*  * 

U  =  U  +  u  .  ,  . 

m  nn  rm 


U  =  U-  +  Ucw  •  •  • 

0JQ  9ni  wia 


•  r 


U] 

[-] 
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V 


(1  - 


c  eH  1 

ni  m  o.^  21  I 

-  -8^  L 


B„(l  +^„)  ^  ]  1-] 

E  0  \  T 

m  m  J  _ 


cos  9 


[.] 


=  -rf (/-  t.  ”  (1  -  1*0(1  rV-  z) 

rm  oE  E  m  m  o-p 

_  m  m 


B  (1  +<r  )  (^-2) 

SI  n 


cos  9 


t] 


sin  0  +  D  rz 
m 


=  -  ._^  r  cos  0  -  cC  r  sin  0  +  _/ 
m  o'm 


=  (c<„  =03  e  +  (y  +  5  z)sin  e 

TTH  m  /"^si  IS  SI 


^  =  -(‘^m  +/-4^)  sin  ©  +  (a  ^  +^^2)  cos  9  +  ^  ^r. 


[-] 
I  “I 
I_32j 

[33  _ 


Tile  Constanta  appearing  in  expressions 
in  [28  J  ,  CsJ  and  w  are  dete: 


^19j  through  1^24^ 


are  determined  explicitly  by  conditions 


M  .  [-]  -  M  and  as  follows; 

Ai  =  — ^  [(‘^2  '^3 


DO 

DO 
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V 


Having  made  these  determinations  the  foUovlng  e:{pres8lon8  may  he 
evaluated: 
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Conditions  1^  [13] ,,  [11.]  and  have  not  been  lapsed  In 

the  determination  of  the  components  of  stress  and  displacement.  Among 
these  r-]  Is  implied  by  j_  3  J  >  [^]  and  by  H  applied  to  the 
components  ^rm^  'diile  the  remaining  conditions, 

[[15J  ,  are  satisfied  Identically  by  the  forms  ^1^  ,  |2^  and  . 

The  constants  appearing  In  expressions  [31]  ,  [32]  »4  [33]  , 

depend  upon  the  method  of  fixing  the  supported  end  of  the  cylinder 
for  their  final  determination.  They  are  interrelated  by  the  following 
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RESULTS  FOR  ISOTROPIC  INNER  AND  OUTER  LAYERS  WITH  CYLINDRICALLY 
AEOLOTROPIC  INTERttePIATE  LAYE^ 

For  the  Isotropic  layers,  1  and  3>  the  stress-strain  relations  are 
given  by  M  ,  and  expressions  H  to  H  ,  inclusive,  m  =  1,  3, 
are  adopted  as  the  forms  of  the  components  of  stress  and  displacement. 


The  stress-strain  relations  and  the  expressions  for  the  components 
of  stress  and  displacement  for  the  intermediate  layer,  2,  are  taken 
from  Reference  1.  Since  the  notation  for  the  elastic  coefficients 
differs  in  form  from  those  for  an  isotropic  material,  the  subscript  2  is 
not  used  with  these  coefficients  and  the  stress-strain  relations  are 
written 

e 


^rz  ^ 

'z22  ■  zz2  Ep  '■  rr2  Kg  ©©2 


"rr2 


^  zrQr  .  ^6r  7 

~  t:'  '  ^  ^  rr?  _ 


cr- 


z9 

®9e2  "  ~K~  222 


Ep  ‘ rr2  Eg 


=6 

©92 


©92 


'ez2  ■  5g 


re 


r©2 


-©z2 


'2r2 


2Gg^  ‘•©22 


1  n 

zr2 • 
zr 


w 


Hie  forms  of  the  components  of  stress  are  taken  as  follows: 

n 


1  =  V  S  F  r'^  z)  cos  e 

‘-ZZ2  n  n  '' 


w 
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In  these  expressions  the  symbols  for  n  =  1,  1  and  j,  Tj^,  T  and 
Dg  are  constants  that  are  to  be  determined  by  the  boundary  and  end 
conditions.  The  remaining  symbols,  together  with  i,  j,  and  k,  are 
defined  under  the  List  of  Abbreviations. 


The  components  of  displacement  for  the  intermediate  layer  are  each 


separated  into  two  parts  as  in  Formulas  [25]  [se]  ,  and  w  • 

In  these  formulas  the  starred  components  are  given  by  expressions 
[31]  ,  [32]  ,  and  [33]  with  m  =  2.  The  primed  components  are 


given  as  follows; 
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z2 


zY 


It. 


s 

n  ^n 


k  „  -k  - -  S  Kn  +  2)F  +  2Gq  hI: 

r  -  T  ,  r  +  k  \  n  L: _ ^  n  9z  pj 

^  (n  +  2)  ^{n  +  2)2  .  k21 


os  6  [-] 


Is^KiiC^ -  z)3 


SnK2n 


n  +  1 


cos  0 


w 


+D^rz  . 


H 


The  constants  appearing  in  the  expressions  for  the  components  of 
stress  and  displacement  for  the  three  layers  are  determined  as  follows 
ly  the  use  of  H  ’  H  ■  0"]  ■  ["]  .  These  directly 
associated  with  the  bending  of  the  inner  and  outer  layers  are  related 


to  the  coefficient  for  the  intermediate  layer  by  the  equations 


=  TT 


Cl 

^  "  ^1^11 


[”] 

H 

H 

Qa] 

H 
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With  the  use  of  these  relations  the  constants  for  n  =  1  and  i  and 
j,  are  evaluated  as  the  solution  of  the  following  system  of  equations: 


.  ^3*3X3,  *  03X, .  C3X3  = _  p.  [so] 

1  +^3)  I  (3  -  4  cr^)r\  -  r\| 


^(1  +^)  J  (3-4 


^(1  +  J  (3  -  4  cT  )  r  +rM 


3'  2  "  3^- 


Y  _iiL_|K  t 

/a.  n  +  1  2n 

<■ 


K  X  iL±i  pf' 
3”  “re  " 


“1  =  “2  =  “3  = 


2^0  1  +GI  +G^  1*1 
1  1  3  3  ez  2J 


The  remaining  constants  in  the  expression  for  the  components  of 
stress  and  displacement,  exclusive  of  those  that  determine  the 
rigid  body  displacements,  are  given  in  terms  of  the  preceding  con- 


stanL-,,  as  follows: 

T 1  ^>3^1^  -  Y  '^i  "'i 


-k-1  2 
1  ^2 
k-1  _k-l 


-^1  '^3"'l'^  ^  -'^3^1  ^2 
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where 


.  ^2(2-m)  kC  2(2-m)  ] 

^  “  Go.  1  ®  1 

©z  \ 

(  ^2(2-m))  2(2-m) 


A  •♦(2-m) 


m  G 


A-  (n  +  2)2  -  k^  1 


(n  +  2)g._^^  ^ 
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4.  RESULTS  FOR  EQUAL  POISSON'S  RATIOS. 

When  the  Poisson's  ratios  of  the  three  layers  are  eqiial,  the 
expressions  for  the  components  of  stress  and  displacement  are  consider¬ 
ably  simplified.  These  simplified  results  are  given  in  this  section 
both  for  the  case  of  an  isotropic  intermediate  layer  and  that  of  a 
cylindrically  aeolotropic  intermediate  layer . 

Since  the  effects  of  varying  the  Poisson's  ratios  are  expected 
to  be  small,  at  least  vd.thln  the  ranges  of  these  ratios  found  in  commonly 
used  engineering  materiails,  the  formulas  of  the  present  section  are  con¬ 
sidered  to  be  applicable  for  obtaining  approximate  results  for  the 
general  cases  of  the  two  preceding  sections. 


w 

[loo] 


Cylinder  with  an  Isotropic  Intermediate  Layer  with 


When  the  Poisson's  ratios  are  equal  formulas  |^3 

4^  through 

39"] 

reduce  to 

1 

^  1 

A  =  A  =  A  =  0, 

1 

123 

1 

fi 

and  B  ■ =  B  =  B  =  0 . 

123 
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Also,  from  w  and  W- 


where 


El  =  E^I^  +  ^  E3I3. 


where  r~  “I 

El  =  E^I^  +  E2I2  +  E3I3.  I  loy 

The  components  of  stress  with  from  >  then  reduce  to  the  forms 

T  .z.  =  ¥ '  H 

^  =^.  =1  =  0.  TlO^ 

^  rrm  00m  r0m  L  J 


-E  P 

=  r  C^-  z)co3  e. 

2zm 


rrm  ^  66in  r9m 


-p  r  '  r^ 

"1  =  ET  Km  +  -  +  ^  (3E  -  2  cro^)  cos  6, 

zrm  Ur  -J 


sin  © 


4 

2GI 


GX  —  G^Xg  4 


and  where 


where  r  ~l 

GX  =  G3^Xj_  4  Ggig  +  G3X3  I  lOd 

and  where  •  J  EX  PI 

H 

The  constants  Jjj^  and  are  determined  hy  |^46  to  1^51"^  ,  inclusive 

with  .p  r~  ~I 

^  KI2J 


^m  ■  P  ■ 
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_  ^  {  2(2-m)  2(2-m) 

I  m  3“)  J) 


^(2-ra)  G„  1  2(2-m)M  ■  T  ~1 

*"4  )J’  H 

and  with  ^  and_rL  as  given  by  55  and  56  ,  respectively, 

in  m  1 _  J  L—  J 

The  components  of  displacement  associated  with  the  preceding  stress 
components  are  given  by  the  expressions 


u'  =  P  2cr_^ 

zm  El  I  2  o  V  a 


ti-  " 


U  -  ~  r  ( ,  _  -  z ) 

m  El  L  6  ^  ~  J 


cos  9 


em  El  L 


ilL  +  a-  ^  ('<-0  si 


sin  S. 


lll^ 

H 


Formulas  W  ’  [^il  and  1^33^  for  the  rigid  body  displacements 


remain  as  given. 


As  an  example  of  the  manner  in  which  the  rigid  body  displacements 


may  be  fixed,  consider  the  conditions 


=  0  for  z  =  0,  r  =  r^,  0  ^  2 

=  0  for  z  =  0,  r  =  r^^,  9  =  f 

Up  =  0  for  z  =0,  r  =  ri,0  =  ± 
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Ihese  conditions  fix  the  plane  of  the  circle  z  =  0,  r  =  and  the 
two  points  on  the  circle  at  6  =  +  •  With  these  conditions,  the 

constants  in  expressions  [31]  ,  [3=]  »d  [33]  take  on  the  values 

u.--i\4-^]  r- 


[U9] 

i.‘  [12^ 

It  is  observed  that  expressions  [us], [lie]  .[us]  ,  [119] 

end  H  are  independent  of  the  index  m.  The  components  u^.  and  U0 
for  a  given  layer  do  not  then  depend  in  a  special  way  upon  the  elastic 
coefficients  for  the  material  of  the  layer  but  rather  upon  the  elastic 


properties  of  the  composite  cylinder.  Therefore,  the  deflection  of 


the  cylinder, 


u  =  Uj.  cos  0  -  Uq  sin  9 


H 


is  given  by  the  formula 


t  -  z)~ 


+  ^  {■'C  -  z )  cos  2 


Mzr  sin  9 
2GI 


H 


with  c>(  and  o  given  by 
1  m  m 


and  119  ,  respectively.  With 


the  use  of  the  latter  formulas  the  deflection  at  the  end  z  =  ^  is 


obtained  in  the  form 
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With  the  substitution  of  |^125  j  and 
inclusive,  it  is  found  that 


M 


into 


75 


to 


79 


-  Bi  =  B3  =  0 


and 


El  -  Eg 


p 

El 


H 

[130] 


The  components  of  stress  for  layers  1  and  3  are  of  the  forms 


105 


to 


108 


inclusive,  with  J  and  L  defined  below  and,  since 
m  ai 


the  constants  are  now  determined  by 


[33]  , 

GI  =  +  Gg^Ia  +  °3^3- 

For  the  intermediate  layer. 


with 


'-zr2  ~  El 


Ez?  / 

- —  r( 

(' 

x:  -  2 

)co6 

^  zz2 

El 

1  = 

II 

"I 

=  0 

rr2 

992 

‘■r92 

+  '^Ik 

-k-1 
r  _ 

(3S2  - 

2cr 

9z ' 


9  -  k‘ 


cos  9 


[131J 

_130 

133J 


rl 


-p 

■ez2  ■  ET 


k-1  -k-1  G  )r2 


9  -  h. 


2GI 


where 


T^EI 


[13] 
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V 


T' 


the  determi-nations  of  these  tvo  constants 


.(313  -  2  ) 

m _ m 

12 


I  _  2(2-m)  .  ,  ^ 

ix.  f  +  3  ^ 


2(2.ra) 

m 


[137] 


vrlth  Tj,  and  T_j,  determined  from  H  and  w  ,  respectively. 


For 


(  3S2  -  2  G92  ) 

+  9  -  k2 


2(2-m)  G  )  2(2-m) 

“  3G  !  m 

zr 


V 


[138] 


and  the  quantities  ~r\_  and  ^  remain  as  given  by  86  and  87  , 

*  t 

respectively.  Having  evaluated  H  and  T'  ,  the  constants  J  and  L 

K  -k  mm 

are  determined  from  the  equations 


1  L,  k+1  -k+l  (3E  -  2i5~ Gft,)  r  [  ,  2  2 

T  -  — n -  <T'r,  +  T'  r,  _  __2 _ 

“^1  ■  r^  -  r  J  k  1  -k  1  ...  2 


1  Ot 

<2  '  k 

L'  =  -J  r  -Nr 
1  10  10 


9  -  k" 


'-KjCr^  +  rg) 


■’a  =  ;r7f  1  X 

7  2  2, 

-  N  (r  +  r  ) 

3  2  3 


'  '2  '  k 

L  =  -J  r  -Nr 
3  3  3  3  3 


T.  (3E,  -  2crGe,)  ^2 

"  -k  *2 - — 72 - 


9  -  k" 


where 


[lUo] 
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■Hie  components  of  displacement  for  layers  1  and  3  are  obtained 
from  [31J  ,  [32]  ,  [33]  ,  [11^  ,  [11^  and 

for  the  intermediate  layer  are  obtained  from  to  inclusive 


and  the  expressions 


u-  =  L.  +] 

22  El  L  2 


^ 


3(9  -  k^)G- 


^^2=  fc  [■  -  "ll  ®  ^ 


[l4^ 

[l4^ 


For  the  case  under  consideration,  conditions 


lead  to  the 


evaluations 


?  p 

.  m  ~  !!rL 2 


3  crrV' 


r  3e,  -  ag-Ce^ 
3(9  - 


/  =  S  =  S  -  €  =  0. 

m  m  m  m 


The  deflection  of  the  cylinder  is  obtained  from 


use  of 


l44  and  145 


H 


with  the 


At  z  =  i  the  deflection  is  given  by 
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1  -  p_ 

\il , 

3E^-2a-=^; 

|M 

Im 

1 

II 

t) 

33' 

(9  - 

/fl  j  I  V  k+1  -k-1  /\  -k+l  ^  \  /T_  > 

£  fl  1-2  +'^3r.i  1-2  )  -  r2  (A  3^+ 

„  ,  r,  >  k-1  -k-1  _n.  A 


-k-1  k-1 

^2 


r.  ) 


r  sin  6 


2GI 


where 


T' 

m 


T  El 
m 

P 


with  T  defined  by 
m 


and  El  by 
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SiETION  IV 


BUCKLING  OF  SANDWICH  CYLINDERS  IN  TORSION-’ - 


By 

H.  W.  MARCH,  Mathematician 
— .  and 

EDWARD  W.  KUENZI,  En’gineer 


Forest  Products  Laboratory,  —  Forest  Service 
U.  S.  Department  of  Agriculture 


Summary 


This  report  presents  a  mathematical  analysis  leading  to  the  critical  stresses 
that  determine  the  buckling  of  cylinders  of  sandwich  construction  in  torsion. 
The  analysis  is  complete  for  cylinders  of  finite  length  having  orthotropic 
facings  and  orthotropic  cores.  Formulas  are  given  for  sandwich  cylinders 
of  finite  and  infinite  length.  Curves  are  given  for  buckling  coefficients  for 
cylinders  with  isotropic  facings  and  orthotropic  cores,  over  a  wide  range  of 
cylinder  sizes  and  material  properties. 


Introduction 


The  purpose  of  this  report  is  to  develop  formulas  by  means  of  which  the 
buckling  stress  can  be  calculated  for  a  sandwich  cylinder  subjected  to  torsion. 


“This  report  is  one  of  a  series  (ANC-23,  Item  57-3)  originally  prepared  and 
distributed  in  June  1953  by  the  Forest  Products  Laboratory  under  U.  S. 
Navy  Bureau  of  Aeronautics  Order  No.  NAer  01336  and  U.  S.  Air  Force 
No.  USAF-33{038)-5  1-4326-E,  Amend.  2(53-131).  This  revision  was 
prepared  under  U.  S.  Navy  Bureau  of  Aeronautics  No.  01778  and  01768 
and  U..  S.  Air  Force  No.  33(616)56-9,  Amend.  A3(57-398).  Further  re¬ 
vision  may  be  made  as  additional  data  become  available. 

2 

“Revised  December  1957  by  Charles  B.  Norris,  Engineer. 

3 

Maintained  at  Madison,  Wis.,  in  cooperation  with  the  University  of  Wisconsin. 
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Methods  have  previously  been  developed  to  determine  the  buckling  stress  of 
homogeneous  isotropic  cylinders  (l)^and  of  plywood  cylinders  (7).  The 
analysis  for  plywood  cylinders  (7)  based  on  an  energy  method,  will  be  used 
in  the  present  report,  with  suitable  modifications  to  take  into  account  the 
effect  of  shear  deformation  in  the  core  of  the  sandwich  shell. 

The  core  and  facings  of  the  cylindrical  sandwich  shell  will  be' taken  to  be  or¬ 
thotropic  with  two  of  their  natural  axes  parallel  to  the  axial  and  circumfer¬ 
ential  directions,  respectively. 

The  facings,  which  may  be  of  either  equal  or  unequal,  thickness,  are  assumed 
to  be  thin,  but  their  flexural  rigidities  are  not  neglected,  as  they  may  be  of 
importance  in  certain  cases.  All  stress  components  in  the  core  are  neglected 
except  the  transverse  shear  components,  and  these  are  taken  to  be  constant 
across  the  thickness  of  the  core.  An  energy  method  will  be  employed  to  de¬ 
termine  the  buckling  stress.  In  contrast  with  the  treatment  of  sandwich 
cylinders  in  axial  compression  (6)  in  which  a  large-deflection  theory  was  used 
the  displacements  will  here  be  assumed  to  be  such  that  a  small-deflection 
theory  is  applicable.  The  use  of  a  small-deflection  theory  instead  of  a  large- 
deflection  theory  constitutes  the  chief  difference  in  the  treatment  of  the  two 
problems.  The  consequence  is  that  in  applying  an  energy  method  a  large  part 
of  the  analysis  for  both  problems  is  the  same  in  its  main  features.  In  order 
to  avoid  repetition  of  the  details  of  exactly  the  same  steps  in  the  present  re¬ 
port  that  are  to  be  found  in  the  previous  report  on  sandwich  cylinders  in  axial 
compression  (6),  reference  will  frequently  be  made  to  that  report.  Sufficient 
details  will  be  presented  here  to  show  the  principal  steps  in  the  analysis  and 
the  differences  that  result  from  using  a  small -deflection  theory  and  a  differ¬ 
ent  form  of  the  buckled  surface  from  that  chosen  for  cylinders  in  axial  com¬ 
pression. 

Large-deflection  theories  for  cylindrical  shells  subjected  to  torsion  will 
yield  only  slightly  lower  critical  shear  stresses  than -small -deflection  theories 
if  the  imperfections  in  the  cylindrical  surface  are  small  compared  to  the  thick 
ness  of  the  shell  (8,  12).  Cylinders  of  sandwich  construction  usually  exhibit 
exceedingly  small  imperfections;  thus,  the  small-deflection  theory  given  here 
seems  adequate. 

Choice  of  Axes  and  Notation 


The  choice  of  axes  is  shown  in  figure  1.  The  coordinate  y  is  measured  along 
the  circumference.  The  components  of  displacement  are  u,  v,  and  w,  re¬ 
spectively,  where  w  is  positive  inward.  Love's  notation  (4)  will  be  used  for 


—Underlined  numbers  in  parentheses  refer  to  numbered  references  under 


Literature  Cited  at  end  of  report. 
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the  components  of  strain  e^^,  etc.  ,  and  for  the  components  of  stress 

X^,  X^,  etc.  The  thickness  of  the  core  will  be  denoted  by  c  and  that  of  each 
facing  by  fj,  and  f^,  respectively,  as  in  figure  2. 

Form  of  the  Buckled  Surface 


As  an  energy  method  is  to  be  used  to  determine  the  buckling  stress,  it  is 
necessary  to  assume  a  suitable  form  for  the  buckled  surface.  The  follow¬ 
ing  form  will  be  chosen: 


—  =  g  sin  0'(y  -  sin  /3x  (1) 

where 


n  =  the  number  of  buckles  in  the  circumferential  direction  and  b  =  the  length 
of  the  cylinder.  The  first  trigonometric  factor  in  equation  (1)  is  suggested 
by  the  form  that  has  been  used  to  represent  the  buckled  surface  of  an  infinite 
plane  strip  (10)  under  uniform  shear.  The  factor  sin  |3x  is  introduced  to  make 
the  deflection  vanish  at  the  ends  of  the  cylindrical  shell.  Further  conditions 
at  the  ends  are  disregarded.  Except  for  very  short  cylinders,  these  condi¬ 
tions  are  not  important. 

It  is  convenient  to  write  equation  (1)  of  the  buckled  surface  in  the  form 

?  ~  §  l^cos  (oy  -  6x)  -  cos  (oy  -  e  x)J  (3) 

where 

6  =  Qy+/3,  (  =  ay  -  p  (4) 


Extensional  Strain  and  Stress  Components 


Expressions  will  now  be  written  for  the  extensional  strains  uniform  across 
the  thickness  of  the  cylindrical  shell  and  for  the  corresponding  mean  mem¬ 
brane  stress  components.  On  these  will  be  superposed  a  system  of.  flexural 
strains,  and  the  energy  of  deformation  associated  with  each  system  of  strains 
will  be  found. 
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The  extensional  strain  components  are  expressed  by  the  equations: 

_  ^ 

®xx  “  3x 

9v  w 

®yy  “  iy  ■  ~ 

_  0u  8v  - 

®xy  “  ‘  8x 


(5) 


In  each  facing  the  corresponding  stress  components  are: 


A. 

“  ~ 

(^XX  + 

^yx  ^yy) 

Y 

(e  + 

yy 

cr  e  ) 

xy  XX 

y 

X. 

x„ 

=  u. 

e 

y 

'^xy 

xy 

(6) 


where  E  and  E  are  Young's  moduli,  p  is  the  modulus  of  rigidity  for  shear- 

^  _ y 

ing  strains  referred  to  the  x  and  y  directions,  cr  -  and  cr  are  Poisson's  ratios 

-  -  yx 

and  X  =  1  -  (T^y.  cTyj^.  All  of  these  quantities  are  elastic  properties  of  the  fac¬ 
ings.  Because  the  stress  components  X^,  Yy,  and  X  are  neglected  in  the 
core,  the  mean  membrane  stress  components  for  the  shell  are: 


Ea  . 

"  \ 

^yx  ^yy 

Eb  , 

(e  +0*  e 

'  yy  xy  xx 

\ 

-  U 

e 

xy 

) 


(7) 


where 


E  =  E  (1  - 
a  x  '  h 


Eb  =  Ey  (1  -J). 


^^m  =  l^xy  -f) 


(8) 


and 


c  +  fi  +  f2 


(9) 


By  using  the  relation 


E.,  cr  =  E,  cr 
a  yx  b  xy 
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where— 

1  1  1 

A  =-^.  B=-=-,  C=- -  - 

%  ►"m 

For  comparison  with  Report  No.  1529  (7)  set  6  =  0  in  that  report  and  inter¬ 
change  A  and  B. 

The  expression  for  w  from  equation  (3)  is  substituted  in  equation  (14).  A 
solution  of  the  resulting  equation  is  found  to  be 


F  =  aj  cos  (oy  -  6x)  -i  a.^  cos  (oy  -  ex) 


where 


,  ,  _  IL_ 

1  "  2Kj’  ^2  -  •  2K^ 


=  A6^  +  +  Ba'^ 

4  ?  2  4 

K2  =  Ac  +  Ca  «  ^  +  Ba 


In  accordance  with  equations  (12)  the  stress  components  can  be  found  from 
the  stress  function  F. 


Extensional  Strain  Energy 


The  strain  energy  associated  with  the  membrane  strains  (5)  is  expressed 
by  the  integral 


b  2i 


%  =  -y  /  f  IXe  +  Y  e  +  X  e  )  dydx 
1  Z  J  J  '  X  XX  y  yy  y  xy'  ' 


With  the  aid  of  equations  (10),  this  integral  becomes 


"The  definitions  of  A  and  B  of  reference  (7)  are  interchanged  here  to  agree 
with  the  definitions  used  in  reference  (6).  It  can  be  shown  that  H  of 
reference  (7)  is  replaced  here  by 
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Zn-r 


W 


.=t// 


,  —  2 

(BX^ 


—  2 
+  AYy 


■o'xy - 


•)  dydx 


By  substituting  the  values  of  Y^,  and  X^,  that  are  obtained  from  equations 
(12)  and  (15)  and  performing  the  integrations,  it  is  found  that 


W,  = 


irrbh 


B  (a ,  +  3,  T  ) 


A'  2.4^ 
A^a,  0  +  a 


2  4v 
>  «  ) 


C(aj^6^ 


4.  =  2  2»  2 
+  a^  €  )  o 

By  using  equations  (16),  (17),  and  (18),  the  expression  for  Wj 
the  simpler  form: 


W, 


irrbhg 

8 


i_+  f_ 

K,  K, 


(21) 

is  reduced  to 

(22) 


Flexural  Energy  of  the  Shell 


As  in  reference  (6),  the  following  expressions  will  be  used  for  the  curva¬ 
tures  and  unit  twist: 


2  2  2 

3^w  a^w  a^w 

ax'^  ay^  axay 


These  expressions  are  exactly  those  used  in  calculating  the  flexural  energy 
of  a  flat  sandwich  plate.  The  approximate  flexural  energy  of  such  a  plate 
was  found  in  references  (5)  and  (2)  by  using  the  tilting  method  of  Williams, 
Leggett,  and  Hopkins  (11).  In  this  method  the  transverse  shear  components 
are  taken  to  be  constant  across  the  thickness  of  the  core.  The  following  ex¬ 
pressions  are  chosen  for  the  components  of  the  displacement  in  the  core: 


Uc  =  -k(C-q)|^ 
Vc  =  -k’  (C-q’)|7 


(23) 


Wj,  =  w  (x,  y) 

The  coordinates^  and  the  distance  q  are  shown  in  figure  2.  Thus  ^  =  q  de¬ 
notes  the  surface  in  which  the  components  of  displacement  in  the  x  direction 
vanish  and  the  parameter  k  describes  the  inclination  in  the  x  direction  to  the 
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normal  of  the  deformed  surface  of  the  panel,  of  lines  that  were  normal  to 
the  undeformed  surface  before  deformation  took  place.  Similarly  q'  and  k' 
are  related  to  the  displacements  in  the  y  direction. 

In  obtaining  the  expression  for  the  flexural  energy  of  deformation,  which 
will  be  denoted  by  W2,  it  is  assumed  that  all  stress  components  in  the  core 
may  be  neglected  except  the  transverse  shear  components.  The  facings  are 
treated  as  thin  plates. 


The  details  of  the  calculation  of  and  its  minimum  value  will  be  found  in 
reference  (6),  which  in  turn  depends  for  some  of  its  details  on  reference  (2). 
The  following  expression  is  found  for  the  minimum  value  of  with  respect 
to  the  quantities  k,  k' ,  q,  and  q' : 


W  = 

2  8 


< 


Aj  +  +  A3’  +  (Aj’A^’  -  A^ 


VA4>  Ag'yi 


Aj'4>  A3’4.  4>^  -  A  ’■ 

1  +  - + -  + - 


4  5 


^4'^' 


+  If  (A^  +  2A2'  +  A3') 


(24) 


where  the  coefficient  of  the  expression  in  brackets  is  different  from  that  in 
reference  (6)  because  the  derivatives  of  a  different  function  w  appear  in  the 
integrals  to  be  evaluated,  and  where 


^1^2 


(^1  +  ^2) 

,3 


(  h  +  c)" 


If  = 


fl3  +  f- 
12 


(25) 

(26) 


4>  = 


II  +^2 


H  =  Aj  ,+  2A^  +  A3 


(6^^  +  €  '^)  +  (6^  +  * 


(27) 

(28) 
(29) 
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(30) 


I  ~  ^  °’yx  ^^^xy  J 

>'4 


+  c") 


^  ■*■  ^l^xy  +  «  ^) 


A4.  =  ^'y  (6^  +  €  ^) 


Cx 


A5  =  2|i  a-' 


(31) 

(32) 

(33) 


In  equations  (32)  and  (33)  and  ji’  ^  denote  the  transverse  moduli  of 

rigidity  of  the  core.  The  elastic  properties  of  the  facings  are  denoted  by  un¬ 
primed  letters  in  equations  (29),  (30),  and  (31). 

It  follows  from  equations  (28),  (29),  (30),  and  (31)  that 


H 


(6^  +  «  ^)  +  2Ey  +  2E  (6^  +  c  ^) 


(34) 


where 


L  =  Ev  O'  +  2\u 

X  yx  '^xy 


(35) 


Work  Done  by  the  Applied  Couple 


The  work  done  by  the  applied  couple  in  producing  the  deformation  of  the 
cylindrical  surface  described  by  equation  (1)  is  given  by  the  following  integral: 

W3  =  .rh^  dydx  (36) 

o  o 


where  t  is  the  mean  shearing  stress  induced  in  the  cylindrical  shell  by  the 
couples  applied  at  its  ends.  The  corresponding  stress  in  the  facings  is  given 
by: 


T 


WADD  m  60-133  127 


Substitution  of  w  from  equation  (1)  in  equation  (36)  yields 

W  - 

3  ~  2 


(37) 


Critical  Stress 


The  critical  shearing  stress  is  determined  from  the  equation 
W3  =  Wj  +  'W2. 

On  substituting  the  expressions  for  and  from  equations  (22), 

(24),  and  (37),  and  solving  the  resulting  equation  for  t,  it  is  found  that 


+  2A2’  +  A3'  +  (Aj’ A3’  -  A^'^)  (-^  +-^ 

A4  A3 


^'4.  A3’<}.  (j>^  (AjA^  -  A^’  ) 


A4  Ag 


+  If  (Af  +  2A2  +  A3  ) 


(38) 


The  mean  critical  shearing  stress  is  the  least  value  of  x  that  satisfies  equa¬ 
tion  (38).  The  justification  for  the  previous  calculation  of  the  minimum  of 
W2  with  respect  to  k,  k  ,  q,  and  q  is  found  in  the  fact  that  these  quantities 

occur  only  in  the  expression  for  ^2' 

It  will  be  convenient  to  introduce  coefficients  K  and  by  the  equations 


h  h 

T  =  KE^  -  and  Xf  =  Kf  E^  - 


(39) 


where  Xf  is  the  shear  stress  in  the  facings  at  the  critical  load  and 
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4',  = 


Jp‘ 


1  ^  2 
4n  y 


(v  +  p) 


(v  +  p)4  +  (-^  -  20-  )(y  +  p)^  +  1 

Ey  <‘„y 


(v  -  p) 


^  (y  -  P)^  +  -  2a  Hy  -  p)2  +  1 


(44) 


For  the  second  term  of  equation  (40)  the  following  additional  parameters  are 
convenient: 


4>Ex 

^  Ik-p'r  rh 


and  S  = - r 


<j>E, 


y  Xp  ,  rh 


(45) 


Substitution  of  the  expression  (^iven  by  equations  (43)  into  equations  (29),  (30), 
and  (31)  yield: 


^4-  It  *4-  II  *4 

Aj  =  E,  A,  ,  A^  .  E^  A,  ;  A,  = 


(46) 


where 


A,"=f  (Y"+6y2p2,p4),_^(^2,p2) 


(47) 


A  ''  T  ^y  1  ->  ^’^y  /  2  ^  Z\ 

A3  =  2  +  2  (y  +  p  ) 

Using  equations  (32),  (33),  (42),  (43),  and  (45)  the  following  relations  are 
found: 


4>  _  TT^X  Sx 

A4’  ■ 

(t>  _  Sy 

A  '  T  2-  4^ 
A5  Jp  O' 


(48) 
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where 


2  ? 

A,  = 


Using  equations  (41)  and  (44)  it  is  found  that 


'772  2 

h  Jp  a 

Substitution  of  the  expressions  given  by  equations  (46),  (48),  and  (50)  into  the 

c 

second  part  of  equation  (40)  and  dividing  by  1  -  yields: 

I  ''  '•  "  II  II  lip  S 

—  Aj  +  2A2  +  Aj  +  (AjAj  -  A^  ) — -  (— rr  +  Sy) 

^  Jp  A^ 

(l-f)4JpS7  ,2^  Aj  ^4^2  a';  a’:  -a”" 

l^H<^x44SyA3)4^S,S  1  3  2 

/  Jp'"  A4  y  j2p4  ^ 


+  (Aj  +  2A2  +  A^) 

h 


Thus  equation  (40)  can  be  written 

Kf  =  4^1  +  4^2 


which  is  substituted  in  equation  (39)  to  obtain  the  shear  stress  in  the  facings 
(t£)  at  the  critical  load.  4"  j  is  a  function  of  J.,  y,  p,  and  the  elastic  proper- 


c 

ties  of  the  facings  and  4^2  i®  ^  function  of  j^,  J,  — ,  - ,  y.  p,  me: 

elastic  properties  of  the  facings.  The  values  of  y  and  p  are  chosen  so  as  to 
obtain  minimum  values  of  K^, 

If  the  facings  are  of  equal  thickness 


I  If 


,  S^,  Sy,  Y,  p,  and  the 
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and  the  number  of  parameters  is  reduced  by  two. 

The  values  given  by  equations  (53)  and  (54)  are  not  greatly  in  error  when 
the  facings  are  not  of  equal  thickness.  They  are  only  ten  percent  too  great 
when  one  facing  is  twice  the  thickness  of  the  other. 


Isotropic  Facings  of  Equal  Thickness 
and  Orthotropic  Core 


For  isotropic  facings  of  equal  thickness  some  simplification  can  be  made 
because  of  the  following  relations; 


E  =  E  =  E 

X  y 

(Tjjy  =  (Tyx  =  <r 

(55) 

E 

^"xy  “  2  (1  +  a) 

Substituting  these  values  in  expression  (44) 


4-1  = 


4Tr^Y 


(v  +  P)^  +  (v  - 

[(y  +  P)^  +  [(y  -  P)^  + 


(56) 


Substitution  of  the  values  given  by  equations  (55)  into  equations  (47)  and 
taking  <r  to  be  1/4  yields: 

A ,  =  7  (y'^  +  +  P'^)  +  TT  +  P^) 


A 


2  ■ 


+  P^) 


(57) 


A 


+  P^) 


Substitution  of  the  values  given  by  equations  (49),  (53),  (54),  and  (57)  into 
equation  (51)  yields 
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obtain  the  least  value  of 


Maximum  Buckling  Stress  Associated 
with  Very  Thin  Facings 


I£  the  facings  are  very  thin,  their  individual  flexural  stiffnesses  may  be  neg¬ 
lected.  The  critical  shear  stress  will  then  be  limited  by  that  obtained  when 
the  cylinder  wall  buckles  into  a  great  number  of  very  short  waves.  This 
critical  stress  is  associated  with  the  instability  of  the  core  in  shear.  Equa¬ 
tions  (41)  and  (42)  show  that  the  value  of  expression 

2 

2  tt  r 
Jp  =  -j- 
n'^h 

becomes  very  small  when  the  number  of  buckles  n  becomes  very  large  and, 
therefore,  the  value  of  may  be  neglected.  Also  only  the  expressions  in 

the  first  term  of  4^2  that  have  the  expression  Jp^  in  their  denominators  need 

be  considered.  Further,  because  the  facings  are  very  thin  the  expression 
c  2 

(1  -  — )  in  the  second  term  of  ip,  is  substantially  zero,  so  this  term  may  be 
h  —  c  2 

neglected  and  the  expression  (1  +  |^)  in  the  first  term  of  4<2  becomes  4.  Equa¬ 


tion  (59)  may  then  be  written 


^  8XrS. 


I  +  9  (y^  4  p^) 


Minimum  values  of  occur  when  p  =  0  and  y 
values  and  the  expression  for  0  in  (60): 


Substituting  these 


f  -  4Xr^ 

Substitution  of  the  expressions  for  and  Sy  given  by  equation  (45)  and  re¬ 
membering  that  Ej^  =  E  and  that  c  is  substantially  equal  to  h 


This  is  the  maximum  critical  shear  stress  that  can  be  obtained  if  the  individ¬ 
ual  stiffnesses  of  the  facings  are  neglected. 
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Very  Long  Cylinders 


When  the  cylinder  is  very  long,  it  is  convenient  to  eliminate  the  parameter 
J  by  means  of  the  relation: 


(63) 


obtained  from  equations  (41)  and  (42).  As  the  cylinder  is  lengthened,  the 
number  of  buckles  n  becomes  two  (10)  and  thus  according  to  equation  (41) 
the  value  of  p  approaches  zero.  Thus  equation  (44)  becomes: 


+  (—  -  +  1 
Ey  ^"xy 


(64) 


Equations  (47)  and  (49)  for  substitution  in  equation  (51)  become: 


"  2  4.^%  2 
.  =  —  V  +  — V 


^1  \  ^ 


(65) 


M  E  u 

y  2 

A3  =  2  ^  +  2 


A4  =  V 


Minimum  values  of  Kj  from  equation  (52),  using  these  new  expressions  for 
4/j  and  are  found  by  choosing  suitable  values  of  y. 


If  the  facings  are  isotropic,  equation  (44)  becomes: 
3 

*1  if 

and  equation  (51)  becomes: 


(66) 
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= 


_2h _  (1  + 


K(1  -£)r 


8r 


,3  n2hS 


1  +■ 


^  +  9  +  |(i  +  ev^)] +  yV 


(67) 


If  the  core  is  also  isotropic,  equation  (66)  is  not  changed  but  equation  (67) 
becomes : 

J_ 

h^  .  If 


2h _  (1  +  Y^)^ 


'*'2  “  .  cx  Y 


Ml  -  th 


1  +11^(1  +  7)  h^ 


(68) 


Values  of  ^  are  chosen  that  make  Kj  =  iJjj  +  a  minimum. 

The  shear  force  per  unit  length  of  the  edge  of  the  sandwich  cylinder  is  given  - 
by: 

,2  2 

N  =  Th  =  (1  -  |)h  =  (1  -  |)KjE  ^  =  (1  -  J)(4.j  +  4;2)E  —  (69) 

The  expression  for  N  obtained  by  substituting  the  expressions  for  ijij  and  ^>2 

given  by  equations  (66)  and  (68),  omitting  the  last  term  in  the  brackets  in 
equation  (68),  can  be  shown  to  be  in  agreement  with  the  expression  obtained 
by  Gerard  (3)  by  another  method.  The  omitted  term  is  associated  with  the 
flexural  rigidities  of  the  facings  and  may  add  considerably  to  the  critical  load 
if  the  core  has  a  small  modulus  of  rigidity  and  the  facings  are  thick. 


Numerical  Results 


Numerical  values  of  the  buckling  coefficient  for  cylinders  having  isotropic 

facings  of  equal  thickness  were  obtained  by  use  of  equation  (59)  and  equations 
(56)  and  (58)  and  are  plotted  in  figures  3  to  8.  These  curves  are  reasonably 
accurate  for  sandwich  constructions  having  unequal  facings  if  the  thickness 
of  one  facing  is  not  more  than  twice  that  of  the  other.  Values  are  plotted  for 
J  ranging  from  10  to  10,000,  from  0  to  2,  for  0  equal  to  0.4,  1.0,  and  2.5 
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c 

and  —  equal  to  1  and  0.7.  The  particular  values  of  Q_  were  chosen  to  agree 

with  the  properties  of  honeycomb  cores  of  hexagonal  cells  and  isotropic 
cores.  Estimates  of  K£  for  other  values  of  0  can  be  obtained  by  interpola¬ 
tion.  Values  of  Kf  for  ^  equal  to  1  apply  to  sandwich  panels  having  facings 

so  thin  that  their  bending  stiffnesses  may  be  neglected.  An  estimate  of  the 
effect  of  facing  thickness  is  given  by  the  values  of  Kr  for  ^  equal  to  0.7. 

In  calculating  the  minimum  values  of  the  buckling  coefficient  (K^),  it  was 

assumed  that  the  values  of  p  were  continuous.  Strictly  speaking,  the  values 
of  p  are  discrete  because  they  contain  an  integral  number  of  buckles  n  in 
the  circumference  of  the  cylinder  as  shown  by  equation  (41).  These  discrete 
values  depend  upon  the  ratio  of  the  radius  of  the  cylinder  to  its  length. 

The  results  of  the  calculation  showed  that  the  values  of  the  buckling  coefficient 
(K^)  were  not  sensitive  to  the  values  of  p;  thus,  the  error  introduced  by  the 
assumption  that  p  is  continuous  is  small. 

As  the  cylinder  becomes  longer  the  number  of  buckles  becomes  less  and  is 
equal  to  two  for  an  infinite  length  (10).  Thus,  for  large  values  of  J  the  num¬ 
ber  of  buckles  approaches  two.  The  results  of  the  calculations  showed  that 
for  J  equal  to  10,000,  values  of  p_  are  about  0.  1  for  substantially  all  of  the 
calculations  made  (p  was  less  than  0.1  for  9  equal  to  2.5  when  was  large). 
Substitution  of  these  values  in  equation  (63)  shows  that  for  two  waves  the 
ratio  of  the  radius  of  the  cylinder  to  its  thickness  is  about  40  and  the  assump¬ 
tion  of  a  greater  number  of  waves  leads  to  a  larger  value.  Thus,  for  large 
values  of  J  these  curves  apply  only  to  cylinders  with  a  ratio  of  radius  to  thick¬ 
ness  in  the  neighborhood  of  40  or  larger.  If  this  ratio  is  less  than  40,  the 
curves  give  conservative  values  of  the  buckling  coefficient. 

Numerical  values  of  Kf  for  infinitely  long  isotropic  cylinders  having  equal 
facings  were  calculated  from  the  sum  of  equations  (66)  and  (68)  using  equa¬ 
tions  (53)  and  (54)  and  are  plotted  in  figure  9  against  S  for  different  values  of 
r  c 

— .  The  solid  lines  apply  to  a  —  value  of  unity.  The  uppermost  solid  line  that 

cuts  across  the  others  indicates  shear  instability  of  the  core  and  is  obtained 
from  equation  (62).  The  dotted  lines  apply  to  a  Rvalue  of  0.7. 
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Conclusions 


The  shear  stress  in  the  facings  at  which  buckling”  of  a  cylinder  of  sandwich 
construction  in  torsion  will  occur  is  given  by: 


T 


f 


=  K^.E 


h 
f  r 


For  orthotropic  facings  and  cores,  is  given  by  equation  (52)  in  connection 

with  equations  (44)  and  (51).  For  isotropic  facings  and  orthotropic  cores, 
equations  (56)  and  (58)  are  used  in  place  of  (44)  and  (51).  Numerical  values 
have  been  calculated  and  are  plotted  in  figures  3  to  8.  Maximum  values  for 
the  critical  shear  stress'  for  cylinders  having  orthotropic  cores  and  very 
thin  isotropic  facings  are  given  by  equation  (62),  and  are  plotted  as  the  upper 
solid  curve  in  figure  9. 


Values  of  for  very  long  sandwich  cylinders  having  orthotropic  facings  and 
cores  is  given  by  equation  (52),  using  equation  (64)  and  equations  (65)  in  equa¬ 
tion  (51)  for  values  of  ipj  and  4<2-  Values  for  very  long  cylinders  having  iso¬ 


tropic  facings  and  orthotropic  cores  are  given  by  the  sum  of  equations  (66) 
and  (67).  Values  for  isotropic  cylinders  are  given  by  the  sum  of  equations 
(66)  and  (68).  Numerical  values  of  this  sum  are  plotted  in  figure  9.  These 
values  are  limited  by  equation  (62)  as  shown  in  the  figure. 


—The  possibility  of  wrinkling  at  a  lower  stress  (similar  to  that  described  in 
Forest  Products  Liaboratory  Report  No.  1810),  (9_)  should  be  investigated. 
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Notation 


^1- 

^2 

defined  by  equation  (16). 

A 

1 

^b‘ 

i  1  »  • 

"^2*  ^3*  ^4*  "^5 

defined  by  equations  (29)  to  (33). 

Af. 

II  II  ti 

•^2’  ^3’  ^4 

defined  by  equations  (47)  and  (49). 

b 

B 


e  ,  e  ,  etc. 
XX  xy 

E  ,  E  .  E 
x’ 


g 

h 

H 

I 

I, 


length  of  cylinder. 

J_ 

Ea' 

thickness  of  core. 
Zcr 


1 


xy 


a 

components  of  strain. 

Young’s  moduli  of  the  facings. 

E.  (.  -  ^). 


thickness  of  the  facings, 
coefficient  in  equation  (1). 
c  +  fj  + 

defined  by  equations  (28)  and  (34). 
defined  by  equation  (25). 
defined  by  equation  (26). 

k! 

rh' 
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k,  k' 
K,  Kf 

Kj.  K 


n 


q.  q' 


r 


S 


y 


u 

V 

w 


W, 


W. 


W. 


^x*  ^y' 


a 


Y 

5 


parameters  introduced  by  equations  (23). 

buckling  coefficients  (see  equations  (39). 

defined  by  equations  (17)  and  (18). 

E  <r  ^  +  2Xu. 

X  yx  '^xy 

number  of  buckles  in  the  circumferential  direction. 

parameters  introduced  in  equations  (23). 

radius  of  middle  surface  of  the  cylindrical  shell. 

defined  by  equations  (45). 

axial  component  of  displacement. 

circumferential  component  of  displacement. 

radial  component  of  displacement. 

extensional  strain  energy. 

flexural  strain  energy. 

work  done  by  applied  couple. 

components  of  stress. 

n 
r  ■ 

JT 

b' 

defined  by  equation  (27). 
inclination  of  buckles, 
defined  by  equation  (4). 
defined  by  equation  (4). 
coordinate  shown  in  figure  2. 
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oai  a 


modulus  of  rigidity  of  facings 
moduli  of  rigidity  of  core. 

t^xy  ■  h^- 

1 

Poisson's  ratios  of  the  facing 

irr 
nb 

mean  shear  stress, 
shear  stress  in  facings. 
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-Buckling  coefficients  of  cylindern  havii 


constructions  in  torsion. 


SECTION  V 


ANALYSIS  OF  LONG  CYLINDERS  OF  SANDWICH  CONSTRUCTION 
UNDER  UNIFORM  EXTERNAL  LATERAL  PRESSURE- 
By 

MILTON  E.  RAVILLE,  Engineer 
2 

Forest  Products  Laboratory,—  Forest  Service 
U.  S.  Department  of  Agriculture 


Summary 


A  theoretical  solution  is  obtained  for  the  stresses  in  long  sandwich 
cylinders  subjected  to  externally  applied  uniform  lateral  pressure. 

The  analysis  is  extended  to  take  into  account  failure  due  to  elastic 
instability,  and  an  equation  is  derived  for  the  determination  of  criti¬ 
cal  loads  on  long  sandwich  cylinders.  The  application  of  this  analysis 
to  long  curved  panels  (portions  of  a  long  cylinder)  is  discussed.  The 
sandwich  construction  is  assumed  to  consist  of  isotropic  membrane  fac¬ 
ings  and  an  orthotropic  core. 


Introduction 


In  this  report  it  is  assumed  that  the  sandwich  cylinder  is  comprised  of 
thin.  Isotropic  facings  of  a  relatively  stiff  material  separated  by  an 
orthotropic  core  of  a  relatively  weak  material.  The  facings  are  assumed 
to  be  so  thin  that  bending  and  shear  in  the  individual  facings  may  be 
neglected.  It  is  further  assumed  that  the  only  stress  components  present 
in  the  core  are  the  normal  stress  on  surfaces  parallel  to  the  facings  and 
the  transverse  shear  stress.  Since  elasticity  theory  is  used  in  regard 
to  the  core,  the  solution  is  not  limited  with  respect  to  core  thickness. 
Results  based  on  the  assumption  of  membrane  facings  are  somewhat  limited 


—This  report  is  one  of  a  series  prepared  by  the  Forest  Products 

Laboratory  under  U.  S.  Navy,  Bureau  of  Aeronautics  Order  No.  01595- 

2 

—Maintained  at  Madison,  Wis.,  in  cooperation  with  the  University  of 
Wisconsin. 
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in  applicability,  but  they  should  be  sufficiently  accurate  for  the 
majority  of  sandwich  panels.  It  is  felt  that  the  core  assumptions 
represent  reasonably  well  the  cores  presently  in  use,  especially  those 
of  the  honeycomb  type. 

The  method  of  analysis  is  based  primarily  on  the  fact  that  the  core  assump¬ 
tions  permit  the  direct  determination  of  the  core  displacement  functions 
insofar  as  their  dependence  on  the  radial  distance  is  concerned.  The  usual 
requirement  of  continuity  of  displacements  at  the  interfaces  is  specified. 
The  stability  analysis  is  based  on  the  method  used  by  Timoshenko^  in  the 
analysis  of  the  stability  of  homogeneous  isotrop^c  cylinders.  The  stresses 
in  the  sandwich  cylinder  are  determined  for  loads  less  than  the  critical 
loadj  and,  in  discussing  buckling,  only  small  deformations  from  this  uni¬ 
formly  compressed  form  of  equilibrium  are  considered. 


Notation 


r,  9,  z 

a 

b 

t 

r 


E 

c 


^9^  '^r9 


radial,  tangential,  and  longitudinal  coordinates,  respectively 

radius  to  middle  surface  of  outer  facing 

radius  to  middle  surface  of  inner  facing 

thickness  of  facings 

modulus  of  elasticity  of  facings 

Poisson's  ratio  of  facings 

modulus  of  elasticity  of  core  in  direction  normal  to  facings 

modulus  of  rigidity  of  core  in  r-9  plane 

intensity  of  uniform  external  lateral  loading 

intensity  of  uniform  external  lateral  pressure,  equal  to  -q 

normal  stress  in  the  radial  direction  in  outer  facing,  inner 
facing,  and  core  respectively 

small  normal  stresses  in  the  radial  direction  in  outer  facing, 
inner  facing,  and  core  respectively 

small  shear  stress  in  the  plane  of  the  middle  surface  of  the 
outer  and  inner  facing,  respectively 


~S.  Timoshenko.  Theory  of  Elastic  Stability,  p.  hh6. 


VADD  TR  b 0-133 


154 


small  transverse  shear  stress  in  the  core 


direct  stress  resultants  In  the  tangential  direction  in  the 
plane  of  outer  and  inner  facing,  respectively 


% 


small  direct  stress  resultants  in  the  tangential  direction 
in  the  plane  of  outer  and  inner  facing,  respectively 


u,  u  ,  radial  displacements  of  outer  facing,  inner  facing,  and 

core,  respectively 

u,  u  ,  u  small  radial  displacements  of  outer  facing,  inner  facing, 

and  core,  respectively 

V,  V  ,  V  small  tangential  displacements  of  outer  facing,  inner  facing, 

^  and  core,  respectively 


^9’  ‘^9’  ^9c 


tangential  strains  in  outer  and  inner  facing,  respectively 

small  tangential  strains  in  outer  facing,  inner  facing,  and 
core,  respectively 


number  of  waves  in  circumference  of  cylinder 


one -half  the  central  angle  of  curved  panel 


log  natural  logarithm 

A,  B,  a',  b',  c',  h',  a  ,  B  ,  C  ,  H  arbitrary  constants 

nnnnnnnn 


Theoretical  Analysis 


The  first  step  in  the  analysis  is  that  of  determining  the  stresses  in  the 
sandwich  cylinder  for  pressures  up  to  the  critical  pressure.  The  cylinder 
remains  circular  and  is  in  a  state  of  uniform  compression.  Throughout  the 
analysis,  q  represents  a  uniform,  lateral  load  acting  in  the  positive  radial 
direction.  The  pressure,  p,  is  equal  to  -q. 
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Equilibrium  of  Core 

As  previously  stated,  It  Is  assumed  that  the  core  can  transmit  only  normal 
stresses  In  the  radial  direction  and  transverse  shear  stresses.  When  the 
cylinder  Is  In  a  state  of  uniform  compression,  the  only  stress  present  In 
the  core  Is  the  normal  stress  In  the  radial  direction,  Considering 

the  equilibrium  of  the  differential  element  of  the  core  shown  in  figure  1,  the 
summation  of  forces  In  the  radial  direction  results  in  the  following  equation; 

dq 

-  CT__  r  dS  dz  +  (a_„  +  dr)(r  +  dr)  d6  dz  =  0 

i  C  I  C  QJ* 


This  reduces  to  the  following  differential  equation  of  equilibrium: 


'i^rc  ®rc 

+  _££  =  0 
dr  r 


(1) 


The  solution  of  equation  (l)  is 

a  =  E  - 
rc  c  r 


(2) 


in  which  E.A  represents  the  constant  of  integration. 


of  the  core, 
tlon: 


is  related  to  the  radial  stress,  a 


rc' 


du. 


q  =  E  •  'T  • 
rc  c  dr 


The  radial  displacement 
by  the  following  equa- 

(5) 


From  equation  (2)  it  follows  that 
dr  ~  r 

Integration  gives 

u^  =  A  log  r/a  +  B  (U) 

The  use  of  log  r/a  instead  of  log  r  merely  alters  the  arbitrary  constant  of 
integration. 


Equilibrium  of  Facings 

With  reference  to  the  differential  elements  of  the  facings  shown  in  figure  2, 
it  is  seen  that,  when  the  cylinder  is  uniformly  stressed  by  the  action  of  the 
external  load,  q,  the  only  stresses  present  in  the  facing  are  the  radial 
stresses,  and  q^,  which  are  exerted  by  the  core,  and  the  tensile  forces 

per  unit  length  of  facing,  Ng  and  Ng .  These  stresses  are  assumed  to  be 

acting  on  the  middle  surfaces  of  the  facings;  an  assumption  which  is  justi¬ 
fied  since  the  facings  are  assumed  to  be  thin.  An  equilibrium  equation 
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can  be  obtained  for  each  facing  by  summing  forces  in  the  radial  direction 
on  each  element.  The  equilibrium  equation  which  pertains  to  the  outer 
facing  is 

qa  d0  dz  -  tJj,a  dS  dz  -  Ng  d0  dz  =  0 

This  reduces  to 

=  qa  -  a  a 
0  r 

or,  since  a  -  (a  ) 

’  r  rc 

r  =  a 

Ng  =  qa  -  a(a^^)  (5) 

r  =  a 


In  a  similar  manner,  the  equilibrium  equation  of  the  inner  facing  is  ob¬ 
tained  as  . 


=  b((7  ) 

rc 


r  =  b 


(6) 


The  requirement  is  now  made  that  the  displacements  of  the  core  and  facings 
be  equal  at  tte'  -respective  interfaces.  It  is  assumed  that  the  core  extends 
to  the  middle  surfaces  of  the  facings.  Thus 


u  =  (u  ) 

r  =  a 


and 


u 


(7) 

(8) 


From  Hooke's  law. 


Ng  =  Et(Cg) 

Since  ,  and,  in  view  of  equation  (7),  Sq  = 


c'r  =  a 


,  then 


Ng  =  Eti  (Uj.) 

a  r  =  a 

By  using  equation  (U)  in  conjunction  with  the  above  equation 
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(9) 


N  =  Et(B) 

^  a 

Also,  from  Hooke's  law, 

K  = 

» 

Since  =  ^  ,  the  above  equation  in  conjunction  with  equations  (8)  and 

(4)  leads  to: 

”0  =  |)  (10) 

Since,  from  equation  (2),  (a„)  =  £„■§  and  (cr^„)  =  equations 

r  =  a  r  =  b 

(5)  and  (6)  may  be  written  as 

Ng  =  qa  -  E^A  (ll) 


and 

Ng  =  E^A  (12) 

After  equating  the  right-hand  sides  of  equations  (9)  and  (ll)  and  the 
right-hand  sides  of  equations  (lO)  and  (12),  two  -equations  are  obtained 
from  which  the  values  of  the  constants  A  and  B  may  be  determined.  These 
two  equations  are 

-  V  =  (15) 


and 


E^A  =  (A  log  b/a  +  B)  (l4) 

From  equation  (l4) 

B  = 

Substituting  the  above  value  of  B  into  equation  (l?)  and  multiplying 
through  by  g~. 


!£! 

Et 


log  b/e 


(15) 
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qa^  E^a 
Et  eT~ 


Et 


-  log  b/a 


A 


from  which 


A  = 


Ej.a 

+  b/a)  -  log  b/a 


or 


_ 1  _ 

(1  +  b/a)  - 

E^a 


(16) 


Substitution  of  the  above  value  of  A  into  equation  (l5)  gives 


B 


gab 

Et 


1  - 


(l  +  b/a) 


Et  log  b/a 


E^b 


Et  log  b/a 


Eca 


(17) 


When  the  value  of  A  given  by  equation  (l6)  is  substituted  into  equations 
(2),  (11),  and  (12),  respectively,  the  following  expressions  for  the 
stresses  in  the  cylinder  are  obtained: 


^rc  = 


(18) 


Ng  =  qa  (1  -  k) 

N  =  qak 
0 


where  k 


_ 1 _ 

(1  +  b/a)  .  .^.°g 

III  a 


(19) 

(20) 


Also,  by  using  the  values  of  A  and  B  given  by  equations  (16)  and  (l7), 
equation  (t)  for  the  core  displacement,  after  some  simplification,  becomes 


qa 

Et 


F+ 

(1  -  k)  +  g-^k  log  r/a 
■'c 


(21) 
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and 


Since  u,  the  radial  displacement  of  the  outer  facing,  equals  (u„) 

I  ^ 

u  ,  the  radial  displacement  of  the  inner  facing,  equals  (u_)  , 

r  =  h 

u  =  g^(l  -  k)  (22) 

and 

=  I?  (25) 

The  analysis  thus  far  is  not  limited  strictly  to  membrane  facings.  The 
only  restrictions  on  facing  thicknesses  are  those  imposed  by  the  usual 
thin-walled-cylinder  theory  plus  the  additional  requirement  that  the  radial 
displacement  of  the  middle  surface  of  the  facings  be  equal  to  the  radial 
displacement  of  the  corresponding  Interfaces. 

The  next  step  in  the  theoretical  analysis  concerns  the  development  of  the 
stability  criteria.  In  this  analysis,  which  follows,  the  facings  are 
assumed  to  be  membrane -type  cylindrical  shells.  As  mentioned  previously, 
it  is  assumed  that  the  stress  situation  in  the  deformed  cylinder  differs 
only  slightly  from  the  stress  situation  which  exists  just  before  buckling. 

A  bar  is  placed  over  the  appropriate  symbol  to  denote  the  small  stresses, 
strains,  and  displacements  which  occur  when  the  cylinder  goes  from  the 
initially  uniformly  stressed  circular  form  to  the  slightly  deformed  con¬ 
figuration.  These  quantities  are  dependent  upon  0  as  well  as  r. 


Equilibrium  of  the  Core 


Since  the  cylinder  is  now  considered  to  be  in  a  slightly  deformed  state, 
the  core  is  also  slightly  deformed.  Figure  5  shows  that,  in  addition  to 

the  radial  stress,  q^,  which  exists  just  before  buckling,  an  additional 

small  radial  stress,  and  a  small  transverse  shear  stress,  must 

be  considered.  These  stresses  result  when  the  core  takes  on  small  deforma¬ 
tions  from  the  initially  uniformly  compressed,  circular  shape,  and  they 
are  dependent  upon  0  as  well  as  r.  The  differential  element  shown  in 
figure  3  is  considered  to  be  in  equilibrium  under  the  stresses  shown.  Sum¬ 
mation  of  forces  in  the  radial  direction  yields  the  following  equilibrium 
equation; 


^rc) 


r  d0  dz  + 


(qfk 


+  '^rc  + 


dCtlfk) 


dr 


■dr  + 


rc  - 


-57^0  (r 


+  dr)  d0  dz 


T  -  dr  dz  + 
r0c 


(T 


r0c 


+ 


9^r0c 

do 


d0) 


dr  dz  =  0 
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If  terns  containing  the  products  of  more  than  three  differentials  are  neg¬ 
lected,  the  above  equation  reduces  to 

a^»dr  d0  dz  +  ~^°r  dr  d0  dz  +  -TjC^Sdr  d9  dz  =  0 

5r  be 


or,  dividing  through  by  r  dr  dS  dz. 


1 

rc  ^  rc  ^  _  r9c 
3r  r  r  d0 


=  0 


(24) 


Summation  of  forces  in  the  tangential  direction  yields  the  following  equa¬ 
tion  of  equilibrium; 


dr 


T.rd0dz+(T„  + 
rSc  '  r0c 


*'®^dr)(r  +  dr)  d0  dz  +  x^g^dr  d0  dz  =  0 


which  reduces  to 


St 


r0c 

Sr 


2t 


r0c 


=  0 


(25) 


The  core  displacements  are  related  to  the  core  stresses  by  the  following 
equations : 


0  =  E  s 

rc  c  dr 


(26) 


and 


'^r0c  “  ®r0 


1  Su^ 
r  S0~ 


r 


(27) 


where  0^  and  are  the  small  radieU.  and  tangential  core  displacements 

from  the  uniformly  compressed  form  of  equilibrium.  Equations  (24  -  27) 
are  sufficient  to  enable  the  direct  determination  of  the  displacement 
functions  u^  and  v^  insofar  as  their  dependence  on  r  is  concerned.  Thus, 

it  is  assumed  that  and  7^  are  expressible  in  the  following  form: 


u  =  f, (r)  cos  n0 
c  i 

V  =  f„(r)  sin  n0 
c  2'  ' 


(28) 

(29) 
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and,  from  equation  (30), 


nA^  ®n 

®rc  ~ 


(33) 


Referring  to  equation  (26),  the  following  equation  may  be  writtein: 


Su  1  nA  B 

c  /  n  ,  nv  - 

-  =  — 

hr  E  r"^  r 

c 


d  fi(r) 

Since,  from  equation  (28),  ^ -  cos  nS, 


d  f, (r)  1  nA'  b' 

dr  E„  r  r 

c 


Integrating 


1  nA 

=  e;  ^ 


Therefore 


“  E 


1  nA  ^  ^ 

( - log  r  +  C^)  cos  n0 


(34) 


Rewriting  equation  (27), 


dv  V  T  1  8u 

_ c  _£  _  rRc _ c 

Sr  ”  r  ~  G  ^  '  r  hd 

re 


Substituting  for  and  u^  their  values  giveji  by  equations  (32)  and  (34), 

7  1  a'  1  n^A'  nB*  log  r  c' 

c  c  j  a  ,  I  n.n®  nv.- 

- - - - ^  sin  n0  +  —  ( - - — I-  -  +  — )  Bin  n0 

Sr  r  G  r^^  E  r^  r  r 

r0  c 
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Using  equation  (29)  for  v  , 


d  f  (r)  f  (r)  1  n2  a'  nB*  log  r  1  c' 

2  2 _ ^ ^  _n  ^  _ n  _  ^ _ r 

r  G^Er^E  r  Er 

r0  c  c  c 


dr 


The  solution  of  the  above  equation  is 

1  <^2 


f^(r)  =  H V  -  ( 


n^  A„  nB„  C„ 

2G  .  ■  r”  ■  E  E“ 

rS  c  c  c 


Therefore 


V  = 

c 


1  n- 
're  c 


2  a’  nB' 


/  _ \  n  n  _  .  n 

■^2G-'2E'r  “E  ‘  E  ■*■ 


sin  ne  (35) 


The  functional  dependence  on  r  of  the  displacements  u^  and  as  expressed 
in  equations  (5*+)  and  (35)  having  been  determined,  it  is  convenient  at  this 


I  r 


point  to  redefine  the  constants  A^,  B^,  C^,  and  H^.  Equation  (34)  may  be 


written  as  follows 
u 


c  =  ®n  7“  ■‘-n^  r/a)  cos  ne 

Then  equation  (35)  must  be  replaced  by 


1  E„  n  a2 
c 


-  nA  a  +  (—  ;; —  -  — )B  —  -nCa(l  +  log  r/a)  +  H  r 
n  '2n  G„2  nr  n'  ^ 

re 


(56) 


sin  n9 
(37) 


The  displacement  equations  are  thus  expressed  in  terms  of  the  simpler  non- 

dimensional  arbitrary  constants  A  ,  B  ,  C  ,  and  H  .  Equations  (32)  and 

n'  n'  n  n  ^ 

(33)  for  the  core  stresses  are  expressed  in  terms  of  the  new  constants  as 
follows : 


'^rec 


E 

c 

-  —  Bn~2 


(38) 


and 
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q  (l  +  €_)  a  de  dz  -  (qk  +  5  )(l  +  e_)  a  d0  dz 


qa  (l  -  k)  +  Ng  (1  +  -  —  -  i  d9  dz  =  0 

a  bo  a  bff 


If  terms  containing  products  of  small  quantities,  l.e.,  products  of  barred 
quantities,  are  neglected  and  the  above  equation  is  divided  by  d0  dz,  the 
result  may  be  written  as 

/  V-  / 

N  =  -  CT  a  +  qa  (1  -  k)e  -  qa  (l  -  k)( - - - o) 

dr  0  a  be  a  b6‘^ 


Since  ~  a  ^  a  above  equation  further  reduces  to 

a  -  /,  ,  s/Q  .  1 


N  =  -  5  a  +  qa  (l  -  k)(—  + - p) 


The  second  equilibrium  equation  of  the  outer  facing  is  obtained  by  summing 
forces  in  the  direction  tangent  to  the  element.  Thus, 

dz  -  ^rg(l  +  Sg)  a  de  dz  =  0 

Again  by  neglecting  the  term  containing  the  product  of  the  barred  quantities 
and  dividing  through  by  iO  dz,  the  above  equation  is  reduced  to 


'Se~  ""  '^rO^ 

The  two  equilibrium  equations  which  apply  to  the  inner  facing  are  obtained 
in  a  manner  exactly  analogous  to  that  used  in  obtaining  equations  (4o)  and 
(4l).  These  equations  are: 

n'  =  a  b  +  qak  (-  +  -  (42) 

^  ^  h  h  be^ 


- 


In  addition  to  the  four  equilibrium  equations  (40  -  h3) ,  the  following  two 
equations,  based  on  Hooke's  law  and  the  strain-displacement  relations,  may 
be  written: 
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(59) 


5  =E  (-B^+C^)  cos  nS 
rc  c  '  nr' 

The  validity  of  the  above  equations  is  easily  verified  by  the  substitution 
of  equations  (38)  and  (39)  into  the  equilibrium  equations,  equations  (24) 
and  (25),  and  further  by  the  substitution  of  the  expressions  given  by  equa¬ 
tions  (36  -  39)  into  the  stress -displacement  relations,  equations  (26)  and 
(27).  It  is  of  interest  to  note  that  the  manner  in  which  the  core  stresses 
and  displacements  vary  with  respect  to  r  is  not  so  simple  as  is  sometimes 
assumed  in  problems  of  this  nature. 


Equilibrium  of  Facings 


Immediately  before  buckling,  the  stress  situation  in  the  facings  is  as 


illustrated  in  figure  2.  Since  a  =  (a  )  and  c  =  (c  )  ,  the 

r  =  b  r  =  a 

stresses  shown  in  figure  2  can  be  determined  from  equations  (18),  (I9),  and 
(20).  As  stated  previously,  the  stress  situation  in  the  facings  when  the 
cylinder  is  in  a  slightly  deformed  state  is  assumed  to  differ  only  slightly 
from  the  situation  which  existed  just  before  buckling.  In  figure  4,  which 
illustrates  the  differential  elements  of  the  deformed  facings,  the  symbols 

No,  No,  T  a,  To,  a  ,  and  a  are  used  to  indicate  this  small  stress  varia- 
9’  6’  r6’  tQ’  r’  r 

tion  which  has  taken  place.  The  assumption  of  membrane  facings  eliminates 
the  necessity  of  considering  bending  and  transverse  shear  in  the  individual 
facings.  In  formulating  the  applicable  equilibrium  equations  of  the  fac¬ 
ings,  account  must  be  taken  of  the  rotation  and  stretching  of  the  facing 
elements  during  buckling.  This  effect  was  found  to  be  negligible  in  regard 
to  the  core.  Due  to  the  difference  in  rotation  between  the  longitudinal 
elements  M  and  N  of  the  outer  facing  and^the  difference  in  rotation  be¬ 
tween  the  longitudinal  elements  M  and  N  of  the  inner  facing,  the  initial 


central  angle,  d0,  becomes 


8^u 


a  a  bo 


bv 


1  S^u' 


)  de  and  (l  +  -  P 

b  80  b  he 
a  result  of  the  tan 


■)  d0 


for  the  outer  and  inner  facings,  respectively.—  As 
gential  strain  in  the  outer  and  inner  facings,  the ^ areas  of  their  differ¬ 
ential  elements  become  (l  +  eg)  ad0  dz  and  (l  +  eg)  b  d0  dz,  respectively. 


u,  V, 


and  e.  represent  small  displacements  and  strains  which 
0 


result  when  the  cylinder  goes  from  the  undeformed  to  the  deformed  state. 


With  reference  to  figure  4  it  is  seen  that  two  equilibrium  equations  can 
be  written  for  each  facing.  Considering  the  outer  facing  first,  the  sum¬ 
mation  of  forces  in  the  direction  normal  to  the  differential  element  of 
the  outer  facing  yields 
4  ^ 

~S.  Timoshenko.  Theory  of  Elastic  Stability,  p.  431. 
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(hh) 


Ne  -  (is)  ■  -JS 

1-H  1-fi  a  a  bd 


vs'  Et  s  Et  /u  ,  1  9v  s 

Ng  =  - - 2  (Cg)  =  - - 5  (-  + - ) 


(^^5) 


1-4“  “  1  -  4  b  b  Se 

The  following  relations  enable  the  right-hand  sides  of  equations  (40  -  4-5) 
to  be  expressed  in  terms  of  the  constants  A^,  C^,  and  which  appear 

in  the  equations  for  tha  core  displacements,  equations  (56)  and  (57): 

3r  =  (3rc)  ^r  = 

r  =  a  .  r  =  b 

or,  on  the  basis  of  equation  (59), 

5^  =  E^  (-  cos  n6  and  5*  =  E  (.  B  \  +  C  #)  cos  n0 


r  c  n  n 

^r0  "  ^^r0c^ 

r  =  a 

or,  on  the  basis  of  equation  (58), 


(46) 


and 


f  =  (t  ) 
r0  r0c 


r  =  b 


r  .  - - B  sin  n0 

rc7  n 

n 


u  =  (u^) 


E  2 

and  =  .  _£  ^  sin  n0  (4?) 

n  b 


and  u'  =  (u  ) 

'  r*  ' 


r  =  b 


r  =  a 

or,  on  the  basis  of  equation  (56), 

u  =  (A^a  +  B^a)  cos  n0  and  u  =  (A^a  +  B^a  —  +  C^a  log  b/a)  cos  n0 

(48) 


V  =  (v  ) 

c 

r  =  a 

or,  on  the  basis  of  equation  (57), 

V  = 


and  V  =  (v  ) 


r  =  b 


1  E„  n 

■  r - p)  -  n  C  a  +  H  a 

n  £in  d  n  n  n 


sin  nO 
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and 


.1 


-  n  A^a  +  (—  - 2^  ~  "  ■‘' 


re 


sin  nO 

(49) 


The  relationships  given  by  equations  (46  -  49)  express  the  assumed  con¬ 
tinuity  conditions  on  the  stresses  and  displacements  at  the  Interfaces . 

The  further  assumption  that  the  core  extends  to  the  middle  surfaces  of  the 
facings  is  implied. 


On  the  basts  of  the  above-listed  continuity  conditions,  equations  (40  -  45) 
can  be  transformed,  respectively,  to  equations  (50  -  55) >  shown  below. 


^0  = 


aNg 

Se“ 


E  a  (B  -  C  )  -  (n^  -  l)  qa  (l  -  k)(A  +  B  ) 
c  n  n  n  n 


E  a 

- ^B„  sin  nS 

n  ^ 


c  os  n0  ( 50 ) 

(51) 


Ne  = 


-  E  a  (B  r  -  C  ) 
c  n  b  n' 

(n^  -  1)  qak  f  (\  +  E  ^n  Z^/a) 


cos  nd 


E„  b  a 

-  =  ~ —  B„  —  sin  n0 

50  n  b^ 


(52) 

(55) 


^0  = 


Ng  = 


Et 


1  - 

Et 

1  -  d" 


/  2  .  /  En  n2  V  2 

-  (n  -  1)  A  +  ( — = - -  +1)B  -  n  C  +nH 

n  2G  2  ”  '■ 

r0 


-  (n^  -  1)  A  -  +  '(-^C - 

b  2Gj,0  2 


2  2 


n^  +  (n^  -  1)  log  b/a 


C  a/b  +  n  H  i 
n  '  n  J 


cos  n0 


The  integration  of  equations  (5I)  and  (55)  yields 


cos  n0 

"(54) 


(55) 


Ec  a 


0  -  n2  “n 


B  cos  n0 


(56) 


and 


E  b  a 
— —  B„  —  cos  n@ 


(57) 
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The  integration  constants  are  zero,  since  Ng  and  Ng  are  dependent  entirely 
on  9,  as  shown  by  equations  (50)  and  (52). 


Four  independent  equations  containing  the  constant  and  as 

well  as  the  load  q,  can  be  obtained  from  equations  (50),  (52),  (54),  (55), 
(56),  and  (57)-  If  the  right-hand  sides  of  equations  (50)  and  (56)  are 
equated,  the  result  is 

O  ^ 

E  a  (B  -  C  )  -  (n^  -  1)  qa  (1  -  k)(A  +  B  )  =  B 

c  n  n'  /IV  ' '  n  n'  2  n 

n 

The  above  equation  may  be  written  as  follows; 


(n^  -  l)(q/E^)(l  -  k)  A^ 


T?  -  1. 


-  (n^  -  l)(q/E  )(1  -  k)  -  (^^— i-i)  B  -  C  =  0 


If  the  right-hand  sides  of  equations  (54)  and  (56)  are  equated,  the  result 
is 


1  -  u 


/  2  \  /  ^  v  2 

{n  -  1)  A  +  i— - +1)B  .  n  C  +  nH 

n  '  n  n  n 

2Gr0  2 


E  a 


which  reduces  to 


(n^  -  1)  A„  + 


h  2G  „  ■  2  ^ 

t9 


Eo  a  (1  - 


n^  Et 


-  n  C  +  n  H  =  0  (5c 

n  n 

If  the  right-hand  side  of  equations  (52)  and  (57)  are  equated,  the  result 
is 


■  ^  (®n  t  •  ^n^  -  -  1)  qak  I  (A^  +  B^  f  log  b/a) 


b  n  n  b  n 


E  a  a 
n  b 


which  reduces  to 
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n2  -  1. 


-  l)(q/E^)k  +  |  (n2  -  l)  (q/E^)k  ^  +  (-  •■  )  | 


+  (rfi  -  l)(q/E^)k  log  b/a  -  ^  = 


Finally,  if  the  right-hand  sides  of  equations  (55)  and  (57)  are  equated, 
the  result  is 


1  -  n' 


r  /  ?  V  a  .  Ec  n2  a^ 

j  J  -  (n^  .  1)  A  -  +  (— S - +  1)  B 

-  (_  b  2G^g  2  b^ 

-  p  +  (n2  -  1)  log  b/aj  f  +  n  - 


which  reduces  to 


-  (n^  -  1)  -  A„  +  (^  .  _  +  1)  r  + 


a  Ec  a  (l  -  4  )  a 


n2  Et 


-  |^n2  +  (n2  -  l)  log  b/aj  ^  +  n  =  0  (6] 

Since  appears  only  in  equations  (59)  and  (61),  it  is  eliminated  immedi¬ 

ately  by  subtracting  equation  (59)  from  equation  (61).  The  result  is 
_  a  f  E^  n^  -  a^  E„  a  (l  - 


(n^  -  1)(-  -  1)  A  +  (-^ - +  1)(-  -  1)  + 

1_  n  r\  ^ 


(b  ®n  -  (b*  -  1)  +  b 


Et 


.]  c„  -  0 


Equations  (5^),  (60),  and  (62)  comprise  a  system  of  three  equations  con¬ 
taining  the  constants  A^,  and  which  occur  in  the  expressions  for 

the  displacement  functions  u^  and  .  A  buckled  form  of  equilibrium  is 

possible  only  if  equations  (58),  (60),  and  (62)  yield  non-zero  solutions 
for  these  constants;  this  requires  that  the  determinant  of  the  coefficients 
of  A^,  and  be  equal  to  zero.  The  equation  used  for  the  determina¬ 

tion  of  the  critical  load  is  obtained  from  this  determinant.  Specifically, 
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(n‘^  -  1)  q/E^  (1  -  k) 


(n2  -  1)  (q/E  )k 
c 


(n2  -  l)(f  -  1) 


(n'^  -  1)  q/E^  (l  -  k) 
-  1 


(n2  -  1)  (q/E  )k  f 


-  1 


V  ^2  „2 

Ej,  n  a 


(57^  -  -+  1)(—  -  1) 


(n2  •-  1)  (q/E  )k  log  b/ 
c 


b 
a 

2  /a 


E^  a  (1  -  ^2)  a 


2 

n  Et 


(-+1)  +(n2-i)^  log  b/a 


The  following  quadratic  equation  is  obtained  from  the  expansion  of  this 
determinant: 


(q2/E^2)  k  (1  .  k) 


r  2  (1  - 


2> 


(b/a)‘ 


E  n2  1  -  b2/a2 

2G^g  2  b2/a2 


+  Ee  a  (1  -  p2)  ^1  +  b/a^ 
n^  Et  b/a 


log  b/a  1  +  q/E  , 

fd  -  b/a)2r 

J  J 

L  (b/a)-  L 

b/a 


-  k  (1  +  b/a)]  +  [1  -  k  (1  -  b/a)] 


n^  -  1 


-  I^b/a  +  k  (1  -  b/a)] 


Ec  .  .1  -  b2/a^ 


(-^  --+!)( 
2Gre  2 


k2/  2 

b  /a 


E-  a  (1  -  p2)  1  +  b/aj]  1  (l  -  b/a) 


2 

n  Et 


■H  r  "2 

b/a  In  b/a 


=  0 


(65) 


Analysis  of  Results 


The  results  of  this  report  are  contained  in  equations  (18),  (l9)j  and  (20), 
the  equations  for  detennlning  the  stresses  in  a  long  sandwich  cylinder  under 
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unlforro  external  loading;  and  In  equation  (63),  the  equation  for  deter¬ 
mining  the  critical  load  on  a  long  sandwich  cylinder.  These  equations 
are  repeated  below; 

a 

a  =  q  -  k 
rc  ^ 


Ng  =  qa  (1  -  k) 


Ng  =  qak 

(q^/E  2)k  (1  -  k) 


(18) 

(19) 

(20) 


^  2  *  b/a)^ 

(b/a)2 


*  P  P  P 

E-  n  1  -  b  /a^ 

(-  +  -)( - ^-^) 

2G  2  b  /a 

r9 


E  a  (1  -  n  )  1  +  b/s 


n^  Et 


b/e 


/• 

1 

1  log  b/a 

\  +  9/Ec  •  ■ 

J 

1  1 

(1  -  b/a)^ 


(b/a)‘ 


^b/e 


-  K  (1  .  Vail  +  ^  [1  -  k  (1  -  b/a)] 


n^  -  1 


jb/a  +  k  (1  -  b/a) 


E-  n^  1  .  b2/a2 

(—  -  -  +  1)( - ;rV) 

2G  2 
re 


v2/  2 

b  /a 


E^  a  (1  -  n  )  ^1  +  b/a^ 
n^  Et  b/a 

where,  in  each  of  the  above  equations, 

1 


1  (1  -  b/a)‘ 

n^  b/a 


=  0 


(65) 


k  = 


(1  t  b/a)  - 

c 


and  log  b/a  signifies  the  natural  logarithm  of  b/a. 

The  application  of  equations  (18),  (l9)>  and  (20)  to  specific  examples  is 

E  a 

very  simple.  If  the  ratios,  _£ _  and  b/a,  are  known,  the  value  of  k  can 

E  t 

be  determined.  Substitution  of  this  value  of  k  into  equation  (18)  yields 
an  expression  for  the  core  transverse  normal  stress  in  terms  of  the  uni¬ 
form  lateral  load,  q,  and  the  ratio  of  the  radial  distance  to  the  middle 
surface  of  the  outer  facing,  a,  to  the  variable  radial  distance,  r.  Tlie 


WADD  IR  60-133 


172 


maximum  nomal  stress  In  the  core  always  occurs  at  the  inside  surface  of 
the  core,  that  is,  at  r  =  h.  Substitution  of  the  k  value  into  equations 
(19)  and  (20)  yields  expressions  for  the  forces  per  unit  length  of  the 

I 

facings,  Ng  and  Ng,-  in  terms  of  the  load,  q,  and  the  radial  distance  to 
the  middle  surface  of  the  outer  facing,  a.  Examination  of  the  expression 

t 

for  k  shows  that,  in  general,  Ng  and  Ng  are  approximately  equal  since  k 

is  approximately  equal  to  one-half  in  the  practical  range  of  dimensions 
and  elastic  constants.  If  the  value  of  the  modulus  of  elasticity  of  the 


core  approaches  Infinity,  k  becomes  equal  to 

'  a  ' 

and  Ng  =  q  —  -■  Thus,  the  statement  that  Ng  and  Ng  are  equal  if  E^, 

is  infinite  and  b/a  1  is  approximately  correct  and  becomes  less  accurate 


as  the  thickness  of  the  cylinder  Increases.  If  E^— >0,  k  — >0;  and  the 

external  load  is  supported  entirely  by  the  outer  facing.  The  solution 
given  by  equations  (18),  (19)^  and  (20)  may  be  compared  to  Reissner's 
solution  of  the  problem  of  the  closed  circular  sandwich  ring  acted  upon 

by  a  uniform  radial  load .5  If  desired,  equations  similar  to  equations 
(18),  (19),  and  (20)  can  be  easily  obtained  for  the  case  of  Internal 
pressure,  or  for  a  combination  of  internal-  and  external  pressure,  by  the 
methods  used  in  this  report. 


Equation  (63),  for  the  determination  of  the  critical  load  on  long  sand¬ 
wich  cylinders,  is  of  practical  interest  only  for  n  =  2.  A  value  of  n  =  1 
represents,  physically,  a  rigid-body  translation  of  the  cylinder,  and 
values  of  n  >  2  result  in  larger  values  of  the  critical  load  than  those 
obtained  if  n  =  2.  If,  in  equation  (65),  the  Integer  2  is  substituted  for 
n  and  the  symbol  p  is  used  to  represent  a  uniform  lateral  pressure,  that 
is,  -q,  the  result  is 


(p2/E  2)  k  (1  -  k) 


r  (i_^^ 

\  (b/a)2 


,  Ec  .A  - 

(— ^  +  2)( - 2^) 

2Gre 


E  a  (1  -  p  )  1  +  b/a 

^ -  ( - ) 


h  Et 


b/a 


log  b/a 


p/Ec-l 


(1  -  b/a)' 


(b/a)' 


[b/a 


-  .  a  .  Va)l .  I  -  MI  -  bW]  -i 


(-^  - 1)( 


1  -  b2/a2^  Eg  a  (1  -  u^)  ^1  +  b/a^ 


V.2/  2 

b  /a 


-)  + 


h  Et 


b/c 


-  k  (1  -  b/a) 
a 

1  (1  -  b/a)2 


=  0 


b/c 


(6M 
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In  the  practical  range  of  the  values  of  the  dimensions  and  the  physical 
constants  of  the  sandwich  cylinder,  equation  (6k)  yields  two  widely  sep¬ 
arated  positive  roots  for  p,  the  lower  of  which  represents  the  critical 
pressure  •  Consequently,  p^_  may  be  obtained  with  very  good  accuracy 

if  the  term  containing  p^  in  equation  (614-)  is  neglected. 


It  is  of  interest  to  examine  the  res'ults  given  by  equation  (64)  in  cer¬ 
tain  limiting  cases.  First,  it  is  noted  that  the  critical  pressure  be¬ 
comes  equal  to  zero  if  either  or  is  set  equal  to  zero.  This  is  to 

be  expected  in  view  of  the  assumption  of  membrane  facings.  li;  is  of 
greater  interest  to  examine  equation  (64)  with  E^  and  taken  to  be 

infinite.  In  this  case,  equation  (64)  becomes 


(1 


i2) 


■  cr 


Et 


3  (1  -  b/a)g 

2,  2 

1  +  b  /a 


(65) 


The  value  of  p  obtained  from  equation  (65)  should  be  comparable  to  the 
cr 

value  obtained  from  the  formula  for  the  critical  pressure  on  a  long  homo¬ 
geneous  cylinder  if  the  moment  of  inertia  used  is  equal  to  the  moment  of 
inertia  of  the  spaced  facings  of  the  sandwich  cylinder.  The, formula  for 
the  critical  pressure  on  a  long  homogeneous  thin  cylinder  is— 


3  E  I 


cr 


5  (1  -  p2) 


(66) 


-o 


If  the  mean  radius,  r  ,  is  taken  equal  to  the  mean  radius  of  the  sandwich 


cylinder,  and  I  is  made  equal  to  the  moment  of  inertia  of  the  spaced 

2  t  ( a  b )  ^ 

sandwich  cylinder  facings,  — -- g- — equation  (66)  may  be  written  as 
follows ; 


p  a  (1  -  p^)  5  (1  -  b/a)^ 

cr 


(1  +  b/a 


(67) 


Comparison  of  equations  (65)  and  (67)  shows  that  they  would  be  equal  if 


(1  +  — ^)  were  equal  to 


(1  +  b^/a^)  «  for  thin 


O  y 

cylinders,  i.e.,  for—^  1.  As  the  b/a  ratio  decreases  from  a  value  of  1, 

7 - - - 

O  _  , 

“S .  Timoshenko.  Theory  of  elasticity,  p.  216. 
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equation  (67)  becomes  less  accurate  since  it  is  derived  on  the  basis  of 
its  being  a  thin  cylinder  with  the  pressure  applied  on  the  middle  surface 
and  thus  fails  to  represent  a  thicker  cylinder  with  the  pressure  applied 
on  the  outer  surface . 

Equation  (65)  may  also  be  applied  to  problems  concerning  the,  stability  of 
long  sandwich  panels  in  the  form  of  a  portion  of  a  cylinder  hinged  along 
the  edges  0=0  and  9  =  2  a  (fig.  5)*  Ifj  in  equation  (63),  n/a  is  sub¬ 
stituted  for  n  and  -p  is  substituted  for  q,  the  smaller  value  of  p  given 
by  equation  (63)  represents  the  critical  pressure  on  a  panel  whose  dimen¬ 
sions  and  properties  are  known.  This  solution  presupposes  the  unsyrametrical 
type  of  buckling  indicated  in  figure  and  therefore  it  may  not  yield  the 
lowest  critical  pressure  for  panels  of  small  curvature,  that  is,  relatively 
flat  panels .  A  relatively  flat  sandwich  panel  may  buckle  symmetrically  with 
no  inflection  points  between  the  supnorts,  as  in  the  case  of  homogeneous 
7 

panels  .— 


Conclusions 


It  is  felt  that  the  solutions  presented  in  this  report  are  accurate  for 
sandwich  cylinders  with  materials  and  dimensions  such  as  to  render  the 
basic  assumptions  applicable.  The  core  assumptions  are  probably  suffi¬ 
ciently  representative  of  all  practical  sandwich  construction.  The  range 
of  applicability  of  the  solutions  could  be  increased  sanewhat  if  the 
effect  of  the  thicknesses  of  the  facings  on  the  overall  stiffness  of  the 
cylinder  were  taken  into  account  in  the  stability  analysis.  This  would 
entail  the  use  of  shell  theory  rather  than  membrane  theory  in  regard  to 
the  facings.  Another  extension  of  the  problem  which  is  of  greater  prac¬ 
tical  importance  would  be  to  obtain  the  corresponding  solutions  for 
cylinders  of  finite  length.  A  supplementary  report  containing  the  solu¬ 
tions  for  cylinders  of  finite  length  is  planned  for  the  near  future. 
Numerical  computations  and  curves  based  on  the  solutions  contained  in 
this  report  have  been  omitted  since  they  will  appear  as  limiting  cases  in 
the  supplementary  report. 


~S .  Timoshenko. 


Theory  of  Elastic  Stability,  p.  230. 
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gure  2. --Differential  elements  of  Inner  and  outer  facings 
of  uniformly  stressed  cylinder. 


Figure  3. — Differential  element  of  core  of  slightly  deformed 
cylinder,  neglecting  changes  in  geometry  of  core. 


Figure  4.— Differential  elements  of  outer  and  inner  facings 
of  slightly  deformed  cylinder. 


-Cross  section  of  long  panel. 


SECTION  VI 
Supplement  to 


AN.ALYSIS  OF  LONG  CYLINDERS  OF  SANDWICH  CONSTRUCTION 
UNDER  UNIFORM  EXTERNAL  LATERAL  PRESSUREj^ 
Facings  of  Moderate  and  Unequal  Thicknesses 

By 

MILTON  E.  RAVILLE,  Engineer 

p 

Forest  Products  Laboratory,—  Forest  Service 
U.  S,  Department  of  Agriculture 


Summary 


A  previous  analysis  of  the  problem  of  sandwich  cylinders  subjected  to 
uniform  external  lateral  loading  Is  extended  In  order  that  results  may 
be  applied  to  sandwich  cylinders  having  relatively  thick  facings  of 
unequal  thickness.  In  the  development  of  the  stability  criteria,  the 
effect  of  the  stiffness  of  the  Individual  facings  on  the  stability  of 
the  composite  cylinder  is  taken  into  account.  "Solutions  are  obtained 
from  which  the  stresses  and  displacements  in  a  stable  sandwich  cylinder 
may  be  determined,  and  an  expression  for  the  determination  of  the  load 
at  which  a  sandwich  cylinder  becomes  elastically  unstable  is  derived. 
The  sandwich  cylinder  Is  assumed  to  consist  of  isotropic  shell  facings 
and  an  orthotropic  core. 


Introduction 


This  report  Is  a  supplement  to  a  previous  report  that  contains  a 
theoretical  analysis  of  sandwich  cylinders  acted  upon  by  uniform  ex¬ 
ternal  lateral  pressure.  Formulas  were  derived  for  the  stresses  in 
the  cylinder,  and  an  expression  was  obtained  from  which  critical 


•^This  report  is  one  of  a  series  prepared  by  the  Forest  Products 

Laboratory  under  U.  S.  Navy.  Bureau  of  Aeronautics  Order  No.  01595 
2  and  U.  S.  Air  Force  Order  No.  A.F.  18(600;-102. 

“Maintained  at  Madison,  Wis.,  in  cooperation  with  the  University  of 
Wisconsin. 
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pressures  may  be  determined.  In  that  report  It  was  assumed  that  the 
facings  of  the  cylinder  were  thin  enough  to  render  membrane  theory 
applicable  and  that  the  facings  were  of  equal  thickness. 

The  purpose  of  this  supplementary  report  is  to  present  solutions  for 
the  stresses  and  for  the  critical  pressures  that  apply  to  sandwich 
cylinders  having  moderately  thick  facings  of  unequal  thickness.  This 
requires  that  the  bending  moment  and  the  transverse  shear  in  the 
individual  facings  be  considered  in  the  development  of  the  stability 
criteria.  This  extension  of  the  previous  work  is  felt  to  be  of  import¬ 
ance  in  view  of  the  fact  that  an  analysis  based  on  membrane  facings  may 
prove  inadequate  as  a  design  criterion  in  cases  where  relatively  thick 
facings  are  used.  It  is  assumed  throughout  that  buckling  takes  place 
at  stresses  below  the  elastic  limit  of  the  sandwich  materials. 

The  method  of  analysis  used  here  follows  closely  the  method  used  in 
the  original  report.  The  same  core  assumption,  namely,  that  only 
transverse  shear  stress  and  normal  stress  on  planes  parallel  to  the 
facings  are  present,  is  used.  It  is  felt  that  this  assumption  applies 
very  well  to  all  practical  sandwich  constructions  because  of  the  rela¬ 
tively  low  load-carrying  capacity  in  the  tangential  direction,  of  the 
core  materials  as  compared  to  that  of  the  facing  materials.  The  facings 
are  assumed  to  be  homogeneous  and  isotropic,  and,  as  Indicated  pre¬ 
viously,  they  are  analyzed  on  the  basis  of  shell  theory  rather  than 
membrane  theory . 


Notation 

r,  0  polar  coordinates 

a  radius  to  middle  surface  of  outer  facing 

b  radius  to  middle  surface  of  inner  facing 

t^  thickness  of  outer  facing 

tj  thickness  of  inner  facing 

t  thickness  of  either  facing  when  t^  =  tj^ 

E  modulus  of  elasticity  of  facings 

H  Poisson's  ratio  of  facings 

Ej,  modulus  of  elasticity  of  core  in  direction  normal  to  facings 

G  -  modulus  of  rigidity  of  core  in  r  -  0  plane 

rt7 
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q  intensity  of  uniform  external  lateral  loading 

CTj,,  Op,  normal  stress  in  the  radial  direction  acting  on  the  outer 

facing,  on  the  Inner  facing,  and  in  the  core,  respectively 

_  _i  _ 

Op,  Op,  Opg  small  normal -stress  in  the  radial  direction  acting  on  the 
outer  facing,  on  thd  Inner  facing,  and  In  the  cpre, 
respectively 

_  _i 

^r0  small  shear  stresses  acting  on  outer  and  inner  facings, 
respectively 


TpQj,  small  transverse  shear  stress  In  the  core 

I 

Ng,  Ng  direct  stress  resultants  in  the  tangential  direction  In 

the  plane  of  outer  and  inner  facing,  respectively 


^0  small  direct  stress  resultants  in  the  tangential  direction 

in  the  plane  of  outer  and  inner  facings,  respectively 


Qqj  Qq 


small  bending  moments  per  unit  length  of  outer  and  Inner 
facing,  respectively 

small  resultant  shear  forces  per  unit  length  of  outer  and 
inner  facing,  respectively 


u,  u 


radial  displacements  of  outer  facing,  inner  facing,  and 
core  respectively 


u,  u',  u;. 


small  radial  displacements  of  outer  facing,  inner  facing, 
and  core,  respectively 


small  tangential  displacements  of  outer  facing,  inner 
facing,  and  core,  respectively 

t 

^0  unit  tangential  strains  in  outer  and  inner  facing, 

respectively 


small  unit  tangential  strains  in  outer  facing,  inner  fac¬ 
ing,  and  cpre,  respectively 


small  changes  in  curvature  of  outer  and  inner  facing, 
respectively 


n 


number  of  waves  in  circumference  of  cylinder 


a 

P 


one-half  the  central  angle  of  curved  panel 

E^a  (1  - 
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7 


qa  (1  -  |i  ) 
Et„ 


2 


12a 


12b 


Eto  (1  -  b2/a2) 
2G^g  b  (1  -  ^2) 


A,  B,  A„,  B„,  C„, 


arbitrary  constants 


Stress  Analysis 


The  sandwich  cylinder  is  assumed  to  be  long  enough  that  the  effect  of 
the  constraints  at  the  ends  is  negligible.  Under  a  condition  of  uni¬ 
form  external  loading,  each  cross  section  remains  circular.  The  dimen¬ 
sions  of  a  cross  section  of  the  cylinder  and  the  positive  directions  of 
the  polar  coordinates,  r  and  6,  are  indicated  in  figure  1.  In  order  to 
avoid  any  confusion  in  regard  to  signs,  the  Intensity  of  the  external 
load  is  assumed  to  act  in  the  positive  r  direction;  obviously,  a  nega¬ 
tive  value  of  £  signifies  compressive  loading.  In  the  analysis  which 
follows,  the  assumption  is  made  that  the  core  extends  to  the  middle 
surfaces  of  the  facings  and  that  the  load  is  applied  to  the  middle 
surface  of  the  outer  facing.  This  assumption  amounts  to  neglecting  the 
half -thicknesses  of  the  facings  as  compared  to  their  radii. 


Equilibrium  of  Core 


As  previously  stated,  the  assumption  is  made  that,  in  general,  the  core 
transmits  only  normal  stresses  in  the  direction  perpendicxalar  to  the 
facings  and  transverse  shear  stresses.  In  the  case  of  uniform  external 
loading,  it  is  noted  that  the  transverse  shear  stress  is  zero  from  con¬ 
siderations  of  symmetry,  and  the  only  stress  present  in  the  core  is  the 
normal  stress  in  the  radial  direction,  •  Considering  the  equilibrium 
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of  the  differential  element  of  the  core  shown  in  figure  2,  the  sum¬ 
mation  of  forces  in  the  radial  direction  results  in  the  following 
equation: 


=  0 


The  differential  element  is  considered  to  be  of  unit  length  in  the 
longitudinal  direction.  The  above  equilibrium  equation  reduces  to 


d0 


rc 


dr 


0 


(1) 


The  solution  of  equation  (l)  may  be  written  as 


a 


rc 


A 

r 


(2) 


where,  for  convenience,  the  constant  of  Integration  is  represented  by 
Eg  A.  Ej,  represents  the  modulus  of  elasticity  of  the  core  in  the  radial 

direction,  and  A  is  an  arbitrary  constant.  The  radial  displacement  of 
the  core,  u^,  is  related  to  the  radial  stress,  by  the  following 

equation: 


"rc 


=  E, 


From  equation  (2)  it  follows  that 

dUp  _  A 
dr  r 


(5) 


Integration  of  the  above  equation  yields 

u^  =  A  log  r/a  +  B  (1+) 

The  use  of  log  r/a  instead  of  log  r  merely  alters  the  arbitrary  con¬ 
stant  of  integration. 


Equilibrium  of  Facings 

With  reference  to  the  differential  elements  of  the  facings  shown  in 
figure  It  is  seen  that,  when  the  cylinder  is  acted  upon  by  an  ex¬ 
ternal  loading  of  uniform  Intensity,  q,  the  only  forces  in  the  facings 
are  the  tensile  forces  per  unit  length  Ng  and  and  in  figure 
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3  represent  the  stresses  exerted  by  the  core  upon  the  outer  and  Inner 
facings,  respectively.  An  equilibrium  equation  can  be  obtained  for 
each  facing  by  summing  forces  In  the  radial  direction  on  each  differ¬ 
ential  element,  considered  to  be  of  unit  length  in  the  longitudinal 
direction.  The  equilibrium  equation  which  pertains  to  the  outer  facing 
is 


qa  dS  -  Oj^a  d0  -  ,Ng  d0  =  0 

This  equation  reduces  to 
Ng  =  a  (q  -  Oj.) 

or,  since  =  (a^.^) 


% 


-  (<^rc)  1 


(5) 


In  a  similar  manner,  the  equilibrium  equation  of  the  inner  facing  is 
found  to  be 


Ng  =  b 


(a  ) 

'  rc ' 


r  =  b 


(6) 


From  the  application  of  Hooke's  law  along  with  the  assumption  that  the 
tangential  stress  in  the  facings  is  uniformly  distributed  through  their 
thicknesses,  the  following  relationships  are  obtained; 

%  =  ^^o 

and 

Ng  =  Et^  eg.  (8) 


If  the  right-hand  sides  of  equations  (5)  and  (?)  and  the  right-hand 
sides  of  equations  (6)  and  (8)  are  equated,  the  following  two  equations 
result; 


(’rc)  1 

r  =  a  j 


and 


=  b 
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since  eg 


I 


U  ^  ,  U 

a  =  b 


these  two  equations  become 


and 


u 


(9) 


(10) 


If  continuity  of  displacements  at  the  interfaces  is  assumed,  then 

u  =  (u  ) 

'  c ' 

r  =  a 


and 


u 


b 


The  above  relationships  in  conjunction  with  equations  (2)  and  (1+)  enable 
equations  (9)  and  (lO)  to  be  written  as  follows: 


and 


I  A  log  b/a  +  B 


E  b 


(11) 


(12) 


Equations  (ll)  and  (12)  may  be  solved  for  A  and  B  with  the  following 
results: 


A  =  |^k 
Ec  • 

and 

B  =  p-  (1  -  k) 


(15) 


(14) 


where 


k 


1  t-b: 

a  t. 


1 _ 

Etr, 

log  b/a 


(15) 
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The  substitution  of  the  value  of  A  given  by  equation  (ij)  into  equation 
(2)  yields  the  following  expression  for  the  radial  stress  in  the  core: 

'^ro  =  <1  f  ^  (16) 

When  the  above  value  of  0^^  is  substituted  into  equations  (5)  and  (6), 
the  following  two  equations  result: 

Ng  =  qa  .(1  -  k)  (17) 


and 


Ng  =  qak  (I8) 

The  substitution  of  the  values  of  A  and  B  given  by  equations  (ij)  and 
(lU)  into  equation  (b)  yields  the  following  expression  for  the  core 
displacement: 


Et, 


c  Et^ 


[1  -  k  (1  -  :pr—  log  r/a)  ] 


Since  u  =  (u  )  and  u  =  (u  )  ,  then 

'^r  =  a  ‘^r=b 


and 


u-  =||^k 


(19) 

(20) 
(21) 


Equations  (16)  -  (2l)  completely  describe  the  stresses  and  displacements 
in  the  sandwich  cylinder  subjected  to  uniform  external,  lateral  loading. 
In  each  of  these  equations  the  value  of  k  is  given  by  equation  (15)- 


Stability  Analysis 


In  discussing  the  stability  of  the  sandwich  cylinder  under  uniform 
external  lateral  pressure,  the  equilibrium  of  a  slightly  deformed 
element  of  the  cylinder  must  be  considered.  It  is  assumed  that  the 
stress  situation  that  exists  in  this  deformed  element  differs  only 
slightly  from  the  stress  situation  that  existed  just  before  buckling. 

The  stresses  before  buckling  are  given  by  equations  (16),  (1T)»  and 
(18).  Since  the  cross  section  of  the  deformed  cylinder  is  no  longer 
circular,  the  small  changes  in  the  stresses  and  the  displacements  result¬ 
ing  from  buckling  will  be  functions  of  0  as  well  as  r.  Following  the 
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system  of  notation  used  in  the  origins.',  report^  a  har  is  placed  over  the 
appropriate  symbol  to  denote  the  small  stresses,  strains,  and  displace¬ 
ments  that  occur  when  the  cylinder  goes  from  the  initial  uniformly 
stressed  circular  form  to  the  slightly  deformed  configuration.  Again, 
it  is  assumed  that  the  core  extends  to  the  middle  surfaces  of  the  fac¬ 
ings  and  that  the  load  q  is  applied  to  the  middle  surface  of  the  outer 
facing.  This  assumption  is  now  somewhat  more  restrictive  fnan  it  was 
in  the  case  of  axial  symmetry,  since  it  now  means  that  the  effect  of 
the  interface  shear  on  the  bending  of  the  facings  and  the  displacements 
due  to  the  bending  of  the  individual  facings  are  neglected.  Hovrever, 
it  is  felt  that,  for  cylinders  with  sheii-type  facings,  the  results 
based  on  this  assumption  should  be  of  sufficient  accuracy. 


Equilibrium  of  the  Core 

Since  the  cylinder  is  now  considered  to  be  slightly  deformed,  the  core 
is  also  slightly  deformed.  A  differential  element  of  this  core  is  shown 
in  figure  h.  In  addition  to  the  radial  stress,  q  -p  k,  given  by  equation 

(l6),  a  small  radial  stress,  and  a  small  transverse  shear  stress, 

^rQc’  buckling  must  be  taken  into  account.  The  differential 

element  is  considered  to  be  in  equilibrium  under  the  action  of  the  stresses 
shown.  Since  the  small  change  in  the  geometry  of  the  core  introduces  only 
small  terms  of  higher  order  into  the  equilibrium  equations,  the  differ¬ 
ential  element  in  figure  k  is  shown  in  its  undeformed  state.  The  analysis 
of  the  core  is  exactly  the  same  as  that  given  in  the  original  report, 
and  for  this  reason  only  the  basic  equations  will  be  repeated  here. 


The  equilibrium  equations  that  apply  to  the  core  are  obtained  by  summing 
forces  in  the  radial  and  tangential  directions  in  figure  k.  These equi¬ 
librium  equations  are: 


^®rc  .  "^rc  ,  ^  ^"'’r0c 


=  0 


and 


^'•'rec  ^  ^'^rGc 
dr  r 


(22) 


(23) 


The  following  stress -displacement  relations  are  also  applicable; 


a 


re 


and 


^rOc  ^r0 


‘l  duc 

.r  dS  dr  r  J 


(2k) 

(25) 
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It  was  shown  in  the  original  report  that  the  small  radial  and  tangential 
displacements  of  the  core  can  be  completely  determined  Insofar  as  their 
dependence  on  r  is  concerned.  Then,  on  the  assumption  that  during  buckl¬ 
ing  the  circumference  of  the  cylinder  subdivides  into  n  waves,  the  dis¬ 
placement  functions  are  written  as  follows; 


cos  nS 


and 


J  +  C„  log  f] 

"  \  J  Sn  Bn  f  -  n  C„  (1  +  log  f )  +  f 


(26) 


sin  n0  (27) 


where 


E  2 

c  c  n 

n  2G  „  ■  2 
rS 


and 


A^,  Bj^,  C^,  and  are  arbitrary  constants, 


After  the  substitution  of  the  expressions  for  u^,  and  given  by  equations 

(26)  and  (27)  into  equations  (24)  and  (25),  the  expresFions  for  the  small 
core  stresses  become 


a 


rc 


cos  n0 


and 


■^rSc 


nr  ^  r 


sin  n0 


(28) 


(29) 


It  may  be  easily  verified  that  equations  (28)  and  (29)  satisfy  the  equi¬ 
librium  equations,  equations  (22)  and  (23). 


Equilibrium  of  Facings 


Figure  5  shows  the  differential  elements  of  the  facings  of  the  slightly 
deformed  cylinder.  In  addition  to  the  forces  that  exist  just  before 
buckling  (given  by  equations  (16),  (17),  and  (18)),  the  small  forces  and 
moments  that  arise  during  buckling  are  shown  in  the  figure.  In  consider¬ 
ing  the  equilibrium  of  the  facing  elements,  account  is  taken  of  the 
rotation,  and  stretching  of  the  facings  which  occurs  during  buckling.  As 
described  in  the  original  report,  the  initial  central  angle,  becomes 


1  8^u 
a  80^' 


d0 


and 


1  8v ' 
b  80 


1 

b  802 


d0 
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for  the  outer  and  inner  facings,  respectively,  and  the  areas  of  the 
differential  elements  of  the  outer  and  Inner  facings  become  (l  +  'eg)  a  d0 

and  (l  +  eg)  b  d0,  respectively.  Three  equations  of  equilibrium  can  be 

written  for  each  facing;  the  differential  elements  are  considered  to  be 
of  unit  length  In  the  longitudinal  direction. 


Considering  first  the  differential  element  of  the  outer  facing,  the  sum¬ 
mation  of  forces  in  the  direction  normal  to  the  surface  yields 

q  (1  +  eg)  a  d0  -  (qk  +  a^)(l  +  €g)  a  de 


_  /  1  ^v 


de  =  0 


—  u  1  dv 

If  the  relationship,  eg  “  a  a  §e’  used  and  small  quantities  of  higher 

order  --  that  is,  products  of  barred  quantities  --  are  neglected,  the  above 
equation  becomes 


3Qg  _  , 

Ng  -  •Kn  =  -  a  CTj,  +  qa  ( 


1  ^^u  \ 

1  -  k)  (  -  + - ) 

Va  a  / 


The  summation  of  forces  in  the  tangential  direction  yields  the  following 
equation:  _ 

r  _  1  r  _  ^Ng  T 

-  qa  (1  -  k)  +  NgJ  +  j^qa  (l  -  k)  +  Ng  +  ^g-  d0 


If  small  quantities  of  higher  order  are  again  neglected,  the  above  equa¬ 
tion  may  be  written  as 

cNfl  —  —  , ,,  , 

—  +  Qg  =  a  Tj.g  (51) 

The  third  equilibrium  equation  of  the  outer  facing  is  obtained  by  equating 
to  zero  the  summation  of  moments  about  point  0,  shown  in  figure  5*  Thus 


_  /_  SMg  N  _ 

-  Mg  +  ^Mg  +  d0y  +  Qg  a  d0  =  0 

The  above  equation  reduces  to 


^  +  a  Qg  =  0 
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Similarly  the  three  equilibrium  equations  which  pertain  to  the  inner 
facing  are  obtained  by  summing  forces  in  the  normal  and  tangential 
directions  and  by  summing  moments  about  the  point  O'  using  the  differential 
element  of  the  inner  facing  shown  in  figure  5.  These  equilibrium  equations 
are: 


ae 


—  t 

b  a 

r 


+  qak 


1  a^u’  ^ 

b  ' 


Sn' 

+  Qe  =  -  t 


and 


5^  +  b  Qe  =  0 


(35)' 

(34) 


(35) 


If  equation  (32)  is  solved  for  Qg  and  the  resulting  value  is  substituted 

into  equations  (30)  and  (3l)>  the  three  equilibrium  equations  of  the  outer 
facing  are  reduced  to  the  following  two  equations: 


and 


1 

a  502 


a  Oj,  +  q  (1  -  k)(  u  + 


aNg  1  aNIg 

Sr*  ■  I  50“  =  ^  ""r© 


(36) 


(37) 


Ii>  like  manner,  if  equation  (35)  is  solved  for  Qg  and  this  value  is  sub¬ 
stituted  into  equations  (33)  and  (34),  the  three  equilibrium  equations 
of  the  inner  facing  are  reduced  to  the  following  two  equations: 


and 


1  a^  ,  a  , 

Nq  + - —  =  b  cf  +  q  _  k 

®  b  502  b 


1  _ I 

50  b  50“  "  ■  ^  '^r0 


(38) 

(39) 


From  the  application  of  Hooke's  law,  the  following  tw9  expressions  relating 
the  small  tangential  forces  per  unit  length,  Ng  and  Ng,  to  the  small  tan¬ 
gential  strains,  Tg  and  eg,  are  obtained: 
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Et 


° 


and 


1  -  |i" 


N0  =  7— S  % 


(40) 


(41) 


1  -  ^ 

Also,  the  following  two  equations  relating  the  bending  moments  In  the 
facings,  Mg  and  Mg,  to  the  changes  In  curvature  in  the  facings,  5<g  andTCg, 
are  applicable:^ 


Et, 


Mg  = 


12(1  -  n^^) 


57  ^ 


and 


^0  ■  to/t  ..2', 


—I 

X, 


12(1  -  Q 


(42) 


(43) 


h 

Since— 


u  .  1  5v 

=  7  7  50 


u 

.  1 

dv' 

b 

b 

50 

/5v 

d^u 

■P 

^0 

d02 

1 

/dv' 

d^u' 

b^' 

U0 

■  d02 

and 


equations  (40)  -  (43)  may  be  written  as  follows: 
_  Etr,  1  5v\ 


_  Etf,  /u  1  dv\ 


(44) 


■Timoshenko,  S.,  and  Goodler,  J.  Theory  of  Elasticity.  New  York,  1951* 

Timoshenko,  S.  Theory  of  Plates  and  Shells,  First  Edition.  New  York 
1940. 
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Ng  = 


Etj  ^  1 

- _  /  _  + - 


1  -  H  Vb  b  Se 


(^+5) 


and 


Mg  =  - 


Mg  =  - 


Eto' 


12  (1  -  V^e  aev 


^2^ 


■ov  d  u\ 


Et.' 

1 


12  (1  -  M^)  b^  '^he  be 


- 1  'kS— t 

'V  o  u 


(46) 


(47) 


The  requirement  is  now  made  that  there  be  continuity  of  displacements  at 
the  interfaces;  t}iat  is. 


n  =  (uj 


u  =  (u J 


=  a 


r  =  b 


V  =  (vj 


(48) 

(49) 

(50) 


r  =  a 


and 


"  =  (vj 


(51) 


The  use  of  equations  (t8)  -  (.51)  enables  equations  (44)  -  (hj)  to  be 
written  as  follows: 


Et 


N0  = 


2 _ ft 

.2^  1. 


^0 


(1  -  b^) 

Et 


be  ) 


r  =  a 


.  ,  bv 

1  ,  c 

u  +  — 


b  (1  -  q^)  \  be 


Et 


and 


^0  =  - 


4  =  - 


r  =  b 

5  b^ 

_ (  c _ c 


12a^  (1  -  p^)  \80  J ^  ^ 


Etj^ 


12b 


d2^ 


flc 


(1  -  p2)  V80  be^  ) 


(52) 


(55) 


(54) 


(55) 


r  =  b 
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(56) 


Et-r 

^  1 

Etf^ 

1  X  *  1 

b  (1  -  h^) 

Vde  be^  / 

12  (1  -  ^2)  b3  ' 

r  =  b 


=  -  b  (t  .  ) 


r  =  b 

(63) 


r  =  b 


If  the  expressions  given  by  equations  (26)  -  (29)  for  the  core  displace¬ 
ments  and  stresses  are  substituted  into  equations  (60)  -  (63),  four  equa¬ 
tions  containing  the  parameters  Aj^,  Bj,,  Cj,,  and  are  obtained.  These 
four  equations  are: 

[(n^  -  1)  7  (1  -  k)  -  (n^  -  1)]  A  +  [(n^  -  l)  7  (l  -  k)  +  5  (l  +  n^  *  ) 

n  no 

+  (1  +  n‘+  -  3  ]  Bn  +  [0  -  n2  (1  +  n2 

+  [n  (1  +  n^  41q)]  Hn  =  0  (64) 

(n^-l)A  +  [-  5  (l+'>)-(l  +  n^i>)+^]B 

'  '  n  ^  n  '  o'  '  o'  ■'  n 

+  [n2  (1  +  O  )]  C  +  [  -  n  (1  +  K  )]  H  =0  (65) 

o'  n  o'  n  '  ' 


[  (n*^  -  1)  7k  ^  -  (n^  -  1) )  A^  +  [(n*^  -  1)  ^  7k  6^  |(l  +  n"^  't'^) 

+  I  (1  +  n**  i*^)  +  3  B^  +  [(n^  -  1)  7k  ^  log  b/a  "  P  ^  ^ 

-  (n^  -  1)  log  b/a  -  rfi  (l  +  n^  ♦^)]  Cj,  +  [n  b/a  (i  +  n^  it*^)]  11^  =  0 


(„2  .  l)  A„  +  (  -  6„  I  (1  *  .1)  -  I  (1  +  n2  .,)  -  ^  1  B„ 

+  [(n^  -  1)  log  b/a  +  (1  +  ♦j^)]  +  [-  n  b/a  (l  +  '>^)]  11^  =  0 

(6?) 

where,  in  each  of  the  above  equations, 

2\ 
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"  2G 


r0 

2 


128*^ 


and 


♦l  = 


12b^ 


Et„ 


Each  of  the  terms  in  equations  (64)  -  (67)  contains  one  of  the  para¬ 
meters  Ajj,  Cfj,  and  that  appear  in  the  displacement  equations  of 

the  cylinder.  A  buckled  form  of  equilibrium  is  possible  only  if  equa¬ 
tions  (64)  -  (67)  yield  solutions  for  these  parameters  which  are  differ¬ 
ent  from  zero.  This  requires  that  the  determinant  of  the  coefficients 
of  these  parameters  must  be  equal  to  zero.  This  determinant  may  be 
written  as  follows: 
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The  use  of  equations  (l6)  throi;gh  (2l)  Is  self-explanatory;  if  the  dimen¬ 
sions  and  material  properties  of  a  given  sandwich  cylinder  are  known,  the 
stresses  and  displacements  in  terms  of  the  external  load,  £,  may  be  easily 
computed.  In  the  case  of  sandwich  cylinders  having  facings  of  equal  thick¬ 
ness,  equations  (l6)  -  (2l)  reduce  to  equations  (l8)  -  (25)  of  the  original 
report.  It  is  of  interest  to  compare  equations  {16)  -  (2l)  with  results 
obtained  by  Reissner.-^  If  Belssner's  equations  are  expressed  in  the  notation 
used  here,  the  following  equations  for  the  cylinder  stresses  result: 


Nfl  =  qn 


^  qa 


1  +  (1  -  b/a) 


^  1  -  b/a) 


2  +  ^  (1  -  b/a) 
L  a 
c 


(at  middle  surface  of  core)  =  q 


g  ^ 

1  +  b/a/  . 


(1  -  b/a) 


Equations  (16)  -  (18),  for  facings  of  equal  tfiickness,  may  be  written  as 
follows; 

I  r  -  b/a  1 


1%  =  qa 


Nfl  =  qa 


a  E^a  ' 

1  +  b/a  -  —  log  b/a 


1  +  b/a  -  log  b/a 

bc^ 


a_  (at  middle  surface  of  core)  =  q 


2 

1  +  b/a 


1  +  b/a  -  —  log  b/a 


Since  the  first  term  in  the  series  expansion  of  log  b/a  is  -  (l  -  b/a),  it 
may  be  noted  that,  for  cylinders  with  b/a  ratios  close  to  1,  the  equations 
of  Reissner  yield  results  that  are  very  nearly  the  same  as  those  given  by 
equations  (16),  (I7),  and  (18),  Relssner's  results  are  based  on  the  assump¬ 
tion  that  a  -  b  <  <  a  and  hence  become  increasingly  less  accurate  as  the 
cylinder  thickness  is  Increased. _ 

eissner,  Eric,  Small  Bending  and  Stretching  of  Sandwich-type  Shells. 
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where 


- 

2G^g  b  (1  -  i^) 

Equal i.,  1  (70),  with  n  =  2,  yields  values  of  within  5  percent  of  the 
values  obtained  from  equation  (69)  for  usual  sandwich  constructions. 

For  cylinders  having  very  thin  facings  (membrane  facings)  of  equal  thick¬ 
ness,  4’q  and  are  assumed  to  be  zero,  and  equation  (70)  reduces  to 


7  =  -  (n" 

'cr 


\  n2  Et  (1  - 

1  +  t  \  1  +  - ^ ^ 


2G_e  b  (1  - 


Equation  (71)}  for  n  =  2,  becomes 


r. 

+  -p)  1  + - ; 


The  value  of  q^^  is  then  determined  from  the  definition  previously  given, 

-  _  Et 

acr - - - oT  /cr 

a  (1  -  q2) 

Equation  (72)  may  be  used  as  a  good  approximation  to  the  expression  ob¬ 
tained  in  the  original  report  (equation  64)  for  the  critical  load  on 
cylinders  with  membrane  facings. 

It  is  of  interest  to  examine  the  results  given  by  equation  (68)  for  certain 
limiting  cases  other  than  that  of  membrane  facings.  Some  of  these  results 
are  given  in  table  1.  In  each  of  these  limiting  cases,  and  are 

neglected  as  compared  to  1. 

Equation  (68)  may  be  used  for  determining  the  approximate  value  of  the 
critical  load  on  long  sandwich  panels  in  the  form  of  a  portion. of  a  cylinder 
hinged  along  the  edges  0=0  aixi  0  =  2a  as  shown  in  figure  (7)-  If,  in 
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equation  (68),  it /a  is  substituted  for  n,  the  smaller  absolute  value  of 
£  obtained  from  equation  (68)  represents  the  critical  load  on  a  panel 
whose  dimensions  and  properties  are  known.  This  solution  applies  to  the 
unsymrnetrlcal  tjrpe  of  buckling  shown  in  figure  (7).  If,  as  in  the  case 
of  relatively  flat  panels,  symmetrical  buckling  with  no  inflection  point 
between  the  supports  occurs,  equation  (68)  is  not  applicable. 


Conclusions 


The  purpose  of  this  report  was  to  extend  the  previous  work  done  in  con¬ 
nection  with  sandwich  cylinders  subjected  to  uniform  lateral  loading  by 
taking  into  account  the  effect  of  the  stiffnesses  of  the  individual  fac¬ 
ings  of  the  cylinder  and  also  by  making  the  results  applicable  to  sand¬ 
wich  cylinders  having  facings  of  unequal  thickness.  The  results  indicate 
that  for  the  majority  of  sandwich  cylinders  the  analysis  based  on  membrane 
facings  is  sufficiently  accurate.  However,  for  cylinders  having  relatively 
thick  facings,  say  in  the  neighborhood  of  one-fourth  of  the  core  thickness, 
the  facing  stiffnesses  have  an  appreciable  effect  on  the  critical  load  on 
the  composite  cylinder.  Since  the  magnitude  of  the  effect  of  the  facing 
stiffnesses  is  dependent  not  only  on  the  dimensions  of  the  cylinder  but 
also  oh  the  mechanical  properties  of  the  core  and  facing  materials,  equa¬ 
tion  (68)  of  this  report  should  be  used  for  computing  the  critical  load 
if  Indications  are  that  this  effect  may  be  important.  In  all  cases  of 
sandwich  cylinders  having  facings  of  unequal  thickness,  the  equations  of 
this  report  should  be  used  for  computing  stresses  prior  to  buckling  and 
for  the  deteiminatlon  of  critical  loads. 
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Figure  5- --Differential  elements  of  inner  and  outer  facings  of  uniformly 


Figure  4. --Differential  element  of  core  of  slightly  deformed  cylinder 
neglecting  changes  In  geometry  of  core. 


OUTER  FACING 


Figure  6, --Configurations  corresponding  to  the  two  eigenvalues 
of  the  stability  determinant. 


BUCKLING  OF  SANDWICH  CYLINDERS  OF  FINIIE  LENGTH 
UNDER  UNIFORM  EXTERNAL  LATERAL  PRESSURE* 

By 

MILTON  E.  RAVILLE,  Engineer 

Forest  Products  Laboratory,  Forest  Service 
U.  S.  Department  of  Agriculture 


Summary 

A  theoretical  analysis  is  made  of  the  problem  of  the  buckling  of 
circular  cylinders  of  sandwich  construction  acted  upon  by  uniform  ex¬ 
ternal  lateral  pressure.  The  solution  obtained  is  based  on  the  assump¬ 
tion  that  the  sandwich  cylinder  is  comprised  of  isotropic,  membrane 
facings  and  an  orthotrop-ic  core.  The  mathematical  solution  of  the 
problem,  which  is  in  the  form  of  a  characteristic  determinant  of  sixth 
order,  is  applicable  to  sandwich  cylinders  of  any  length  and  of  any  core 
thickness.  Numerical  results  are  obtained  for  various  values  of  the 
parameters  that  enter  the  problem,  and  cuir/eB  are  included  which  illus¬ 
trate  how  ^he  critical  load  varies  as  the  values  of  these  parameters  are 
varied. 


*This  report  is  one  of  a  series  prepared  by  the  Forest  Products 

Laboratory  under  U.  S.  Navy,  Bureau  of  Aeronautics  Order  No.  01 593- 
♦^^aintained  at  Madison,  Wis.,  in  cooperation  with  the  University  of 
Wisconsin. 
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Introduction 


Sandwich  construction,  as  it  is  usually  employed,  consists  of  two 
relatively  thin  sheets  of  a  material  that  has  comparatively  high 
strength  and  stiffness  properties  separated  by  and  bonded  to  a  relative¬ 
ly  thick  layer  of  a  lightweight  material  that  has  comparatively  low 
strength  and  stiffness  properties.  The  two  outer  sheets  are  commonly 
called  "facings,"  and  the  inner  layer  is  commonly  called  the  "core." 
Sandwich  construction  results  in  a  composite  structural  element  having  a 
higher  strength -weight  ratio  than  can  be  obtained  through  the  use  of  a 
single  homogeneous  material.  For  this  reason,  its  major  fields  of  ap¬ 
plication  are  in  structures  in  which  weight  is  a  prime  factor,  such  as 
aircraft  and  guided  missiles.  However,  the  development  of  improved  ad¬ 
hesives  and  the  Improvement  of  methods  of  fabrication  are  making  its  use 
practical  in  an  Increasingly  large  number  of  different  types  of 

structures . 

I 

With  the  Increased  use  of  sandwich  construction  has  come  an  in¬ 
creased  need  for  reliable  design  data.  Much  work  of  both  a  theoretical 
and  an  experimental  nature  has  been  done  in  an  effort  to  provide  this 
needed  information.  The  physical  properties  of  many  of  the  different 
materials  used  in  sandwich  construction  have  been  determined,  and  analy¬ 
ses  of  a  number  of  stress  and  stability  problems  have  been  performed. 
Progress  has  been  somewhat  retarded  because  of  the  many  variables  that 
must  be  taken  into  account  in  sandwich  analysis  and  the  fact  that  the 
analysis  of  a  layered  system  is  Inherently  a  difficult  problem.  A  rath¬ 
er  complete  bibliography  of  the  part  of  this  work  that  pertains  particu¬ 
larly  to  shells  and  shell -like  structures  is  contained  in  a  recent  U.  S. 
Navy  publication  (^) . 
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The  purpose  of  this  thesis  Is  to  obtain  the  solution  for  the 
critical  load  on  circular  sandwich  cylinders  of  finite  length  subjected 
to  uniform  external  lateral  pressure.  The  critical  load  Is  defined  as 
the  intensity  of  pressure  at  which  the  cylinder  buckles  due  to  lack  of 
stiffness.  This  particular  problem  was  chosen  becairse  it  represents  a 
fundamental  problem  In  sandwich  analysis  that  has  not  previously  been 
solved.  The  need  for  such  an  analysis  has  arisen  in  connection  with  the 
design  of  certain  component  parts  of  aircraft  and  guided  missiles.  It 
Is  believed  that  the  method  of  analysis  used  here  is  somewhat  more  rig¬ 
orous  than  methods  that  have  been  used  to  obtain  solutions  to  other 
problems  in  sandwich  analysis.  For  this  reason,  the  general  method  of 
approach  may  be  of  use  In  obtaining  solutions  to  other  new  problems  as 
well  as  in  obtaining  better  solutions  to  some  problems  that  have  already 
been  solved.  The  solution  contained  in  this  thesis  represents  an  exten¬ 
sion  of  a  previous  solution  of  the  problem  of  the  buckling  of  long  sand¬ 
wich  cylinders  (2) . 

In  this  thesis  It  is  assumed  that  the  sandwich  cylinder  is  com¬ 
prised  of  very  thin,  Isotropic  facings  of  a  relatively  stiff  material 
separated  by  an  orthotropic  core.  The  analysis  of  the  facings  is  based 
on  membrane  theory.*  The  core  Is  considered  to  have  such  a  low  load¬ 
carrying  capacity  in  the  tangential  and  longitudinal  directions  as  com¬ 
pared  to  the  facings  that  the  normal  stress  In  the  core  in  these  direc¬ 
tions  and  the  shear  stress  In  the  core  on  planes  perpendicular  to  the 
facings  and  in  these  directions  may  be  neglected.  This  core  assumption 

*An  analysis  taking  into  account  the  stiffnesses  of  the  individual  fac¬ 
ings  may  be  made  along  the  lines  illustrated  in  Ref.  5* 
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has  been  widely  used  In  sandwich  analysis  and  is  known  to  represent  ac¬ 
tual  seuidwlch  construction  very  well.  The  fact  that  no  further  simpli¬ 
fying  assumptions  in  regard  to  the  core  are  needed,  other  than  the  as¬ 
sumption  that  it  behaves  as  an  elastic  continuum,  is  emphasized.  It  is 
felt  that  some  analyses  of  sandwich  construction  are  not  sufficiently 
accurate  in  certain  ranges  of  physical  properties  and  dimensions  because 
of  additional  simplifying  assumptions  made  in  regard  to  the  core.  The 
action  of  the  core  and  facings  is  related  by  the  assumption  that  their 
displacements  are  equal  at  the  interfaces  between  the  core  and  facings. 
These  interfaces  are  assumed  to  be  at  the  middle  surfaces  of  the  facings. 

The  method  used  for  establishing  the  stability  criterion  is  simi¬ 
lar  in  concept  to  that  used  by  Timoshenko  in  the  .analysis  of  the  buck¬ 
ling  of  homogeneous  cylinders  of  finite  length  subjected  to  uniform  ex¬ 
ternal  lateral  pressure.*  The  assumption  is  made  that,  for  pressures 
less  than  the  critical,  the  circular  cylinder  remains  circular  and  the 
only  stresses  present  are  a  uniform  circumferential  stress  in  the  fac¬ 
ings  and  a  uniform  radial  stress  in  the  core.  In  discussing  the  buck¬ 
ling  of  the  sandwich  cylinder,  only  small  deflections  from  this  uniform¬ 
ly  compressed  form  of  equilibrium  are  considered;  thus,  the  stresses  in¬ 
duced  in  the  cylinder  as  it  goes  from  the  circular  to  the  slightly  de¬ 
formed  configuration  may  be  considered  small  as  compared  to  the  pre- 
buckllng  stresses 


*See  Ref.  4,  Art.  83. 
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Notation 


r,  9,  z 


radial,  tangential,  and  longitudinal  coordinates, 
respectively 

radius  to  middle  surface  of  outer  facing 

radius  to  middle  surface  of  inner  facing 

thickness  of  each  facing 

length  of  cylinder 

modulus  of  elasticity  of  facings 

Poisson's  ratio  of  facings 

modulus  of  rigidity  of  facings 

modulus  of  elasticity  of  core  in  direction  normal  to 
facings 

modulus  of  rigidity  of  core  in  r£  plane 
modulus  of  rigidity  of  core  in  rz  plane 
intensity  of  uniform  external  lateral  loading 


Et  log  k 

1+1 - - - - 

a  a 


small  normal  stress  in  core  in  radial  direction 


r0'  rz 


u  ,  V  ,  w 
c  c  c 


small  transverse  shear  stresses  in  core 

small  reidial,  tangential,  and  longitudinal  core  dis¬ 
placements,  respectively 


e  ,  ,  e 

r  e  z 
c  c  c 


’'re  '  ’rz  '  ’ez, 
c  c  c 


small  radial,  tangential,  and  longitudinal  core 
normal  strains,  respectively 

small  radial,  tangential,  and  longitudinal  core 
shearing  strains,  respectively 


number  of  waves  in  circumference  of  buckled 
cylinder 
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X 


xa 

i 


Eq  n2 


small  normal  strains  In  tangential  direction  In 
outer  and  inner  facings,  respectively 


e 

z 


> 


€' 

Z 


small  normal  strains  In  longitudinal  direction  In 
outer  and  Inner  facings,  respectively 


u,  V,  w 


small  radial,  tangential,  and  longitudinal  dis¬ 
placements  of  outer  facing 


u' ,  V ,  w' 


small  radial,  tangential,  and  longitudinal  dis¬ 
placements  of  Inner  facing 


N  ,  n 

6  z  Bz 


normal  forces  and  shear  force  i)er  unit  length  of 
outer  facing 


3 


a 


normal  forces  and  shear  force  per  unit  length  of 
Inner  facing 

E  a  (1  -  u2) 
c 

Et 

qa  (1  - 
Et 


log 


natural  logarithm 


arbitrary  constants 
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Mathematical  Analysis 


In  the  analysis  of  the  sandwich  cylinder,  cylindrical  coordinates 
£,  £,  and_z  are  used;  the  dimensions  of  the  cylinder  and  the  positive 
directions  of  the  coordinates  are  Indicated  In  figure  1.  The  radii  to 
the  middle  surfaces  of  the  outer  and  Inner  facings  are  denoted  by  a  and 
b,  respectively,  aind  the  thickness  of  each  facing  Is  denoted  by  t.  The 
origin  of  the  coordinate  system  Is  placed  at  the  middle  cross-section  of 
the  cylinder  whose  length  Is  denoted  by  Since  £  is  considered  to  be 
positive  when  It  acts  in  the  positive  r-directlon,  buckling  occurs  at  a 
negative  value  of 

As  mentioned  previously,  the  cylinder  at  the  insteint  before  buck¬ 
ling  Is  assumed  to  be  In  a  state  of  uniform  compression  with  the  same 
stress  condition  existing  in  each  cross-section.  This  assumption  amounts 
to  saying  that  the  ends  of  the  cylinder  are  not  supported  until  after 
the  cylinder  has  been  initially  compressed  by  a  uniform  load  Just  less 
than  the  critical  load.  The  expressions  for  the  pre-buckling  stresses 
are  known,*  and  they  are  shown  in  figure  2.  The  radial  stress  in  the 
core  Is  equal  to  q  S  k  and  the  circumferential  force  per  unit  length  in 
the  outer  and  inner  facing  is  qa  (l  -  k)  and  qak,  respectively,  where  £ 

Is  the  intensity  of  uniform  external  lateral  loading  in  the  positive  r 
direction  and 

k  =  _ i _ 

b 


b  Et  log  i 
a  ■  Ec  a 


*See  Ref.  2. 
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In  discussing  the  equilibrium  of  the  cylinder  after  it  has  buckled  into 
a  slightly  deformed  shape,  it  is  assumed  that  the  existing  stress  con¬ 
dition  differs  only  slightly  from  the  stress  condition  shown  in  figure  2. 
The  analysis  of  the  equilibrium  of  a  differential  element  of  the  core  of 
the  deformed  cylinder  and  of  differential  elements  of  the  facings  of  the 
deformed  cylinder  is  considered  next. 


Equilibrium  of  the  Core 

A  differential  element  of  the  deformed  core  is  shown  in  figure  3. 

The  stresses,  0  ,  t  ,  and  t  ,  resulting  from  buckling  are  assumed  to 
r  rQ  rz 

be  small  as  compared  to  the  pre-buckllng  stress,  q  ^  k.  As  indicated 

previously,  the  assumption  is  made  that  the  core  and  facing  materials 

are  such  that  the  stresses,  a  ,  a  ,  and  t  ,  in  the  core  may  be  neglected. 

d  Z  V  z 

The  effect  of  the  changes  in  geometry  in  the  core  was  found  to  be  negli¬ 
gible;  therefore,  for  simplification,  the  equilibrium  equations  are  writ¬ 
ten  on  the  basis  of  the  original  geometry  of  the  element.  The  summation 
of  forces  in  the  radial  direction  results  in  the  following  equilibrium 
equation; 


d  q  S  k  <^0 

(.  q  5  k  -  0  )  r  d0  dz  +  (q  f  k  +  Op  +  - - —  dr  +  v—  dr) 

r  r  ^  dr  dr 


d  T  d  T 

(r  +  dr)  de  dz  +  -  d0  dr  dz  +  — : — ^  (r  +  dr)  dr  d9  dz  =  0 

de  ^ 


If  terms  containing  the  product  of  more  than  three  differentials  are  neg¬ 
lected,  the  above  equation  may  he  written  as 


1  ^^r9 

r  de 


+ 


0 


(1) 
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The  sunimatlon  of  forces  In  the  tangential  direction  yields 


St 

T  r  d9  dz  +  (t  +  ^  dr)  (r  +  dr)  d9  dz  +  t  -  dr  d0  dz  -  0 

Tu  ^  o  r  ^ 


or,  neglecting  differentials  of  higher  order, 


Sr  r 


(2) 


Finally,  the  summation  of  forces  in  the  longitudinal  direction  yields 

St 

-T  rdSdz  +  (T  +  — dr)  (r  +  dr)  dS  dz  =  0 
rz  rz  Sr 


or 


(3) 


If  it  is  assumed  that  during  buckling  the  generators  of  the  cyl¬ 
inder  deflect  into  a  half-wave  of  a  sine  curve  and  the  circumference 
subdivides  into  n  waves,  the  core  displacements  take  the  following  form: 


u 

c 

=  f^{r) 

COS 

n0 

cos 

«z 

£ 

V 

c 

=  f2(r) 

Bin 

n0 

cos 

nz 

£ 

w 

c 

cos 

n0 

sin 

nz 

£ 

where  u^,  v^,  and  w^  represent  the  small  displacements  of  the  core  in 
the  r,  S,  and  z  directions,  respectively.  At  the  ends  of  the  cylinder 
S^u 

u  euid  — ^  are  both  zero,  which  represents  the  conditions  of  simply- 
—  Sz^ 
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supported  edges.  The  displacements  u^,  v^,  and  are  related  to  the 
core  strains  by  the  following  equation^ 


du  dv  V 


bu  bw 


(5) 


dv 

( 

dz 


Also,  since  Cq  and  are  assumed  to  be  zero. 


du 

a_  -  E^e  =  E„ 

^c  ar 


In  addition  to 


and 


au  aw 

T  =  G  7  =0 

rz  rz  rz^  rz  'dz  dr 


(6) 


(7) 


(8) 


If  the  expressions  for  the  displacements  given  by  equations  (U)  are  sub¬ 
stituted  into  equations' (6),  (7),  and  (8)  and  the  resulting  expressions 
for  the  stresses  are  required  to  satisfy  the  equilibrium  equations (l) , 


*See  Ref.  5*  P-  505* 
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(2),  and  (3),  the  functions  of  r  that  appear  In  equations  (4)  may  be 
determined  explicitly,  thus  enabling  equations  (4)  to  be  written  as 
follows : 


^c  “  ^  ®n®^  +  "^n®  ^08  §  +  Dj^r]  cos  nS  cos  J2. 

A  2 

^c  =  ^  "  V  -  “  V  (1  +  s)  + 

n  D  r  log  £  ■(■  L  r]  sin  nfl  cos  ££ 
nan  / 


(9) 


(10) 


Wc  “  log  £  +  X  Bj,r  +  X.  C^^r  (log  £-!)+>. 

^  log  |)  +  RjjS]  cos  n0  sin  ^ 

where  X  =  ££ 
i 


and  A  ,  B  ,  C  ,  D  ,  L  ,  and  R  are  arbitrary  constants, 
n  n'  n  n  n  n 


(11) 


If  equations  (9),  (lO),  and  (ll)  are  substituted  into  equations  (6),  (7), 
and  (8),  the  expressions  for  the  core  stresses  become 


«r  "®c 


(-  2_+C  ^  +  D  )  cos  n0  cos  — 

'  "n  ^2  "  r  n'  / 


(12) 
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(15) 


A  2 

T  =  -  E  —  ^  sin  n9 
rd  c  n  2 


cos 


■nz 

I 


D 

T  =  -  E  —  —  cos  n0  sin  (l^) 

rz  c  \  r  z  \  ' 

i 

It  may  be  easily  verified  that  equations  (l2),  (13),  and  (lU)  satisfy 
the  equilibrium  equations,  equations  (l),  (2),  and  (5)«  Furthermore, 
because  of  the  manner  in  which  they  were  derived,  it  may  be  stated  that 
the  right-hand  sides  of  equations  (9)>  (lO),  and  (ll)  are  unique  func¬ 
tions  insofar  as  their  dependence  upon  r  is  concerned. 


Equillbrluoi  of  the  Facings 

Differential  elements  of  the  facings  of  the  deformed  cylinder 
are  shown  in  figure  1.  The  quantities  N.,  n',  N  ,  n',  N.  ,  and  N'  ' 

u  p  Z  Z  uZ  uZ 

represent  normal  and  shear  forces  per  unit  length  of  the  facing  upon 

which  they  act.  The  core  is  assumed  to  extend  to  the  middle  surfaces 

of  the  facings,  and  the  stresses  (a  )  ,  (a  )  ,  ,  (t  „)  ,  (t  )  , 

r'r^Si’  ^  rg'^r=a^  ^  rz'^r=a' 

(t  _)  ,  ,and  (t  )  ,  are  the  stresses  exerted  by  the  core  on  the  fac- 

^  rg^r=b^  rz'^r=b 

Ings  as  a  result  of  buckling.  All  of  the ’above  quantities  are  assumed 
to  be  small  in  comparison  to  the  pre-buckling  stresses  shown  in  figure 
1  as  functions  of  the  intensity  of  loading,  a.  The  bending  moments  and 
Transverse  shear  forces  in  the  individual  facings  are  neglected,  in  ac¬ 
cordance  with  membrane  theory.  In  formulating  the  equilibrium  equations 
that  apply  to  the  facings,  the  effects  of  the  rotation  and  stretching  of 
the  facings  must  be  taken  into  account.  Because  of  the  stretching  of 
the  facings,  the  areas  of  the  differential  elements. of  the  outer  and 
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Inner  facings  become  (l  +  e^)  (l  +  e^)  a  d0  dz  and  (l  +  e^)  (l  +  e^) 
b  d0  dz,  respectively.  Since  the  middle  surface  strains  are  small  quan¬ 
tities,  their  products  may  be  neglected,  and  the  differential  areas  of 

the  outer  and  inner  facings  are  expressed  as  (l  +  e  +  e  )  a  d0  dz  and 

_ 8  z _ 

(l  +  e'  +  e '  )  b  de  dz,  respectively.  Tiie  difference  in  the  amount  of 
_ 9  z _ 

rotation  of  the  sides  ^  and  CD  with  respect  to  the  £-axls  of  the  cyl¬ 
inder  causes  the  central  angle  of  the  outer  element  to  change  from  dS 


to  (l  +  —  - ii)  d0.  Similarly,  the  central  angle  of  the  inner 

a  00  a  ^  2 

be 

- .  2 

element  changes  from  ^  to  (l  +  ^  ^  dS  •  Also,  the  differ- 

be^ 


ence  in  the  amount  of  rotation  of  the  sides  ^  and  CD  with  respect  to 
2 

the  r-axis  is  (^7  +  ;  the  corresponding  difference  in  rotation  of 

d0dz  dz 

2 

the  sides  A'B'  and  CD'  of  the  inner  element  is  ( ■-  ^ -  +  ^H_)  .  The 
-  -  bebz  bz 


foregoing  expressions  for  the  differences  in  rotation  are  as  given  by 
Timoshenko*’  and  verified  by  Osgood  and  Joseph.-**  Since  and  are 
small  quantities,  the  relative  rotation  of  the  sides  upon  which  they  act 
may  be  neglected;  the  Inclusion  of  this  effect  leads  only  to  small  terms 
of  higher  order. 

The  necessary  changes  in  the  geometry  of  the  facings  having  been 
established,  it  is  now  possible  to  write  three  equilibrium  equations 
for  each  facing.  Considering  first  the  differential  element  of  the 
outer  facing,  the  summation  of  forces  in  the  direction  perpendic\ilar  to 
the  surface  of  the  exement  results  in  the  following  equation: 

■*See  Ref.  4,  Art,.  79. 

**See  Ref.  6. 
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-  (qk  +  ®2^  a  d9  dz  +  q  (l  +  +  e^)  a  d0  dz 

-  [qa  (1  -  k)  +  R.]  {1  +  €,)  (1  +  1  ^  -  1  d©  dz  -  0 

®  2  a  ae  a  ^2 

If  the  relationships  €g  =  H  +  i  ^  and  ^  are  used  and  terms  con¬ 

taining  the  products  of  the  small  stresses  and  displacements  are  neg¬ 
lected,  the  above  equation  of  equilibrium  reduces  to 

Ng  =  -  a  +  q  (1  -  k)  (u  +  (15) 

The  summation  of  forces  in  the  tangential  direction  perpendicular  to  the 
rz  plane  yields 

SN 

-  [qa  (l  -  k)  +  Ng]  (1  +  €^)  dz  +  [qa  (l  -  k)  +  ^  dP] 

SN- 

(1  +  ej  dz  -  (1  +  €g)  a  d0  +  ^  dz)  (l  +  c^) 

a  d0  -  (1  +  +  G^)  a  d0  dz  =  0 


If  terms  containing  products  of  the  small  stresses  and  strains  are  neg¬ 
lected,  the  above  equation  reduces  to 


1  ''"0  ^ 

150"  + 


r0 ' r=a 


The  summation  of  forces  in  the  tangential  direction  parallel  to  the 
z-sucls  yields 
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&N 

-  N^(l  +  Eg)  a  d0  +  (N^  +  (1  +  Eg)  a  dfl  -  Ng^  (l  +  E^)  dz 

dN 

*  <"»»  *  (1  *  .,)  1«  (1  *  -s  *  s) 


2 

a  dfl  dz  -  qa  (l  -  k)  ^  d9  dz  ■  0 

ddoz  oz 


If  small  quantities  of  higher  order  are  again  neglected,  the  above  equa¬ 
tion  of  equilibrium  reduces  to 


^’’z  ^  1  ^"ez 

a  ~5fl~ 


(^rz^r=a  *  <l 


/d^v 

5s5r 


+ 


du\ 

51^ 


(17) 


The  three  equilibrium  equations  pertaining  to  the  inner  facing  are  ob¬ 


tained  in  a  similar  manner  and  may  be  written  as  follows; 


0  r  r=b  b  ^2 


b  de  dz  re^r=b 


and 


dN' 

z 


dz  b  dS 


)  „  +  q  2.  k 

rz^r=b  ^  >,  ^ 


(18) 

(19) 


(20) 


On  the  basis  of  Hooke's  law,  the  following  stress-strain  relations  in 
the  facings  are  applicable: 


Et 

l-p2 


(e„  +  pE^) 


(21) 
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Et 


(23) 


=  Gt  7^ 
9z  9z 


2(l-Hi)  62 


^9  = 


Et 


1-U 


2  ^^9  ^  ^^z^ 


(24) 


1-4 


2  ^  ^^9'> 


and 


Et 


^flz  ^9z  2(l-Hi)  '6z 


71. 


(25) 


(26) 


where  E,  G,  and  ji  represent  the  modulus  of  elasticity,  the  modulus  of 
rigidity,  and  the  Poisson's  ratio  of  the  facings,  respectively.  Since 
the  strain-displacement  relations  given  by  equations  (5)  are  applicable 
to  the  facings  if  r  is  replaced  by  a  for  the  outer  facing,  r  is  replaced 
by  b  for  the  inner  facing,  and  the  corresponding  displacements  of  the 
facings  are  substituted  for  the  core  displacements,  equations  (21) 
through  (26)  may  be  written  as  follows: 


N 

9 


Et 

o 

1-4“ 


(H  +  1  ^  +  n 

a  a  de  dz 


(27) 


(28) 
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a 

N'  =  b  (o  )  ,  +  q  §  k  (u  +  - 2.) 

0  r'r»b  ^  b  '  c  ^2  ■'r=b 


(U2) 


b  d0  da 


- 


rd'r=b 


(i*?) 


end 


dz  b  d6 


du 


^'^rzKr=h  *  S  ^  ^dedz  *  dz 


(hk) 


If  the  expressions  for  the  facing  stresses  given  by  equatlcns  (33)  through 
(5Q)  are  substituted  Into  the  above  equations  and  then  a  further  substi¬ 
tution  is  made  of  the  expressions  for  the  core  displacements  and  stresses 
given  by  equations  (9)  througlj  (lU),  six  equations  containing  the  para¬ 
meters  A^,  B^,  C^,  D^,  L^,  and  that  appear  in  the  equations  for  the 
displacements  of  the  core,  equations  (l*),_are  obtained.  These  equations 
may  be  written  as  follows: 
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C\J  lOJ 
^1  (d 


OJ  Icy 

^1  <a 


OJ  |OJ 

.o|  (d 


OJ  lOJ 
^1  <d 


OJ  jcvj 

^1  d 


K^  lfr\ 

od 


fO  {tr\ 
td 


CVJ 

0 

kj 

< 

n 

cy  Icy 

1  ■ 

^1  d 

xi|d 

1 

0 

<0 

1— t 

d. 

>...✓ 

1 

J3|  d 

^  1^ 

usi  as 

:■■-{ 
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Equations  (45)  through  (50)  are  satisfied  if  the  constants  B^,  C^, 
L^,  and  are  all  equal  to  zero.  This  represents  the  uniformly  com¬ 
pressed  circular  form  of  equilibrium  of  the  cylinder.  A  buckled  form  of 
equilibrium  is  possible  only  if  equations  (45)  through  (50)  yield  non¬ 
zero  solutions  for  the  constants.  This  requires  that  the  determinant  of 
the  coefficients  of  these  constants  be  equal  to  zero.  This  determinant  ^ 
is  shown  on  the  following  page. 
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The  lowest,  real,  negative  value  of  a  for  which  this  determinant  equals 

zero  is  proportional  to  the  critical  load  on  the  sandwich  cylinder;  this 

value  of  a  will  be  referred  to  hereafter  as  a  .  After  this  value  is 
—  cr 

determined,  the  value  of  the  critical  load  is  obtained  from  the  defini¬ 


tion  of  a  given  previously,  that  is. 


Et 


Numerical  Computations 

Since  the  determinant  in  equation  (5^)  contains  the  eigenvalue,  a, 
in  four  of  the  six  rows,  the  expansion  of  this  determinant  would  result 
in  a  fourth  order  algebraic  equation  in  a.  Only  the  lowest  of  the  four 
roots  of  this  equation,  the  root  that  corresponds  to  the  type  of  buckling 
illustrated  in  figure  5>  is  considered.  Actually,  for  very  short  cyl¬ 
inders  the  root  that  corresponds  to  a  face  wrinkling  type  of  buckling  may 
become  the  lowest;  a  satisfactory  analysis  of  this  type  of  buckling  re¬ 
quires  that  the  stiffnesses  of  the  individual  facings  be  taken'  into  ac¬ 
count.  Since  this  analysis  is  based  on  the  assumption  of  membrane  fac¬ 
ings,  no  attempt  is  made  to  de'termine  the  critical  loads  for  this  type  of 
buckling. 

A  literal  expansion  of  the  determinant  in  equation  (51)  was  found 
to  be  impractical.  It  was  decided  that  numerical  solutions  obtained  by 
the  use  of  machine  methods  would  be  the  most  practical  way  to  obtain  de¬ 
sign  data.'*  Originally,  the  sixth  order  determinant  as  expressed  in 

■••The  facilities  of  the  Numerical  Analysis  Laboratory  of  the  University  of 
Wisconsin  were  utilized  for  this  purpose. 
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equation  (51)  was  used  as  a  basie  for  numerical  solution.  However,  due  to 
inherent  difficulties  involving  the  subtraction  of  large  numbers  of  almost 
equal  magnitude,  it  was  found  that  sufficient  accuracy  could  not  be  ob¬ 
tained  from  this  approach.  The  sixth  order  deteralnant  was  then  reduced 
in  literal  form  to  the  fourth  order  determinant  shown  on  the  following 
page.  A  procedure  was  devised  from  which  the  desired  eigenvalue  can  be 
determined  from  this  determinant  by  trial  and  error.  For  a  given  set  of 
parameters  containing  the  dimensions  and  physical  constants  of  the  cyl¬ 
inder,  the  assumed  value  of  the  eigenvalue  is  varied  until  the  value  of 
the  determinant  becomes  equal  to  zero.  In  this  procedure  each  of  the 
twenty-four  terms  in  the  expansion  of  the  fourth  order  determinant  is  com¬ 
puted  independently  and  their  sum  is  accumulated  subsequently.  This  elim¬ 
inates  the  previous  difficulties  in  connection  with  losu  of  accuracy. 

Tile  numerical  results  obtained  are  shown  In  Tables  1  through  7« 

The  values  are  believed  to  be  correct  to  three  significant  figures.  More 
significant  figures  can  be  obtained  by  simply  increasing  the  number  of 
trials  used.  Curves  based  on  the  values  given  in  the  tables  are  shown  in 
figures  6,  T ,  and  8. 
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V. 


Discussion  of  Results 


The  curves  shown  In  figures  6,  J ,  and  8  are  sufficient  to  illus¬ 
trate  how  the  value  of  the  critical  load  varies  when  the  parameters 

representing  the  dimensions  and  physical  constants  of  the  sandwich  cyl- 

E  E 

inder  are  varied.  The  ratios,  — •  1*  and  =>  10,  were  chosen  since 


G, 


r0 


rz 


they  represent  average  values  of  these  quantities  for  honeycomb  cores, 
the  type  most  commonly  used  in  sandwich  construction.  Computations  were 
also  made  with  the  values  of  these  ratios  interchanged  since  either  situ¬ 
ation  is  possible  depending  upon  the  orientation  of  the  core  in  fabrica¬ 
tion.  The  numerical  results  shown  in  Tables  1  through  k  and  the  curves 
shown  in  figures  6  and  7  are  based  on  the  value  of  p  =■  1.  This  value 
represents  a  cylinder  that  has  a  weak  core;  that  is,  a  core  having  rela- 

The  numerical  results  shown  in 
Tables  5  through  7  and  the  curves  shown  in  figure  8  are  based  on  the 


tively  low  values  of  E  ,  Q_^,  and  G 

c  rw  rz 


value  of  3  =  1,000.  This  value  of  3  represents  a  cylinder  having  a 

E 


stiff  core.  Since  3  depends  not  only  upon  the  ratio,  but  also  upon 

—  E 


the  ratio,  the  practical  range  of  values  of  3  is  quite  wide. 


Most 


constructions  have  dimensions  and  physical  properties  such  that  they  lie 
in  the  range  between  3=1  and  3  =  1,000.  All  of  the  computations  are 
based  on  the  value  of  .5  for  the  Poisson's  ratio  of  the  facings.  Re¬ 
sults  were  obtained  for  only  two  different  ratios  of  the  radius  of  the 

b 


inner  facing  to  the  radius  of  the  outer  facing;  that  is, 


,95  and 


—  =  .98.  The  procedure  for  machine  computation  was  set  up  so  that  anj’  or 

Cl 

all  of  the  valvies  of  these  parameters  may  be  changed  without  any 
difficulty. 
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The  results  Indicate  that,  for  cylinders  having  relatively  weak 
cores.  It  Is  advantageous  to  have  the  higher  value  of  the  modulus  of  ri¬ 
gidity  of  the  core  In  the  tangential  direction  rather  than  In  the  longi¬ 
tudinal  direction.  The  effect  of  the  Interchange  of  these  values  Is  evi¬ 
dent  from  a  comparison  of  the  curves  shown  In  figures  6  and  7*  This  ef¬ 
fect  Is  negllgihle  in  the  case  of  cylinders  having  strong  cores  as  evi¬ 
denced  by  the  fact  that  the  values  given  In  Tables  5  6  are  essentially 

the  same.  The  fact  that  In  all  cases  the  overall  stiffness  properties  of 
the  core  become  relatively  more  important  as  the  length  of  the  cylinder 
decreases  may  also  be  noted.  As  the  length  of  the  cylinders  is  increased, 
the  values  obtained  for  the  critical  loads  approach  the  values  previously 
obtained  on  the  assumption  that  the  cylinder  length  is  infinite.  This  Is 
to  be  expected  because,  If  X  is  allowed  to  approach  zero,  the  character¬ 
istic  determinant  for  this  problem  reduces  to  the  determinant  obtained  on 
the  basis  of  the  assumption  of  infinite  length.*  As  the  cylinder  de¬ 
creases  In  length,  the  value  of  the  critical  load  approaches  the  value 
(l  -  — )  and  then  falls  off  rapidly  to  zero.  In  any  particular  case, 

if  is  set  equal  to  Infinity,  the  critical  load  becomes  equal  to 

G  Q  (l  -  — )  when  the  length  becomes  zero.  This  limiting  value  for  the 
rt?  a 

critical  load  is  characteristic  of  stability  analyses  of  sandwich  con¬ 
struction  ill  general,  if  is  assumed  to  be  infinite.  Examination  of 
the  curves  shown  in  figures  6,  7,  a-nd  8  shows  that  the  curves  represent¬ 
ing  the  cylinder  having  a  weak  core  reach  a  maximum  in  the  neighborhood 

*See  Ref.  2. 
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of  •£.  =>  1  while  those  representing  the  cylinder  having  a  strong  core  con- 
a _ 

tlnue  to  rise  until  the  cylinder  becomes  very  short.  Indications  are 
that  analyses  based  on  the  assumption  of  E^,  =  ®  are  sufficiently  accurate, 
in  most  ranges,  but  for  relatively  short  cylinders  having  weak  cores  such 
an  assumption  may  result  in  serious  error. 

Conclusions 

The  stability  analysis  of  sandwich  cylinders  subjected  to  uniform 
external  lateral  pressure  presented  here  is  believed  to  be  sufficiently 
accurate  for  use  in  design.  The  assumption  of  membrane  facings  limits 
the  range  of  applicability  of  the  results  somewhat  since  the  effect  of 
the  stiffnesses  of  the  individual  facings  may  be  appreciable  in  some  sand¬ 
wich  cylinders,  particularly  short  cylinders  having  relatively  thick  fac¬ 
ings.  However,  for  cylinders  of  usual  sandwich  construction  and  with  a 
length  equal  to  or  greater  than  the  radius,  this  theory  is  believed  to 
be  adequate.  The  analysis  of  the  problem  with  the  stiffnesses  of  the 
facings  taken  into  account  may  be  performed  without  a  great  amount  of 
additional  work,  and  such  an  analysis  is  planned  for  the  near  future. 
Additional  curves  giving  the  values  of  the  critical  loads  for  a  greater 
number  of  values  of  the  pareuneters  entering  the  problem  are  desirable 
for  design  purposes.  Such  curves  may  be  easily  prepared  in  the  future  if 
the  time  and  expense  involved  in  their  preparation  seem  justifiable.  The 
Forest  Products  Laboratory  is  at  present  planning  a  series  of  tests  for 
the  purpose  of  correlating  the  results  with  the  theoretical  results  ob¬ 
tained  here. 
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Table  1.  Values  of  Critical  Pressure  Expressed  in  Terras  of  d. 


cf 

X 

b 

a 

Icr  “ 
Et 

IT  a 

J 

=  .95 

p  =  1 

^  4 

□re 

Ec 

^  =  10 
°re 

X 

:  n 

:  -cf 

:  X 

:  n 

:  -rf 

0 

2 

0.00273 

2.62 

8 

0.0114 

0.25 

It 

.00292 

3.15 

II 

.0119  - 

.40 

tt 

.00373 

3.15 

9 

.0117 

.50 

tt 

.00496 

3.15 

10 

.0116 

.50 

3 

.00564 

3.15 

11 

.0115 

.60 

rt 

.00588 

3.15 

12 

.0115 

.80 

fi 

.00673 

3.15 

13 

.0115 

1.00 

fi 

.00826 

3.15 

14 

.0115 

.80 

4 

.00768 

3.50 

13 

.0115 

1.00 

fl 

.00808 

3.50 

14 

.0115 

1.20 

If 

.00869 

3.50 

15 

.0115 

1.40 

ff 

.00936 

4.00 

15 

.0114 

1.60 

It 

.0107 

4.00 

16 

.0114 

1.40 

5 

.00950 

4.50 

16 

.0113 
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.00993 

4.50 

17 

.0112 
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ft 
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.0104 

4.50 
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.0106 
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.0115 
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23 
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Table  8.  Values  of  Critical  Pressure  Iltprossed  in  Terras  of  d 


E 

-  “  .98  — £  =  10 

a  G 
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n 
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Table  Values  of  Critical 
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3ure  Expressed  In  Terms  of  eC . 
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Table  5»  Values  of  Critical  Pressure  Expressed  in  Terms  of  a 
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Table  6 


.  Valurja  of  Critical  Trassiu'c  lirin'cariel  in  n/’  <i. 
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Table  7 


Values  of  Critical  Pressure  Expressed  in  Terras  of  a 
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Figure  U.- -Differential  elements  of  deformed  facings. 
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Figure  6. —Critical  preesure  in  terms  of  versus  the  length  to  radius  ratio. 
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Figure  7. --Critical  presexire  in  terms  of  versus  the  length  to  radius  ratio. 


SECTION  VIII 

DESIGN  CURVES  FOR  THE  BUCKLING  OF  SANDWICH  CYLINDERS 

1 

OF  FINITE  LENGTH  UNDER  UNIFORM  EXTERNAL  LATERAL  PRESSURE- 


i ' 

By 


CHARLES  B.  NORRIS,  Engineer 
JOHN  J.  ZAHN,  Engineer 

Forest  Products  Laboratory,  —  Forest  Service 
U.  S.  Department  of  Agriculture 


This  report  contains  curves  and  formulas  for  the  calculation  of  the  critical 
external  pressure  of  sandwich  cylinders  of  finite  length.  The  facings  of  the 
sandwich  are  equal  and  isotropic  and  their  individual  stiffness  is  not  taken  in¬ 
to  account.  The  core  is  isotropic  or  orthotropic  having  natural  axes  in  the 
axial,  tangential,  ,  and  radial  directions  of  the  cylinder.  Curves  are  given  for 
isotropic  cores  and  for  orthotropic  cores  having  certain  relative  elastic  proper 
ties.  If  the  cores  are  very  rigid,  the  method  yields  results  that  are  substan¬ 
tially  those  of  von  Mises. 


Introduction 


This  report  presents  design  curves  for  the  critical  external  pressure  on  sand¬ 
wich  cylinders,  calculated  according  to  the  formulas  developed  in  Forest 


—This  progress  report  is  one  of  a  series  (ANC-23,  Item  57-3)  prepared  and 
distributed  by  the  Forest  Products  Laboratory  under  U.  S.  Navy,  Bureau 
of  Aeronautics  Order  Nos.  NAer  01898  and  U.  S.  Air  Fotce  Contract  No. 
DO  33(616)58-1.  Results  here  reported  are  preliminary  and  may  be  re¬ 
vised  as  additional  data  become  available. 

2 

"Maintained  at  Madison,  Wis.  ,  in  cooperation  with  the  University  of  Wiscon¬ 
sin. 
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Products  Laboratory  report  1844-B  (9).  ““  The  sandwich  cylinders  have  iso¬ 
tropic  facings  and  orthotropic  or  isotropic  cores.  The  natural  axes  of  the 
orthotropic  cores  are  axial,  tangential,  and  radial.  These  formulas  reduce 
substantially  to  those  developed  by  von  Mises  (7,  10,  14)  when  the  core  is 

very  rigid  so  that  the  stiffness  of  the  spaced  facings  of  the  sandwich  (no  re¬ 
duction  of  the  stiffness  due  to  shear  strains  in  the  core)  can  be  used  as  sug¬ 
gested  in  3.1.5  of  ANC-23  (11). 

A  great  deal  of  investigative  work  has  been  done  on  isotropic  cylindrical  shells 
subjected  to  external  pressure  since  report  1844-B  giving  theoretical  analysis 
of  sandwich  cylinders  was  written.  It  has  been  found  that  experiment  some¬ 
times  yields  critical  loads  that  arc  less  than  those  predicted  by  von  Mises' 
theory  (13).  This  has  been  attributed  to  two  causes.  First,  the  experimental 
cylinders  contained  imperfections  that  lowered  the  critical  load  (1,  2,  3,  8, 

12).  Second,  lower  energy  levels  are  associated  with  post-buckling  configura¬ 
tions  of  the  cylinder  than  with  those  just  at  buckling,  the  former  being  reached 
without  the  necessity  of  pAssing  through  (snap-through)  the  latter  or  the  energy 
necessary  for  passing  through  the  latter  being  supplied  by  vibration  or  shocks 
(4,  5,  6,  8). 

The  curves  published  in  this  report  do  not  consider  snap-through  buckling  or 
cylinders  with  imperfections.  Sandwich  cylinders,  howev'er,  are  much  more 
perfect  than  their  solid  counterparts  because  they  are  thicker  and  the  effect 
of  an  imperfection  is  in  proportion  to  the  ratio  of  its  amplitude  to  the  thick¬ 
ness  of  the  cylindrical  shell.  Also,  the  curves  neglect  the  stiffnesses  of  the 
individual  facings.  These  stiffnesses  add  a  considerable  amount  to  the  criti¬ 
cal  loads  when  the  cylinders  are  short,  and  it  is  for  short  cylinders  that  "snap- 
through'!  is  likely  to  occur  (6).  Thus  it  seems  that  these  curves  are  useful  in 
the  design  of  sandwich  cylinders. 

/  * 

Development  of  Formula  from  which 
Design  Curves  Were  Calculated 

The  critical  external  pre.<5sure  is  found  by  placing  the  determinant  on  page  23 
of  report  No.  1844-B  (9)  equal  to  zero  and  solving  for  a.  This  determinant 
can  be  simplified  if  the  transverse  modulus  of  elasticity  of  the  core  (E^)  is 
assumed  to  be  infinite.  For  most  core  materials  except  possibly  for  low 
density  foams  is  sufficiently  large  so  that  this  assumption  yields  values  of 

the  critical  pressure  that  are  only  very  slightly  too  great. 

3  ^ 

Underlined  numbers  in  parentheses  refer  to  the  references. 
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Before  is  allowed  to  approach  infinity,  the  first  and  fourth  columns  of  the 

Gr2 

determinant  are  multiplied  by  - and  the  third  column  is  divided  by  E^. 

^c 

Then  when  E^,  approaches  infinity,  the  expressions  in  rows  3,  4,  5,  and  6  in 
column  3  approach  zero. 


The  expressions  in  each  row,  excepting  the  first,  are  replaced  by  new  ex¬ 
pressions,  as  indicated  by  the  following  formulas  in  which  R  represents  the 
expression  in  the  row  designated  by  its  subscript  and  in  some  column.  The 
primed  values  are  the  new  ones  to  be  substituted  for  the  old.  These  substi¬ 
tutions  can  be  made  without  changing  the  value  of  a  because  of  the  well-known 
properties  of  determinants  and  because  the  determinant  is  equated  to  zero. 


^2 

=  ^2  ~2  +  ^ 

1 

R3' 

=  R3  + 

^4’ 

=  R4  +  R^’ 

^5’ 

=  2R5  +  (n^ 

+  3X^)  R^’ 

=  2Rb  +  (n^ 

+  3  ^  K^)  R2 

These  substitutions  cause  the  expressions  in  column  3  and  rows  2,  3,  4,  5, 
and  6  and  those  in  column  6  and  rows  3,  4,  5,  and  6  to  become  zero,  and  the 
determinant  is  readily  reduced  (by  minors)  from  a  6-by-6  to  a  4-by-4  deter¬ 
minant.  This  determinant  is  simplified  slightly  by  replacing  the  second  row 
of  expressions  by  the  second  row  minus  the  first  row. 

After  the  determinant  was  written  in  this  form,  a  change  in  parameters  was 
made  using  the  following  nomenclature: 

R  -  -  mean  radius  of  the  sandwich  cylinder 

h  --  thickness  of  the  sandwich 


c  --  thickness  of  the  core  of  the  sandwich 

n  --  number  of  half  waves  in  the  circumference  of  the  cylinder 
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V  = 


E 


t 


(1  -  h 


h  q(l  -  h 

K  _  a-- - g  - 


where  G^^  and  Gj^q  are  the  moduli  of  rigidity  of  the  core  associated  with  the 
radial  and  axial  directions  and  with  the  radial  and  tangential  directions;  E, 
and  t  are  the  modulus  of  elasticity,  Poisson's  ratio,  and  thickness  of  the  fac¬ 
ings;  and  q  is  the  external  critical  pressure  on  the  cylinder. 


The  radii  a  and  b  were  eliminated  by  the  following  equations  obtained  from  the 
geor  etry  of  the  cylinder;  '  i 


a  =  R  +  -  (h  +  c) 


b  =  R  -  1  (h  +  c) 


and  the  following  substitutions  made: 


h  +  c 


The  expressions  in  the  final  determinant  are: 


Row  1,  column  1 


n  -  1  $  +  4) 


n^  $  -  4) 


3>  +  4) 


/iCy  + 

3n^  \f  /  \  ^ 


Row  2,  column  1 


1  I  8^  ^ 

n^$-4) 


Row  3,  column  1 


-5i_  L2  +  3,2  4.2v  [i 

\l/  \  ^  /  ^  \i/  _3 


iT^  /r\^  /$  +  4>V 


n^  \l 
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Row  4,  column  1 


Row  1,  column  2 


Row  2,  column  2 


2(f> 

3>  +  (j) 


Row  3,  column  2 


+  (n^  -  1) 


Row  4,  column  2 


2  v  4>  +  4>  ^ 


$2-4,2 


2  X  ^ 


K  +  3  IT' 


R\ 


^  -  <{>  ,  2  /R\^  (^  -  «i>)^  X  fn2 

4>+  4>  \i/  $($  +  4>)  ^ 


2  ,,  ($+  4>)^  +  4.^ 

4$  -  4,^^ 

Row  1,  column  3 
-^+ 

$+4)  3r$+4)  ^$+4) 

Row  2,  column  3 


2  X  ^,2 


K|  n'^  +  3tt 


RY  ^  - 


-  -i-  V  $  -  4» 

3r$+4>  °®$+4) 
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Row  3,  column  3 


-  2  +  _ii_  +  3^2 

$  +  4.  ^  +  3Tr 


.  3.^  (f)  (i 

(fj  (^/l  ■ 

'  _ 


2  /  v2l 

'  |±±i )  .  jL  V 
\  ^  J  3r  $  + 


$  -  4) 
$  +  4> 


Row  4,  column  3 


-'Hf/ Ml  ■».  fvf 


_4  V 
3r  $  +  4) 


Row  1,  column  4 

Row  2,  column  4 

.  llL.  4> 


3  \  i  ;  $ 

Row  3,  column  4 

fyf  n^  f-  -  ^  kI  (^' 

Vv  V  $  /  "  3  16$  K  ‘^iTy 


-^(ff  Mr 


,  2  /„x2  2  2  r 

2tt  /  R\  $  +  4)  - 

'  ‘  3  It/  — ? —  " 

J  ^  \^/  $2 


Row  4,  column  4 


ilf  M)'  [t  -  r  -  ¥  (f)'  P 


Mf 


Thia  det.,™„.„t  was  placed  ,q„al  ,o  zero  and  values  of  K  determined  by 


electronic  computer  for  various  values  of  — 


R’  T' 
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The  limits  of  these  curves  for  very  long  and  very  short  cylinders  may  be 
found  from  Forest  Products  Laboratory  Reports  1844-A  and  1844-B  (9).  For 
infinitely  long  cylinders  having  membrane  facings  and  for  v/hich  the  modulus 
of  elasticity  of  the  core,  in  the  radial  direction  is  infinite,  equation  (72)  of 
Report  -1'844^A  becomesr^ 


K  = - 

■‘2  ($2  +  4,  2) 


3  ($-<)))  ^ 

-  4,^  _ Et 

16<|)  Grq  (1  -  p^)  h 


For  very  short  cylinders  having  membrane  facings  and  for  which  the  modulus 
of  elasticity  of  the  core  in  the  radial  direction  is  infinite,  the  equation  on  page 
28  of  report  1844-B  becomes: 

“  *  +  9  Et 

From  this  equation  the  critical  hoop  compression  per  unit  length  of  cylinder 
is  found  to  be: 

~  z"  ^R0 


which  is  the  usual  limit  imposed  on  the  edge  compression  of  sandwich  con¬ 
structions  with  membrane  facings  by  the  shear  instability  of  the  core. 


Description  of  Design  Curves 


The  design  curves  (figs.  1  to  6)  apply  to  sandwich  cylinders  having  equal  iso¬ 
tropic  facings  and  isotropic  cores  or  orthotropic  cores  having  thin  natural 
axes  parallel  to  the  axial,  tangential,  and  radial  directions  of  the  cylinders. 

I 

They  are  plots  of  K  against  —  where  the  critical  pressure  is  given  by 

XV 


Et 

(  1  -  h 


K 


and  E,  p,  and  t  are  the  modulus  of  elasticity,  Poisson's  ratio,  and  thickness 
of  the  facings;  h  the  thickness  of  the  sandwich;  I  and  R  the  length  and  mean 
radius  of  the  cylinder.  Six  curve  sheets  are  presented,  two  for  each  of  three 

values  (10,  50,  and  100)  of  — .  One  of  each  pair  of  curve  sheets  applies  to 


( 
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sandwich  having  very  thin  facings  for  which  the  ratio  of  the  thickness  of  the 
core  to  the  thickness  of  the  sandwich  (c’/h)  is  substantially  unity.  The  other 
applies  to  sandwich  for  which  this  ratio  is  0.7.  Each*curve  applies  to  sand¬ 
wich  having  a  particular  value  of: 


where  is  the  modulus  of  rigidity  of  the  core  associated  with  the  radial  and 

tangential  directions  of  the  cylinder.  It  was  found  that  the  modulus  of  rigidity 
of  the  core  associated  with  the  radial  and  axial  directions  (Gpj^^)  has  very  little 
influence  on  the  critical  pressure.  It  does  not  enter  the  formulas  for  the  criti¬ 
cal  pressure  of  very  long  or  very  short  cylinders.  Calculations  were  made  for 
cylinders  having  three  values  of 


r 


G 

R 


Rz 

Re 


These  values  were  2.5,  1.0,  and  0.4  to  agree  with  the  values  appropriate  for 
some  honeycomb  and  isotropic  core  materials.  Many  of  the  curves  are  not 
affected  by  the  use  of  these  values.  Some  of  them  are  affected  slightly.  This 
is  shown  in  the  figures  by  the  use  of  three  adjacent  curves;  the  greatest,  in¬ 
termediate,  and  least  values  of  K  are  associated  with  the  greatest,  intermedi¬ 
ate,  and  least  values  of  r. 
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Figure  1.  --Values  of  K  for 
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